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This  report  is  a  user's  manual  for  ALBERT,  a  package  of  four 
Caltech  computer  programs  which  calculate  analytic  expressions  for 
general  relativistic  curvature  tensors  and  equations  of  motion. 

The  input  data  for  ALBERT  are  analytic  forme  for  the  c  ova  riant  and 
contravariant  components  01  the  metric  teneor,  gafe  and  g  ,  and  fox 
the  mixed  components  of  the  stress -energy  tensor,  .  The  output, 
variable  at  the  discretion  of  the  user,  Includes  analytic  exp res- 
8 ions  for  the  equations  of  motion,  ,b  *  0;  for  the  Chrlstoffel 
symbols,  I*  and  r#bc;  for  the  Rlemann  tensor,  R#bcd;  for  Ricci 
tensor,  R*b;  for  the  Einstein  tensor,  C*b;  and  for  the  Ricci  scalar, 

R.  An  option  is  included  for  expanding  the  answers  in  a  power  series 
in  some  parameter,  c,  which  appears  in  the  input  tensors;  and  for 
keeping  only  terms  of  zero,  first,  and  second  order.  ALBERT  is  written 
in  the  IBM  language  FORMAC,  which  is  quite  similar  to  FORTRAN;  but 
which,  unlike  FORTRAN,  provides  for  algebraic  manipulation  of  symbols 
and  for  analytic  differentiation  of  functions.  ALBERT  should  be 
usable,  without  modification,  on  any  computer  which  has  a  FORMAC 
costlier.  A  complete  list  ag  of  the  programs  in  ALBERT  is  Included 
'.in  this  report. 


in  part  by  the  Rational  Science  Foundation  [OP-5391] 
the  Office  of  Naval  Research  [lonr-220(A7) ). 


Supported 


4. 


'Alfred  P.  Sloan  Foundation  Research  Fellow. 


Jtp'&'cr 

TMjHWP' - --- 

THIS  TECHNICAL  REPORT  IS  ONE  OF  THE  GRAPE  A1V  SER\tfii(g  't' 

PpSCIENTITK  UEMORANPA,  TECHNICAL  REPORTS  »  MI  SCE  IE 


WRITINGS. 


August  196? 


Blank  Page 


TABLE  OP  CONTENTS 


rage 

I.  Introduction  and  Motivation .  1 

*11.  Brief  Description  of  ALBERT  . .  2 

III.  How  to  Use  the  Progress  in  ALBERT . . . .  .  4 

*A.  Brief  Description  of  the  FORMAC  Language .  4 

■ft  _ 

B.  EINSTEIN,  the  Program  for  Calculating  Ricci  and 

Einstein  Tensors  . .  9 

1.  Input  for  EINSTEIN  . .  9 

2.  Complete  Listing  of  EINSTEIN .  13 

3.  Output  for  EINSTEIN  . .  21 

it.  Sample  Problems  for  EINSTEIN .  25 

C.  RIEMANN,  the  Program  for  Calculating  Riemann  Tensors  ...  & 

1.  Input  for  RIEMANN . 34 

2.  Complete  Listing  of  RIEMANN .  35 

3.  Output  for  REMANN  . .  40 

4.  St;  pie  Problems  for  RIEMANN .  41 

D.  MOTION,  the  Program  for  Calculating  the  livergence 

of  the  Stress-Energy  Tensor .  42 

1.  Input  for  MOTION .  42 

2.  Complete  Listing  of  MOTION . 43 

3.  Output  tar  MOTION .  49 

4.  Sample  Problems  for  MOTION  . . 51 

E.  SUBSTITUTE,  the  Program  for  Making  Changes  of  Variables  .  .  St: 

1.  Input  for  SUBSTITUTE .  56 


*The  reader  interested  in  setting  up,  as  rapidly  as  possible,  a  working  program 
for  calculating  Einstein  and  Ricci  tensors  need  read  only  these  sections. 


Page 


2.  Complete  Hating  of  SUBSTITUTE .  59 

3.  Output  of  SUBSTITUTE  . .  62 

4.  Sample  Problems  for  SUBSTITUTE .  64 

IV.  A  Brief  Guide  to  FOHMAC  Programing  . .  69 

A.  The  FORMAC  System . 69 

B.  FORMAC  Executable  Statements  . .  70 

C.  Ways  to  Avoid  Exhausting  the  Computer's  Memory  ......  75 

V.  References .  76 


il 


I.  INTRODUCTION  AND  MOTIVATION 


The  initial  step  In  ma ay  general-relativity  research  projects  is  to 
select  a  coordinate  Bystem  and  corresponding  analytic  forms  for  the  metric 
tensor,  g#b,  and  the  stress-energy  tensor,  W^;  and  then  to  confute  the  equa¬ 
tions  of  motion,  W); :  ,b  »  0,  and  the  Einstein  tensor,  0^,  for  use  in  the  field 
equations  0  b  »  8*  W  At  later  stages,  one  also  has  occasional  need  for  the 
Riemann  and  Ricci  tensors,  Rabcd  and  Rab,  and  for  the  Ricci  scalar,  R. 

So  long  as  one  deals  with  simple  systems  which  possess  a  high  degree  of 
symaetry,  the  computation  hy  hand  of  Wa^  b,  0^,  R*b,  end  R  ie  aot  too 

tedious.  However,  in  recent  years,  the  systems  under  study  have  become  more 
and  more  conplex;  and  the  task  of  computing  the  necessary  equations  of  motion 
and  curvature  tensors  has  became  more  and  more  formidable.  Recent  computations 
of  curvature  tensors  have  sometimes  required  many  weeks  of  tedious  labor,  and 
the  possibilities  of  mistakes  of  minuB  signs  and  factors  2  in  the  final  results 
have  been  very  disturbing. 

In  response  to  this  situation,  we  have  developed  the  computer  programs 
described  here  for  calculating  curvature  tensors  and  equations  of  motion. 

These  programs  have  been  used  extensively  over  the  last  seven  months  in  rela¬ 
tivistic  astrophysics  research  at  Caltech.  (See,  e.g.,  the  analysis  of  the 
effects  of  rotation  on  the  structure  of  general  relativistic  stellar  models  by 
Hartle  (1967);  also,  the  analysis  of  nonradlal  pulsations  of  relativistic 
stellar  models  by  Thorne  and  Caapolattaro  (1967).) 

After  these  programs  were  completed  and  in  use  at  Caltech,  we  learned  that 
Fletcher  (1965,  1966)  and  Clemens  and  Matzner  (1967)  have  developed  independently 
computer  programs  for  doing  the  same  Job.  The  three  different  systems  are 
described  Tie  fly  in  Fletcher,  Clemens,  Matzner,  Thome,  and  Zimmerman  (1967). 
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ALBERT  la  a  package  consisting  of  four  independent  computer  progress  — 
EINSTEIN,  RIEMANN,  MOTION,  and  SUBSTITUTE  —  which  perform  four  independent 
tasks. 

EINSTEIN  is  a  program  which  takes  a  metric  tensor,  g^,  and  calculates 
for  it  the  Christoff  el  symbols,  and  r^0,  the  mixed  components  of  the 
Ricci  tensor,  R#  ,  the  Ricci  scalar,  R,  and  the  mixed  components  of  the 
Einstein  tencor,  G#b  -  rJ*  -  J  R  6gb.  The  input  to  EINSTEIN  are  functional 
forms  far  gab  and  gab;  and  the  output  are  printed  expressions  for  all  the 
above  quantities . 

RIEMANN  la  a  program  for  calculating  the  covarlant  components,  R  ... 
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of  the  Rlemann  curvature  tensor  far  any  given  metric  tensor,  g^.  The  input 
to  RIEMANN  are  functional  farms  for  g^b  and  gabj  and  the  output  are  printed 
expressions  for  V^,  r^0,  and  R^. 

MOTION  la  a  program  far  calculating  the  divergence,  V#b.b,  of  a 
s  tress -energy  tensor,  W#  •  (We  denote  the  stress -energy  tensor  by  W  b  rather 
than  T#b  because  the  symbol  T  is  reserved  for  the  tin*  coordinate.)  The  input 
to  MOTION  are  functional  forms  for  g>b,  gab,  and  W#b{  and  the  output  are 
printed  expressions  for  r  r  . c,  and  W  b  .  . 
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SJBSTITOTE  is  a  program  for  making  changes  of  variables  in  long,  ccqpU- 
cated  expressions .  The  input  to  SUBSTITUTE  are  the  changes  of  variables 
desired,  e.g. 

X  »  -  in  (l-2*^X); 

and  the  expressions  in  which  those  changes  are  to  be  made,  e.g. 


e-x  f  _L  +  x. 

lx2  WV 
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The  output  ere  printed  expressions  for  the  desired  quantities  in  terns  of  the 
new  variables,  e.g. 


The  programs  EINSTEIN,  KEEMANN  and  MOTION  are  severely  limited  by  core 
storage  problems:  They  are  Incapable  of  calculating  curvature  tensors  and 
equations  of  motion  for  metrics  and  stress-energy  tensors  which  are  too 
complicated .  (Examples  of  input  which  will  and  will  not  coqpute  are  given  in 
44  III.B.4,  HI.C.4,  III.D.4,  and  UI.E.4.)  However,  complicated  problems  can 
be  solved  by  putting  into  EINSTEIN,  R3334ANN,  and  MOTION  simplified  metric  forms, 
e.g. 

g^  -  a(x,y,t) 

rather  than 

-  ev(x)  (1.  +  A^XjT)  P^cos  Y)  +  Ag(X,T)  Pg(coe  Y)}} 

and  by  then  using  SUBSTITUTE  to  make  a  change  of  variables  in  the  resultant 
curvature  tensors  and  equations  of  motion,  e.g. 

A(X,Y,T)  -  eVM  (1.  +  AX(X,T)  P^coe  Y)  +  Ag(X,T)  Pg(con  Y)}. 

An  option  is  Included  in  EINSTEIN,  RHMANN,  MOTION,  and  SUBSTITUTE  for 
expanding  the  output  in  powers  of  a  small  parameter  and  retaining  only  terms 
of  zero,  first,  and  second  order. 

EINSTEIN,  RIB4A1IN,  MOTION,  and  SUBSTITUTE  are  all  written  in  the  IBM 
language  FOHMAC.  P0H4AC  is  an  extension  of  FORTRAN  IV,  hut  unlike  FORTRAN  IV, 
FQRMAC  can  currently  be  used  only  on  IBM  7090  end  7094  computers .  Hopefully 
an  improved  version  of  FORMAC  will  be  available  for  the  IBM  360  system  some¬ 
time  in  196S. 
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in.  RCW  TO  USX  THE  PROGRAMS  XH  AUERS 
A.  Brief  Description  of  the  FORMAC  language 

In  this  section  ve  present  only  enough  information  on  FORMAC  to  eneble 
the  reader  to  use  the  programs  in  ALBERT.  A  mere  detailed  description  of 
FORMAC,  intended  to  help  the  reader  vho  vlshes  to  write  slarle  FORMAC  programs 
far  himself  or  to  modify  ALBERT  to  suit  his  own  individual  needs,  will  he 
found  in  5  IV  of  this  report.  For  a  complete  description  of  the  FORMAC  system 
and  for  guidance  in  writing  complicated  FORMAC  programs,  the  reader  Is  re* 
f erred  to  the  standard  DM  FORMAC  i  (IBM  1964). 

The  FORMAC  system  le  capable  of  performing  formal  algebraic  manipulations, 
and  differentiations  of  algebraic  expressions .  It  can  also  evaluate  numerically 
the  algebraic  expressions  vhich  it  creates,  but  no  such  numerical  evaluations 
are  used  in  the  programs  described  in  this  report. 

FORMAC  uses  notation  vhich  ia  almost  identical  to  that  of  FORTRAN  IV: 

1.  All  algebraic  quantities  —  constants,  parameters,  spacetime  coordi¬ 
nates,  functions,  algebraic  expressions  —  are  represented  by  'words" 
vhich  are  constructed  from  6  or  fever  Latin  letters  and  digits.  For 
example,  in  ALBERT  ve  use  X, Y,Z,  and  T  to  represent  the  4  spacetime 
coordinates;  and  ve  use  BP  to  represent  a  small  parameter  in  terms 
of  which  power  series  expansions  are  made.  As  further  examples,  in 
the  sample  problem  given  in  {  IIX.B.4,  ANU  represents  the  function 
v(X)  and  HI  represents  the  function  H^(X,T).  As  a  final  example, 
with  the  above  definitions  of  X,  Y,  Z,  T,  BP,  ANU,  and  HI,  the  state¬ 
ment 

LET  PAM3  -  T  *  ANU  +  Hl/l, 

tells  the  cosputer  to  use  the  word  PAM3  to  represent  the  algebraic 
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expression  T  v(x)  +  Hj(X,T)/Y.  (When  defining  en  algebraic  expres¬ 
sion,  one  suit  always  use  the  word  LET  in  front  of  the  definition# 
Otherwise,  the  counter  will  atteapt  to  evaluate  the  expression 
nunerlcally,  in  the  usual  FORTRAN  manner.) 

2,  Algebraic  quantities  can  be  subscripted  In  the  same  Banner  as  is  done 
In  FORTRAN#  Far  example,  we  denote  the  covariant  components  of  the 
metric  tensor,  g.  .,  by  the  words  Gl(l,J);  0l(l,3)  Is  g^,  Ol(2,2) 

ie  $g2*  etc# 

3.  The  computer  distinguishes  constants  and  poraaaters  from  functions 
and  from  algebraic  expressions  by  means  of  "ATfSfflC"  and  "DEPEND" 
statements,  which  appear  near  the  beginning  of  the  program.  (These 
statements  are  similar  to  the  E3QERKAL  statement  of  FORTRAN  IV.) 

a.  The  ATOMIC  statement  identifies  all  algebraic  quantities  ("atomic 
variables")  which  do  not  represent  internally  constructed  alge¬ 
braic  expressions.  All  nonnumerlcal  constants,  parameters, 
spacetime  coordinates,  and  functions  must  appear  in  the  ATOMIC 
statement.  For  example,  for  a  program  involving  only  the  quanti¬ 
ties  described  in  1.  above,  all  quantities  except  PAM3  appear  in 
the  ATOMIC  statement,  and  the  atomic  statement  takes  the  form: 

ATOMIC  X,  Y,  Z,  T,  EP,  ANU,  HI. 

b.  The  DEPEND  statement  identifies  all  functions  which  are  not 
algebraic  expressions  —  l.e.,  which  are  atomic  and  it  states 
the  functional  dependences  of  the  functions.  For  example,  one 
identifies  v(x)  and  H1(X,T)  as  atomic  functions  by  the  above 
ATOMIC  statement  plus  the  DEPEND  statement 

DEPEND  (ANU/X),  (Hl/X,T). 

Note  that  PAMZ  as  defined  above  is  not  an  atomic  function;  it 
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is  an  algebraic  expression  made  up  of  atomic  functions  and  other 
atonic  variables,  and  it  therefore  does  not  appear  in  either  an 
ATOMIC  statement  or  a  DEPEND  statement. 

4.  In  constructing  algebraic  expressions  one  can  perform  the  elementary 
opersti'wjs  of  addition,  subtraction,  multiplication,  division,  and 
exponentiation  (l.e.,  raising  to  a  power).  The  notation  far  these 
operations  is  Identical  to  that  of  FORTRAN  IV: 

+  denotes  addition 

-  denotes  subtraction 

*  denotes  multiplication 

/  denotes  division 

**  denotes  exponentiation. 

As  in  FORTRAN  IV,  the  order  in  which  these  operations  are  performed 
is  **,  followed  by  *  and  /,  followed  by  +  and  -.  For  example,  con¬ 
sider  the  following  FOHMAC  statements  and  their  translation  to 
ordinary  algebraic  language: 


FORMAC  LANGUAGE 

CONVENTIONAL  NOTATION 

1ST  FAM3  -  T»ANU»Hl/l 

this  sets  PAM3  equal  to 

Tv^ 

1ET  SS  -  PAM3**2.+3.*(X-Y)**T 

this  sets  SS  equal  to 

In  constructing  algebraic  expressions  one  can  use  numbers  which  are 
floating  point  (l.e.  have  decimals  like  2.  and  3.  above)  or  which 
are  Integers  (e.g.  2  and  3).  When  doing  floating  point  manipulations 
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the  computer  generally  outputs  exponential  notation  --  e.g.  "2.5E-1" 
(2.5X10-1)  instead  of  l/4.  When  doing  integer  algebraic  manipula¬ 
tions,  FORMAC  uses  2**(-2)  instead  of  2.5E-1.  When  integer  (i.e. 
rational)  expressions  are  mixed  together  with  floating  point  expres¬ 
sions  by  the  user,  the  computer  converts  them  all  to  floating  point 
before  manipulating. 

5.  The  computer  recognizes  and  knows  how  to  differentiate  the  following 
elementary  functions: 


FORMAC  LANGUAGE 

CONVENTIONAL  NOTATION 

FMCEXP(ANU) 

v 

exponential:  e 

FMCSIN(2.«T) 

sine:  sin(2T) 

IMCC)j6(Hl**2) 

cosine:  cos^  ) 

fmcl#Sq(y) 

natural  logarithm:  log(Y) 

FMCATN(EMCEXP{  ANU) ) 

arctangent:  tan  1(ev) 

FMCtfTN(3*Z+7) 

hyperbolic  tangent:  tanh(3Z+7) 

These  functions  can  be  used  freely  in  algebraic  expressions. 

6.  FORMAC  is  capable  of  differentiating  any  algebraic  expression  con¬ 
structed  according  to  the  above  ruleB.  If  A  is  a  function  of  X,  Y, 
T,  but  not  of  Z 

DEFEND  (A/X,Y,T), 

then  various  derivatives  of  A  are  denoted  as  follows 


FORMAC  LANGUAGE 

CONVENTIONAL  NOTATION 

fmcdif(a,(x,  l)) 

dA/dX 

IMCDIF(A,(X,2)) 

d2A/dX2 

IMCDIF(A, (x,i) , (Y,l) ) 

s2A/ax9Y 

FMCDIF(A,(X,2),(T,1)) 

53a/^x23t 
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In  constructing  algebraic  expressions,  one  can  use  FORMAC  differentia¬ 
tion  freely,  and  the  requested  differentiation  will  always  be  performed 
to  completion.  For  example,  consider  the  following  sequence  of 
connands: 


FORMAC  CCMtANO 

RESULT  OF  C0M4AND 

AStklC  X,  Y,  Z,  T,  EP,  ANU,  HI 

Defines  the  atomic  variables  for  the 
program. 

DEFEND  (ANU/x),(H1/X,T) 

Defines  the  atomic  functions  and  their 
dependences . 

LET  FAM3  -  T*ANU*Hl/Y 

Bets  PAMS  equal  to  TvH^/Y. 

LET  Q1  -  FMCDIF(PAM3,T,l) 

Performs  the  required  differentiation 
on  PAM3,  yielding 

Q1  -  vHj/y  +  (Tv/YjbKj/dT. 

LET  CC  -  FMCDIF(FMCSIN(Hl)+T*X, 
(T,1),(X,1)) 

Performs  the  required  differentiation 
yielding 

CC  «  (cos  ^Hd^/Sr&X) 

-  (sin  HjKdHj/dsKdHj/dX) 

+  1 

I£T  CD  «  FMCDIF(ANU,Y) 

Performs  the  required  differentiation 
yielding 

CD  -  0. 

There  are  connands  in  the  FORMAC  language  which  remove  parentheses  from 
algebraic  expressions,  which  factor  terms  out  thereby  inserting  parentheses, 
which  truncate  power  series,  and  which  perform  other  useful  manipulations. 
However,  none  of  these  connands  are  necessary  to  the  use  of  the  programs  in 
ALBERT,  bo  we  shall  delay  describing  them  until  §  IV. B. 
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1.  Input  for  EINSTEIN 


EINSTEIN  Is  the  program  vhich  calculates  Ricci  and  Einstein  tensors. 

In  order  to  use  EINSTEIN,  one  must  supply  the  foilmrinp;- 

1.  An  AT)^tIC  statement  identifying  all  atomic  variables  for  the 
problem.  The  ATOMIC  statement  takes  the  form 

ATOMIC  X,  Y,  Z,  T,  EP,  .  .  .  . 

The  user  inserts  after  EP  all  functions  and  parameters  which 
appear  in  his  metric  tensor. 

2.  A  DEPEND  statement  identifying  all  atomic  functions  which  appear 
in  the  metric  tensor  and  specifying  their  functional  dependences. 

For  the  form  of  the  DEPfMD  statement,  see  page  5. 

3.  A  statement  declaring  the  value  of  the  integer  EFTERM.  EPTEHM  is 
used  to  tell  the  computer  whether  or  not  to  expand  the  curvature 
tensors  to  second  order  in  the  parameter  EP  and  truncate  thereafter 
(i.e.,  keep  only  terms  in  EP**0,  EP**1,  and  EP**2). 

If  EFTERM  -  0  terms  in  all  powers  of  EP  will  be  retained. 

If  EPTEHM  /  0  only  terms  of  order  0, 1,  2  in  EP  will  be  retained. 

4.  "IE T"  statements  defining  all  nonzero  covariant  components  of  the 
metric  tensor,  g^,  in  terms  of  atomic  variables,  g^  is  denoted 

in  ALBERT  by  <  i(l,j).  The  metric  components  1,  2,  3,  4  refer  to  the 
spacetime  coordinates  X,  Y,  Z,  T  respectively. 

5.  "LEST"  statements  defining  all  nonzero  contravariant  components  of 
the  metric  tensor,  g*^,  in  terms  of  atomic  variables,  g*^  is 
denoted  in  ALBERT  by  G2(l,j).  The  matrix  G2(I,J)  must  be,  of  course, 
the  algebraic  inverse  of  Gl(l,j)  -  at  least  up  to  second  order  in 
the  expansion  parameter  EP. 
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These  5  statements  are  Inserted  by  the  user  directly  into  the  body  of  the 
program  EINSTEIN  at  the  points  dlllneated  by  eoament  cards.  (See  listing  In 
the  next  section). 

For  exaagile  (another  example  Is  given  In  }  HI.B.4),  to  calculate  Wed 
and  Einstein  tensors  for  the  standard,  static,  spherically  symmetric  metric 

d«2  *  -e^r^dr2  -  r2(d02  +  sln20  dq>2)  +  ev^dt2  (l) 

one  must: 

1.  Mentally  replace  r  by  X,  0  by  Y,  <p  by  Z,  t  by  T,  v  by  ANU,  and  A  by 
AIM) 

2.  Supply  the  following  cards  at  the  required  points  In  EINSTEIN 

AlfiMIC  X,  Y,  Z,  T,  EP,  ANU,  ALAM 
DEPEND  (ANU,ALAH/X) 

EPTEHM  -  0 

LET  01(1,1)  -  -FMCEXP(ALAM) 

LET  01(2,2)  -  -X**2 

LET  01(3,3)  -  -(X*IMCSIN(Y))*»2 

LET  01(1,4)  -  IMCEXP(ANU) 

LET  02(1,1)  -  -FMCEXP(-ALAM) 

IET  02(2,2)  -  -l/x**2 

LET  02(3,3)  -  -l/(X*IMCSIN(Y))**2 

LEI  02(4,4)  -  PMCEXP(-ANU) 

As  in  FORTRAN  all  cards  supplied  by  the  user  Bust  have  their  first 
entry  In  or  after  column  7  and  their  last  entry  in  or  before 
72.  Long  statements  can  be  continued  over  as  many  as  nineteen 
continuation  cards.  Any  card  with  a  nonzero  entry  In  colunu  6 
considered  a  continuation  card. 
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3.  Run  the  program  on  an  IBM-7090  or  7094  computer  using  a  IBSYS-FORMAC 
operational  tape  (obtainable  from  IBM),  which  contains  the  FORMAC 
compiler . 

The  user  is  free  to  give  the  functions  and  variables  entering  into  hie  metric 
tensor  any  names  constructed  from  6  or  fewer  characters,  except  the  following: 


NAMES  WHICH  MUST  NOT  BE  USED  AS 
VARIABU3  NAMES  IN  ALBERT 


!•  X,  Y,  Z,  T  except  for  use  as  spacetime  coordinates. 

2.  EP  except  for  use  as  a  constant  (not  functional)  expansion  parameter; 
and  EPTERM,  except  for  use  to  control  expansions  in  EP. 

3.  01  or  02  except  for  use  as  metric  tensors. 

4.  W  except  for  use  as  stress-energy  tensor. 

5.  Any  name  beginning  with  a  digit  (0,  ...,  9);  i.e.,  all  names  must  begin 
with  Latin  letters. 

6.  Any  name  beginning  with  the  following  (these  are  forbidden  in  all 


F0RMAC  programs): 

ALGC0N 

DOUBLE 

I^GICA 

ATOMIC 

END 

MATCH 

AUTSIM 

ERASE 

|6rder 

BCDC0N 

EVAL 

PARAM 

CENSUS 

EXPAND 

PART 

C0MP12 

FIND 

REAL 

C0DEM 

EMC 

SUBST 

C0EFF 

IUNCTI 

SUBRKiU 

DEFEND 

IMC 

SYMARG 

DIMENS 

INTEOE 

LET 

11 
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7.  The  following  hums,  which  are  used  for  various  other  purposes  by  ALBERT 


ADS,  AA,  AAA 
BB0I1I,  BLANK 
CARD,  CCS 
Dl,  D2,  KEFF 
at,  exp 

I,  J,  K,  L,  M,  H 

II,  JJ,  JJJ,  LL 
LIKE,  LBL 

MCS,  MU,  HU 


HSUB8,  HXXP 

FP,  Q,  B 

Rica,  RTEN,  no 

8,  SO,  81,  85 

SCAL 

SUM 

Tl,  T2,  T3,  IT 
ZO,  Zl,  Z2 


FOntAC  requires  a  huge  amount  of  core  storage  for  its  operations.  Con¬ 
sequently,  if  one  inputs  metric  tensors  which  are  too  complicated,  he  will 
exhaust  the  memory  capacity  of  the  computer  before  finishing  the  computation. 
When  this  happens,  the  only  way  of  using  the  ALBERT  programs  as  they  stand  1b 
to  simplify  the  form  of  the  metric  tensor.  The  simplified  Ricci  and  Einstein 
tensors  which  result  therefrom  can  be  reexpressed  in  terms  of  the  more  compli- 
cated  metric  by  means  of  the  program  SUBSTITUTE .  For  examples  of  metrics  which 
did  and  did  not  compute  without  simplification,  and  of  the  simplifications 
used  for  certain  complicated  metrics,  see  {}  III.B.4,  HI.C.4,  IH.D.4,  and 
HI  .E.4. 

One  of  the  moot  effective  ways  to  avoid  exhausting  the  computer's  memory 
is  to  avoid  using  conplicated  denominators  in  the  expressions  for  the  metric 
tensor.  FQRMAC  does  not  handle  denominators  efficiently.  For  example 


01(1,1)  -  W(i.-2*a/x), 

where  04  depends  on  X,  is  much  more  likely  to  produce  trouble  than 
01(1,1)  -  FMCEXP(ALAM), 


where  A1AM  depends  on  X. 
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2.  Complete  Listing  of  EINSTEIN 

In  the  listing  of  Einstein  which  follows,  and  In  all  other  program 
listings  in  this  report,  all  control  cards  are  omitted.  (Control  cards 
are  cards  which  precede  and  follow  each  program  and  which  tell  the  computer 
such  things  as  who  the  user  of  the  program  is  and  what  language  the  program 
is  written  in.)  Because  control  cards  vary  in  format  from  one  computer 
installation  to  another,  the  user  should  consult  his  friendly  computer 
representative  for  special  instructions  on  their  use.  The  user  should 
warn  his  computer  representative  that  he  is  using  FOHMAC,  which  requires 
somewhat  different  control  cards  than  other  languages. 
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ilBFHC  ONE  NODECK 
SYMARG 
C 

C  EINSTEIN,  A  PROGRAM  TO  CALCULATE  THE  RICCI  AND  EINSTEIN  TENSORS 


C 

C 

c 

c 

c 

c 

c 

c 

c 


ATOMIC  X,Y, Z,T,EP, ANU.ALAM 
DEPEND  ( ANU, AL AM/X ) 


} 


Uaeri  Replace  these  carts 
with  your  own  input. 


ATOMIC  AND  DEPENO  STATEMENTS  MUST  BE  SUPPLIED  BY  USER 
IMMEDIATELY  PRECEDING  THIS  COMMENT. 

THE  ATOMIC  STATEMENT  MUST  LIST  ALL  FUNCTIONS  AND  INDEPENDENT 
VARIABLES  USED  IN  A  PARTICULAR  PROBLrM. 

THE  DEPEND  STATEMENT  DEFINES  ALL  THE  FUNCTIONAL  DEPENDFNCIES 
FOR  THE  PROBLEM,  I.E.  WHICH  FUNCTIONS  DEPEND  ON  WHICH  VARIABLES. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 

c 

c 

c 

c 

c 

c 

c 


REAL  MCS ( 4, 4, 4 1 
LOGICAL  Q 

DIMENSION  CARD112),LINEI12),CCSI4,4,4),G1I4,4),G2I4,4),IH4) 

DIMENSION  RICCI (4  I 
INTEGER  EPTERM 


IF  EPTERM  =  0  TERMS  CONTAINING  ALL  POWERS  OF  EXPANSION 
PARAMETER  WILL  BE  RETAINED. 

IF  EPTERM  .NE.  0  TERMS  CONTAINING  THE  EXPANSION  PARAMETER 

,  r.„r  f r\A  T0  A  P3HER  GREATER  THAN  2  WILL  BE  DISCARDED. 

A  CARD  CONTAINING  THE  INTEGER  VALUE  OF  EPTERM  MUST  FOLLOW 
THIS  COMMENT. 


EPTERM  =  0 


User:  Replace  this  card 
i  with  your  own  input. 


DO  I  I  =  1,4 
00  1  J  =  1,4 
LET  GIII.JI  =  0. 
LET  G2I l,J)  =  0. 
CONTINUE 


ALL  ELEMENTS  OF  THE  COVARIANT 
CUNTRAVAR I  ANT  METRIC  TENSOR, 
TO  ZERO. 


METRIC  TENSOR,  G1(I,J),  AND  THE 
G2II.J),  HAVE  BEEN  SET  EQUAL 


USER  MUST  SUPPLY  "LET"  STATEMENTS  DEFINING  ALL  NON  ZERO 
ELEMENTS  OF  THESE  TENSORS  IMMEDIATELY  FOLLOWING  THIS  COMMENT. 


LET 

Gl( 

1 

rl) 

LET 

Oil 

2, 

>2) 

LET 

Gl( 

3 1 

>3) 

LET 

G1I 

4, 

,4) 

LET 

G2I 

li 

.1) 

LET 

G2( 

2, 

,2) 

LET 

G2I 

3, 

,31 

LET 

G2I 

4, 

,4) 

-FMCEXPI AL  AM) 

-X**2 

-IX*FMCSINIY))**2 
FMCEXPI ANU) 
-FMCEXPI -ALAM) 
-1/X**2 

-1/IX*FMCSIN(Y))**2 
FMCEXPI -ANU) 


l  User:  Replace  theee  cards 
f  with  your  cm  input. 


J 
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THE  COVARIANT  METRIC  TENSOR  --  OUTPUT 
WR l TE ( 6  *  503 ) 

503  FORMAT (  51H1C0VARIANT  METRIC  TENSOR  NON  ZERO  ELEMENTS  (INPUT)  I 
DO  2  I  *  I » 4 
DO  2  J*l,4 

LET  Q  *  MATCH  ID,  G1(I,J),0. 

IF  ( Q  )  GO  TO  2 
MR  I  T  E ( 6 ,800  I  I » J 
800  FORMAT ( 1H  I5.IHL  I5.1HL  ) 

BEGIN  *  O. 

46  LET  BEGIN  *  BCOCON  01 <  I ,  J > , L INE, 12 

WRITE(6,200)(LINE(L),L=2,12> 

IFIBEGIN.NE.O. )G0  TO  46 
2  CONTINUE 

THE  CONTRAVAR I  ANT  METRIC  TENSOR  —  OUTPUT 
WR  I  TE  I  6 »  700  I 

700  FORMAT! 54HOCONTR AVAR  I  ANT  METRIC  TENSOR  NON  ZERO  ELEMcNTS  (INPUT)) 
DO  47  I  =  1,4 
DO  47  J  =  1 , 4 

LET  Q  »  MATCH  ID,  G2(I,J),C. 

IF ( Q ) GO  TO  47 
WR  f  TE ( 6,500 )  I,J 
500  FORMAT ( 1H  I5,1HU  I5,1HU  ) 

BEGIN  *  0. 

48  LET  BEGIN  =  BCOCON  02  (  I ,  J ) , L I NE , 12 
WR I TE (6,200)1  LINE (L),L  =  2, 12) 

IFIBEGIN.NE.O. 1 00  TO  48 

47  CONTINUE 

THE  COVARIANT  CHR l STOFFEL  SYMBOLS 
WRITE(6»201) 

201  FORMAT! 34H1THE  COVARIANT  CHR 1 STOFF EL  SYMBOLS  ) 

00  3  I  =  1,4 

DO  3  J  *  I  ,  4 

00  3  K  =  1,4 

LET  Tl  =  0. 


GO 

TO 

(31,32,70,711,1 

31 

LET 

Ti 

=  FMCDIF(G1<J,K),X,1> 

GO 

TO 

33 

32 

LET 

Tl 

=  FMCDIF(G1(J,K),Y,1) 

GO 

TO 

33 

70 

LET 

Tl 

=  FMCOIF(Gl(J,K),Z,l) 

GO 

TO 

33 
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71  LET  T1  =  FMCDIF(G1!J,K), T,l) 

33  CONTINUE 

LET  T2  =  0. 

GO  TO( 35 « 36* 72» 73)  *J 

33  LET  T2  *  FMCDI F ( Gl 1 1 , K I  *  X , 1 ) 

GO  TO  34 

36  LET  T2  =  FMCD I F ( GI ( I , K J ,  Y, 1 ) 

GO  TO  34 

T2  LET  T2  =  FMCDI F < Gl (I • K ) , Z , l ) 

GO  TO  34 

73  LET  T 2  *  FMCD I F I Gl ( I , K> , T , 1 ) 

34  CONTINUE 

LET  T3  =  0. 

GO  TO  (38, 39, 74,75), K 

38  LET  T3  =  FMCDI F ( Gl ( I , J > , X , 1 ) 

GO  TO  37 

39  LET  T3  =  FMCD I F < Gl I  I , J)  ,  Y , l ) 

GO  TO  37 

74  LET  T3  *  FMCDI F( Gl( I , J) ,  Z , 1 ) 

GO  TO  37 

75  LET  T3  =  FMCDIF < Gl (  I ,  J ) ,  T ,  11 

37  LET  CCS (  I , J,K )  -  0. 5* ( T1 *T2- T3 ) 

LET  CCS  I J, I ,K )  *  CCS ( I , J ,K ) 

LET  Q  *  MATCH  ID,  CCSI I , J , K ) , 0. 

I F  (  Q ) GO  TO  3 
HR  ITE(6, 203)1, J,K 

203  FORMAT ( 1H  I5.2HLS  I3.2HLS  I4.1HL) 

LET  ANS  =  EXPAND  CCS(I,J,K) 

LET  ANS  *  ORDER  ANS,!NC,FUL 
BEGIN  *  0. 

5  LET  BEGIN  *  BCDCUN  ANS,LINE,12 
WRITE(6»200)(LINE(L)»L*2,12) 

200  FORMAT! 1H  5X  12A6) 

IFIBEGIN.NE.O. )G0  TO  5 
ERASE  ANS 
3  CONTINUE 

ERASE  Ti,T2,T3 

THE  MIXED  CHR I STOFFEL  SYMBOLS 
WRITE! 6,205 ) 

205  FORMAT!  30H1THE  MIXED  CHRISTOFFEL  SYMBOLS  ) 
DO  6  I  =  1,4 
DO  6  J  *  1,4 
DO  6  K  »  1,4 
LET  MCS! I, J,K)  =  0. 

DO  7  L  =  1,4 

LET  TT  =  G2(K,L)*CCS! I,  J,L) 

LET  TT  =  EXPAND  TT 
IF(EPTERM)506,507,506 
507  LET  SUM  =  TT 
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ERASE  TT 
GO  TO  514 
506  CONTINUE 

LET  SI  *  COEFF  TT,EP**0,R 

LET  S2  *  COEFF  TT,EP**1,R 

LET  S3  =  COEFF  TT,EP**2.R 

ERASE  TT 

LET  SUM  =  SI  ♦  S2*EP  ♦  S3*EP*EP 
ERASE  SI , S2 «  S3 
514  CONTINUE 

LET  MCSII.J.K)  =  MCSII.J.K)  ♦  SUM 
ERASE  SUM 
7  CONTINUE 

LET  MCSI I , J.K  )  =  EXPAND  MCSII.J.K) 

LET  MCS(J.I.K)  =  MCSII.J.K) 

LET  Q  =  MATCH  If).  MCSI  I  .J.K)  ,0. 

IFIOIGO  TO  6 
WRITEI6.769)  I, J.K 

769  rORMATIlH  15.2HLS  n,2HLS  I4.1HU  ) 

LET  ANS  *  ORDER  MCi ! I , J ,K) , INC . FUL 
BEGIN  *  0. 

B  LET  BEGIN  =  BCOCON  ANS, LINE. 12 
WRITE(6»?00)(LINEIL),L=2»12) 

IF  I  BEGIN. NE.O. ) GO  TO  B 
ERASE  ANS 
6  CONTINUE 
C 
C 

00  350  1=1,4 
DO  350  J= 1 , 4 
DO  350  K= 1 , 4 
ERASE  CCSII.J.K) 

350  CONTINUE 

C  THE  RICCI  TENSOR 

C 

WRITEI6.217) 

217  FORMAT  1 36H1 THE  RICCI  TENSOR,  MIXED  COMPONENTS  ) 


DO 

50 

N  =  1,4 

DO 

50 

K  *  1,4 

LET 

SI 

=  0. 

DO 

51 

I  =  1,4 

DO 

51 

M  =  1,4 

LET 

T 1 

=  0 . 

GO 

TO 

(?1»22,23,24),I 

21 

LET 

T 1 

=  FMCDIFIMCSIN.M,  I  (,X,1) 

GO 

TO 

19 

22 

LET 

T 1 

=  FMCDIFIMCSIN.M,  I  I.Y. 1) 

GO 

TO 

19 

23 

LET 

T 1 

=  FMCDIFIMCSIN.M,  I  ),Z, 1) 

GO 

TO 

19 

24 

LET 

Tl 

=  FMCDIFIMCSIN.M,  I  >,T,1) 
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19  CONTINUE 

LET  T1  -  EXPAND  T1 
LET  T2  =  0. 

GO  TO  ( 52(53* 54, 55  I ,  N 

52  LET  T2  =  -FMCO IF ( MCS I  I , M , I) ,  X  ,  1 ) 

GO  TO  56 

53  LET  T2  *  -FMCDIF ( MCS ( I ,M, I ) , Y ,  1 ) 

GO  TO  56 

54  IET  T2  *  -FMC0IF(MCS(I,M,I),Z,1» 

GO  TO  56 

55  LET  T2  =  -FMCDIFIMCSI 1,M, I) ,T,1J 

56  CONTINUE 

LET  T2  *  EXPANO  T2 

LET  TT  =  G2IK,MI*ITH-T2  > 

ERASE  Tl,T2 
LET  TT  *  EXPAND  TT 
IF(EPTERM)504,513,504 
513  LET  SI  =  SI  *  TT 
ERASE  TT 

LET  SI  ■  EXPAND  SI 
GO  TO  51 
504  CONTINUE 

LET  SO  =  COEFF  TT,EP**D,R 
LET  S3  =  COEFF  TT,EP**1,R 
LET  S4  =  COEFF  TT,EP**2,R 
ERASE  TT 

LET  SI  =  SI  ♦  ISO  «-  S3*6P  ♦  S4*EP*EPI 
ERASE  S0»S3*S4 
LET  SI  =  EXPAND  SI 
51  CONTINUE 

LET  S2  =  0. 

DO  84  I  =  1,4 

DO  84  J  =  1,4 

DO  84  M  ■  1,4 

LET  TT  =  G2IK,M)*(-MCS(J,M,I)*MCSII,N,JH-MCS<M,N,!>* 

1  MCS(I,J,J>) 

LET  TT  =  EXPAND  T1 
IFIEPTERM 1510*511*510 
511  LET  S2  =  S2  +  TT 
ERASE  TT 

LET  S2  =  EXPAND  S2 
GO  TO  84 
510  CONTINUE 

LET  SO  *  COEFF  TT,EP**0,R 
LET  S3  =  COEFF  TT,EP**1,R 
LET  S4  *  COEFF  TT,EP**2,R 
ERASE  TT 

LET  S2  =  S2  ♦  (SO  *  S3*EP  ♦  S4*EP*tP) 

ERASE  SO ,  S3 , S4 
LET  S2  =  EXPAND  S2 
84  CONTINUE 
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LET  RTEN  =  SI  ♦  S2 
ERASE  SI.S2 
LET  RTEN*  EXPANO  RTEN 
IF  (N-K)  983,984,983 
984  LET  RICCI (N)  *  RTEN 

983  CONTINUE 

LET  0  =  MATCH  ID, RTEN, 0. 

IF (01  GO  TO  50 

LET  ANS  =  ORDER  RTEN, INC, FUL 
HR  I TE ( 6 ,216 ) K,N 

216  FORMAT ( 1H  I5.1HU  I5,1HL  I 

501  FORMAT (215) 

BEGIN  =  0. 

58  LET  BEGIN  *  BCDCON  ANS, LINE, 12 

WRITE(6,200)(LINE( I L ) , I L  *  2,12) 

502  FORMAT ( 12 A6 ) 

IF(BEGIN.NE.O.)GO  TO  58 
ERASE  ANS 
50  CONTINUE 
92  CONTINUE 

THE  RICCI  SCALAR 

LET  SCAL  =  RICCII1)  ♦  RICCK2I 
LET  SCAL  =  EXPAND  SCAL 
LET  SCAL  *  SCAL  ♦  RICCI (3) 

LET  SCAL  *  EXPAND  SCAL 
LET  SCAL  =  SCAL  ♦  RICCI  (4) 

LET  SCAL  *  EXPAND  SCAL 

LET  ANS  =  ORDER  SCAL, INC, FUL 

HRITE  (6,2331 

233  FORMAT  (  17H1THE  RICCI  SCALAR) 
regin  =  n. 

234  LET  BEGIN  =  BCDCON  ANS, LINE, 12 

HRITE (6,2001 (LINE ( IL),  I  L  =  2,12) 

IF  (  BEGIN. NE  >0.  I  Gil  TO  234 
ERASE  ANS 

THE  EINSTFIN  TENSOR,  DIAGONAL  COMPONENTS 


HR  I  T  E ( 6, 22 1 ) 

221  FORMAT!  5BH1THE  EINSTEIN  TENSOR,  MIXED  COMPONENTS, 
ITS  ) 

DO  516  L=1 ,4 

LET  ANS  =  RICCI (L  )  -  0.5*SCAL 
LEI  ANS  =  EXPANO  ANS 
LET  ANS  =  ORDER  ANS,  INC,  FUL 
WRITE(6, 219)1, L 
219  FORMAT ( 1H  15, 1HU  I5.1HLI 
BEGIN  =  0. 

342  LET  BEGIN  =  BCDCON  ANS,  LINE,  12 


DIAGONAL  ELEMhN 
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WRI TE ( 6 *200 )  ( L I NE ( 1 1 ) » 1 1»2  > 1 2 ) 
IF( BEGIN. NE .0.1  GO  TO  342 
ERASE  ANS 
516  CONTINUE 
93  CONTINUE 
STOP 
END 
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3.  Output  for  BUteTEiN 

The  printed  output  for  KXNSTTEIR  includes:  the  metric  c  exponents,  g^j 
and  which  were  input  hy  the  user;  the  Chrlstoffel  symbols,  which  ere 
calculated  from  the  formulas 


rijk  “  $  *®ik,J  +  8Jk,i  " 


(2) 


t  m 

ij 


(3) 


(note  that  the  Chrlstoffel  symbols  are  symmetric  In  the  first  two  indices); 
tht  mixed  components  of  the  Ricci  curvature  tensor,  calculated  from 


l 


00 


the  Ricci  scalar. 


(5) 


and  the  diagonal  mixed  components  of  the  Einstein  tensor, 

0*,  =  R*  -  i  R  »k  .  (6) 

n  n  *  n 

(Recall  that  the  off-diagonal  components  of  the  Einstein  tensor  are  Identical 
to  the  off-diagonal  components  of  the  Ricci  tensor.)  Components  which  are 
aero  are  not  printed  out;  and  when  two  conponenta  are  equal  by  symnetry  con¬ 
siderations  (e.g.,  *  r^i^) ,  only  one  of  them  is  printed. 

The  format  of  the  printed  output  is  self-explanatory  except  for  the 
following:  Before  each  component  of  a  quantity  the  computer  prints  the 
numerical  values  of  the  indices,  plus  symbols  to  identify  whether  the  indices 
sire  contravariant  (U  for  upper)  or  covariant  (L  for  lower).  For  example,  in 
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the  Ricci  tensor 

UJ  3L  «— »  R^j  . 

Moreover,  in  the  Cbrlstoffel  symbols  an  8  is  used  to  Identify  the  indices 
which  are  syantric: 

1I£  5U3  4L  «-»  r^,  symmetric  under  interchange  of  1  and  3. 

For  metrics  which  are  too  ccnflicated,  the  computer  will  exhaust  its 
memory  capacity  before  collating  the  computations .  When  this  occurs,  the 
computer  prints  out  "ERROR  TRACE.  CAU£  IN  REVERSE  ORDER.",  followed  by  a 
list  of  routines  being  used  at  that  moment,  and  then  followed  by  "FOHMAC 
ERROR  NUMBER  '7537'  HAS  OCCURRED.  EXECUTION  TERMINATED."  The  computer  then 
terminates  the  confutation.  Other  error  messages  any  occur  as  a  result  of 
faulty  input.  In  order  to  interpret  these  error  messages  the  reader  should 
refer  to  pp.  213-222  of  the  FOHMAC  manual  (IEM  1964). 

As  an  example  of  successful  output  from  a  confutation  by  EEI8TEIN,  we 
give  below  the  results  far  the  static,  spherics  lly-syinetric  metric  of  equa¬ 
tion  (l).  The  input  which  genera'  Jd  these  results  was  given  in  {  m.B.l 
and  also  in  the  listing  of  EINSTEIN  in  ft  III.B.2. 
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COVAHI ANT  METRIC  TENSOR  NON  ZERO  ELEMENTS  I  INPUT) 

1L  1L 

-FMCEXPI ALAMJt 
2L  2L 

-X**2.0t 
3L  3L 

-X**2.0*FMCSIN(YI**2.0t 
4L  4L 

FMCEXPI ANU ) $ 

CONTRAVARI ANT  METRIC  TFNSOR  NON  ZERO  ELEMENTS  I  INPUT I 
1U  1U 

-FMCEXPI -ALAM ) t 
2U  2U 

-X**(-2.0lt 
3U  3U 

-X**(-2.G)*FMCSlN(Y)**(-2.0)t 
4U  4U 

FMCEXPI-ANUI t 


THE  COVARIANT  CHR I  STfJFFEL  SYMBOLS 


1LS 

1LS 

1L 

-5. 

OE-l«TMCEXPl ALAM  l*FMCDI FI ALAM, 

<X,I))S 

US 

2LS 

2L 

-xt 

1LS 

3LS 

3L 

-X*FMCSlN(YI**2.0t 

US 

4LS 

4L 

5.0H-1*FMCEXPIANU>*FMCDIFI ANU, IX, 

lilt 

2LS 

2LS 

1L 

Xt 

2LS 

3LS 

3L 

-X**2.0*FMCSIN( Y)*FMCCOS( Ylt 
3LS  3LS  1L 

X*FMCSlNIY)**2.Gt 
3LS  3LS  2L 

X**2.0*FMCSIN| Y)*rMCCOS(Y)t 
41 S  4LS  1L 

-5.0E-1*FMCFXPI ANU)*FMCDIFI ANU, IX, l)  It 


the  mixed  christoffel  symbols 

1LS  US  1U 

5.0C-1*FMCDIF I ALAM, (X,ll It 
1LS  2LS  2U 

X**(-1.0lt 
1LS  3LS  3U 

X**(-l .0) t 
1LS  4LS  4U 

5 • Ot- 1 *F  MCDI F I ANU, IX,1  I  It 
2LS  2LS  UJ 

-X*FMCEXPI-ALAM)t 
2LS  3LS  311 

FMCSIN(Y)**l-l.f  I *FMCCOS  I Y ) t 
3LS  3LS  IU 

-X*FMCSIN(YI**2.0*FMCEXP(-ALAM)t 
3LS  3LS  ?U 

-FMCSINIYI*FMCCI)S(Ylt 
4LS  4LS  1U 

5.0b-l*FMCEXPI-ALAM+ANUI*FMC0IFIANU, IX, II  It 
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THfc  RICCI  TENSOR,  MIXED  COMPONENTS 
IU  1L 

-X**I-1.0)*FMCEXPC-ALAMI ♦FMCDIFI ALAM, (X, l  II  -  2 . 5E-KFMCEXP I -AL  AMI  * 
FMCDIFIALAM,  (X,l>  I  *F  MCOI F  (  ANU,  (  X  *  1  I  )  *2 . 5E-  1*FMCEXP< -ALAMI  *FMCDIF  I 
ANU, IX, II l**2.0+S.0L-l*FMCEXP(-ALAH)*FMCDIFC  ANU, (X, 2  I It 
2U  2L 

X**(-2.0)*FMCEXPI-ALAM)-X**I-2.0)-5.0E-1*X**(-1.0I*FMCEXP|-ALAMI * 
FMCDIFIALAM,  IX.II I *5 .OE-1 *X**( -1 .0 1 *FMCEXP< - ALAM I *FMCDI F I ANU, I  X, 1 1 
It 

3U  3L 

X**(-2.0I*FMCEXP(-ALAM)-X**{-2.0I-5.0E-1*X**(-1.0)*FMCEXP(-ALAMI* 
FMCDIFIALAM,  (X,  11 I *5 .0E-1*X** I -l .0 )*FMCEXPI - ALAM I *FMCDIF I ANU, I  X, 1 1 
It 

4U  4L 

X**(-1.0I*FMCEXP|-ALAMI*FMCDIF(ANU,(X,1 )}-2. 5E-1 *FMCE XP < - AL AM  I • 
FMCOIF I ALAM,  ( X,  1 1 ) *FMCD l F ( ANU, I  X , l 1 1 +2 . 5E-1*FMCE XP( -ALAM) *f  MCOIF I 
ANU, tX,l))**2.0+5.0E-l*FMCEXP(-ALAM)*FMCDIF(ANU»IX,2l It 


THE  RICCI  SCALAR 

2.0*X**(-2.0I*FMCEXP(-ALAM}-2.0*X**<-2.0I-2.0*X**(-1.0)*FMCEXP(- 
ALAM I *FMCDIF (ALAM, (X, II  I +2 .0*X** I - 1 .0  I *FMCEX PI -AL AM  I *FMCD I F I ANU, IX 
,11 I-5.0E-1*FMCEXP(-ALAM1*FMCD1F(ALAM, I  X, 1 1 I *FMCDIF< ANU,  I  X , 1 ) }  + 
5.0t-l*FMCEXP(-ALAMI*FMCniFIANU,IX,ll l**2.0«-FMCEXP(-ALAMI  *FMC01F < 
ANU, (X, 2)1$ 


THE  EINSTEIN  TENSOR,  MIXED  COMPONENTS,  DIAGONAL  ELEMENTS 
1U  1L 

-X**I-2.0)*FMCEXP(-ALAM|*X**|-2.0I-X**I-1.0I*FMCEXP(-ALAM1*FMCDIF( 
ANU, IX, lilt 
2U  2L 

5.0E-m*»(-1.0l*FHCEXP(-ALAMI*FMCDlF(ALAM,  (X,  1 1  1-5. OE-l *X ** I -1 . G I 
*FMCEXP(-ALAM)*FMCt)IKANU,  IX,1 1  I  *2. 5E-1  *FMCE  X  P( -ALAM}  *FMCD  I F  ( ALAM , 
(X,l I  I *FMCDI  F I ANU , (X,ll I -2 .5E- 1 *FMCE XP( - AL AM  I *FMCDIF I ANU,  I  X , 1 }| ** 
2.0-5.0L-1*FMCEXP(-ALAM»*FMCOIF( ANU,  I  X,2» It 
3U  3L 

5. OE-l *X **(-!. 0»*FMCFXP(-AL AM l*FMCDIF(ALAM, I  X  ,  1 1  I -5.0E-1*X** | -1 . 0  I  | 
*FMCEXP(-ALAM)*FMCD1F1  ANU,  (X,l  I  I  «-2 . 5E- 1  *FMCE  X  PI  -  ALAM  I  ♦FMCO  IF  I  ALAM,  j 
I  X , 1 1  I *FHCDI F I ANU, I  X, 1  I  I -2 . 5E-l*FMCEXP I -AL AM > *FMCOIF I ANU,  I  X , 1 1  1 ** 
2.0-5.0E-l*FMCEXP(-ALAM)*FMCDIF(ANU, (X,2I It 
4U  AL 

~X**I-2.0IPFMCEXP(-ALAM)*X**I-2.0>+X**(-1.0)*FMCEXP(-ALAM)*FMCD!FI  i 
ALAM, (X, lilt  \ 
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k. _ Staple  Problems  for  EXNSTEIM 

Input  end  output  for  one  sample  problem,  the  etatie  spherically  symnetric 
metric  of  equation  (l),  were  given  In  the  last  two  sections.  As  a  second  and 
such  more  complicated  example,  consider  the  metric  for  a  star  In  even-parity, 
nonradlal  pulsation,  as  discussed  by  Thorne  and  Campolattaro  (1967): 


ds 


2 


ev  (1  +  Hq  P^dt2  +  2  Hx  Pi  dtdr  -  eX  (l  -  Hg  P|)  dr2 
-  r2  (1  -  K  Pi)  (dfi2  +  sln20d(p2). 


(7) 


Here  v  and  X  are  functions  of  r;  HQ,  H^,  Hg  and  K  are  functions  of  r  and  tj 
and  Pj(cobS)  Is  the  Legendre  polynomial  of  order  t. 

To  prepare  Input  for  computing  the  Ricci  and  Einstein  tensors  of  the 
metric  (2),  we  make  the  following  changes  to  FORMAC  notation: 


t  -*•  T 

r  ♦  X 

9  ♦  y 

<p  ♦  Z 


(these  are  the  standard  coordinates  as  accepted 
by  ALBERT). 


X  ♦  ALAM 
V  ♦  AMU 
Hq-  H° 

Hx  +  HI 
H2  -  H2 
K  -  AK 


AL 


(these  names  were  chosen  arbitrarily,  subject 
only  to  the  rules  outlined  In  §  III.B.l.) 
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Also,  we  make  the  decision  to  let  the  conputer  treat  P^(cosO)  •  AL(Y) 
as  an  unspecified  function.  Finally,  because  we  want  the  Bled  and  Einstein 
tensors  only  to  first  order  (linear)  in  the  perturbation  from  equilibrium, 
we  attach  the  expansion  parameter  EP**2  to  the  metric  perturbations,  AK,  HO, 
HI,  and  H2,  wherever  they  appear.  The  resultant  input  and  output  are  as 
follows: 

INTO!: 


ATOMIC  X  »Y»7»T»FP*AL»ANU»ALAM*AK*H0*Hl*H2 
PFPFND  ( AMU* ALAM/X ) *  IAL/Y) 

OFPFND  (HO. HI .H2.AK/X.T) 


FPTERM  =  3 


LET 

Gin.i ) 

= 

LET 

Gin. 4) 

S 

LFT 

G  1  ( 4 . 1  ) 

S 

LET 

G1  (2,2 ) 

= 

LFT 

G1  (3.3 ) 

= 

LFT 

G 1 (4 .4  ) 

= 

LFT 

G2  ( 1 , 1  ) 

= 

LFT 

G2I1.4) 

s 

LET 

G2  ( 4 » 1  ) 

= 

LFT 

G2I2.2 I 

e 

LET 

G2 (3 . 3 ) 

= 

LET 

G2 ( 4 .4 ) 

= 

-  ( -AL*EP**2*H2+1 * ) aFMCEXP ( ALAM) 
AL*EP*»2»H1 

Gltl. 41 

-(l.-AK*AL*EP**2  >*X*«2 
-n.-AK*AL»EP*»2)*X**2*FMCSIN(Y)**2 
( AL*CP**2*H0+1. ) *FMCEXP( ANU) 
-(AL*EP*»2*H2+1. ) *FMCEXP ( -ALAM ) 
AL*FP**2*H1*FMCEXP ( -ALAM-ANU) 

G2i 1.4) 

-( AK#AL*EP**2  +  1. )*X**(-2) 

-  ( AK*AL*EP**2+1 . ) * ( X*FMCS IN  (  Y ))  •*  ( -2  I 
( -AL*EP*»2*H0+1 . ) *FMCEXP ( -ANU ) 


POTTO: 

See  following  pages 
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COVARIANT  METRIC  TENSOR  NON  ZERO  ELEMENTS  I  INPUT) 

IL  IL 

-|-AU*EP**2.0*H2+I.O )*FMCEXP( ALAM) t 
IL  AL 

AL*EP**2.0*Hlt 
2L  2L 

-(-AK*AL*EP**2.0+1.O)+X**2.OJ 
3L  3L 

-(-AK*AL*EP**2.O+1.0)*X**2.0*FMCSIN(Y)**2.Oi 
AL  IL 

AL*EP+*2.0*Hlt 
AL  AL 

( AL*EP**2.0*H0+l.O)*FMCEXP( ANU )  t 

CONTRAVAKI ANT  METRIC  TENSOR  NON  ZERO  ELEMENTS  (INPUT) 
1U  1U 

-( AL*EP**2.0*H2+1 .0) *FMCEXP(-ALAM)$ 

1U  AU 

AL*E P**2. 0 +H 1 *FMCfcXP< -ALAM- ANU )$ 

2U  2U 

-( AK*AL*EP**2 .0+1 .0) *X**(-2.0 I* 

3U  3U 

-<AK*AL*EP**2.0+1.0)*X**(-2.0)*FMCSIN<  Y)**(-2.0)» 
AU  iU 

AL*hP**2.U*Hl*FMCtXPI-ALAM-ANU)$ 

AU  AU 

(-AL*EP**2.O*H0+l.O)*FMCEXP(-ANU)$ 


THE  COVARIANT  CHR1ST0FFEL  SYMBOLS 
1LS  1LS  IL 

5.0E-1*AL*EP**2.G*H2*FMCEXP( AL AM  I *FMCDI F ( ALAM, (X, 1)  ) +  5. OE- 1  * AL*EP 
**2.0*FMC£XP( ALAM)*FMC0IF(H2. (X, 1 ) )-5 . OE-1 *F MCE XP < ALAM) *FMCDIF ( 
ALAM, (X.I  )  )t 
1LS  IL  S  2L 

-5. jE-l*EP**2.0*H2*FMCEXP( ALAM) *FMCDIF ( AL,  (  Y , 1 ) )$ 

1LS  1LS  AL 

-b.OE-I*AL+EP**2.0*FMCEXP(ALAM)*FMCDIF(H2,  ( T , l ) ) ♦ AL*cP**2 . O+FMCD I F 
(HI, (X, nit 
1LS  2LS  IL 

5 .Ot-l*EP**2 .0*H2*FMCEXP( ALAM) *FMCDIF ( AL, ( Y, 1 ) ) t 
1LS  2LS  2L 

AK*AL*EP**2.0*X+S.0£-I*AL*EP**2.0*X**2.0*FMCDIF (AK,(X,1))-Xt 

1LS  2LS  AL 

5.0E-1*EP**2.0*H1*FMCDIF( AL, (Y,l) )* 

1LS  3LS  3L 

AK*AL*EP**2.0*X*FMCSIN(Y) **2. 0+5.0E-1* AL*EP**2. 0*X**2 .0*FMCSIN( Y ) 
**2.0*FMCDIF( AK, (X, 1 ) >-X*FMCS  IN ( Y ) **2 .0 1 
1LS  ALS  IL 

5 ,0E-1*AL*EP**2.0*FMCEXP( ALAM)*FMCDIF ( H2, 1 T,  1 ) ) % 

LLS  ALS  2L 

-5.0E-1*EP**2.0*HL*FMCDIF( AL.IY.ll )$ 
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1LS  4LS  4L 

5.0E-1*AL*EP**2.G*H0*FMCEXP< ANU  >  *F MCI) I F  I  AMU,  (X,  1 ) ) ♦ 5. Ob- i *AL*E P** 
2.0*FMCEXPIANU)*FNCDIF(H0,<X, 1) i+5.UC-l*FMCtXP(ANUI*rMCDIF(ANU,  IX, 
1  )  )  V 

2LS  2LS  1L 

-AK*AL*tP**2.Q*X-5.0F-l.*AL*EP**?.G*X**2.0*FMCOIF(  AK,  I  X, 1 1 ) »Xt 
2LS  2LS  2L 

5.0t-l*AK*EP**2.0*X**2.0*FMCDlF( AL  , I Y,  1  )  )  t 
21 S  2LS  4L 

-5.0E-l*AL*EP**2.0*X**2.C*FMCniFIAK,( 1,1)1* 

2LS  3LS  3L 

AK*AL*EP**2.0*X**2.0*FMCSINI  Y  I  *FMCC(IS  1  Y  ) +t>.  0  F-l*  '  K*E  P  **  2  .  Q*X**2 . 0  ♦ 
FMCSIN(Y)**2.0*FMC01FIAL,  (Y, 1 ) ) -X**2 . 0*F MCS I N I Y I ♦FMCCOS ( Y )  l 
2LS  4LS  1L 

5.0L-1*EP**2.0*H1*FMCDIF(AL»(Y,1 ) I t 
2LS  4LS  2L 

5.0E-l*AL*EP**2.r.*X**2.C*FMCQ!FIAK»  I  1 ,  1 )  )  t 
2LS  4LS  *»L 

5.0L-l*tP**2.0*HC*FMCEXPI ANU)*FMCDIF( AL, I Y,1  ) I* 

3LS  3LS  1L 

-AK*AL*EP**2.0*X*FMCSINt  Y) **2.0-5. OE-  l*AL*EP**2. 0*X**2. 0*FMCSINI  Y I 

**2.0*FMC()IF ( AK, (X,l ) )+X*FMCSlN( Y I **2.0* 

3LS  3LS  21 

-AK*AL*EP**2 .0*X**2.U*FMCSINIY )*FMCCOS I Y) -5. OE-l*AK*FP* *2. i*X* *2.0 
*FMCSIN(Y)**2.0*FMCDIF(AL, IV, 1 ) I  *X**2. 0*FMCS  I N I  Y  )  *FMCC11S  <Y)  t 
3LS  3LS  4L 

-5.GE-1*AL*EP**2.G*X**2.0*FMCSIN( Y  1**2. )*FMCniF(  AK, IT , 1 )) t 

3LS  4LS  3L 

5.0E-1*AL*EP**2.0*X**2.0*FMCSINI Y ) **2 . 0*FMCD I F  (  AK , IT,  1) )t 
4LS  4LS  1L 

-5. OE- 1* AL*EP **2. 0*H0*FMC EXP I ANU ) *FMCOIF I ANU, IX, l ) )-i.OE-l*AL*EP** 
2.0*FMCEXPIANUI*FMCDIFIH0, (X, 1 ) l ♦ AL*EP**2 .0*FMCU I F ( HI , I T, l  I  >- 
5.0E-1*FMCEXP( ANU )*F MODI  FI ANU, IX, l  )  It 
4LS  4LS  2L 

-5.GE-l*fcP**2.0*HO*FMCEXP(ANU)*FMCDIF(AL, IY,  1) )t 
4LS  4LS  4L 

5.0E-1*AL*EP**2.0*FMCEXP( ANU )*FMCOIF1HO, I T ,  1  ) )* 


THE  MIXED  CHRISTC1FFEL  SYMBOLS 
1LS  1LS  1U 

-5.0E-i*AL*EP**2.0*FMCDIF(H2, IX,1 I ) +5 . OE- 1 *FMCDI  F I AL AM, CX, 111* 

1LS  1LS  2U 

5.0E-l*fcP**2.0*H2*X**(-2.0J*FMCFXPI  AL AM) *FMCDI F ( AL , IY,1 )  I  * 

1LS  US  4U 

-5.0E-1*AL*EP**2.0*H1*FMCEXPI-ANU)*FMCDIF( ALAM, IX, l I l-5.0E-l*AL*EP 
**2.0*FMCEXPI ALAM-ANU)*FMCDIFIH2, I T,  1 1 » ♦AL*EP**2.0*FMCEXPI-ANU»* 
FMCDIFIH1, IX, l  II* 

1LS  2LS  1U 

-5.0E-1*EP**2.0*H2*FMC0IFI AL, IY,1 I )* 

1LS  2LS  2U 

-5.0E-1*AL*EP**2 .0*FMC0IF  I  AK, I  X, 1 ) 
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1LS  2LS  4U 

!>.0C-l*tP**2.0*Ml*FMC.EXP(-ANU)*FMCDlF(AL,(Y,  X  >  It 
US  31 S  31) 

-5.0fc-l*AL*EP**2.0*PMCDlF( AK,(X,ll  »♦***! -1.01* 

1LS  4LS  1IJ 

5.0t-l*AL*EP**2.0*Hl*rHC6XP<-ALAHI*FMCDlF<ANU, (X,l) >-5.0b-.*AL*tP 
**2.0*FHC0IFIH2t  I  f,  l  ))S 
1LS  4LS  21) 

5.0l-1*£P**2.0*H1»X**(-2.0  )*FMCDfFIAL, ( Y. 1 ) ) t 
1LS  4LS  4U 

b.OE-l*AL*FP**2.0*FMr.OIF(HO,  (X,  1 II  ♦  5 . 0E- 1*FHCDI  F  I ANU,  ( X ,  1  ))  t 
2LS  2LS  1U 

AK4AL*EP**2.0*X*FMCE <P ( -AL AM ) -AL*l P**2 .0*H2* X*FMCEXPI -AL AM ) ♦ 

5.0E-] *AL*lP**2.0  +  X**2.0*FMCEXP<-ALAM|*FMCDIF(AK,  «X,1 ))-X*I HCEXPI- 
ALAH) t 

2LS  2LS  2U 

-5.0E-14AX*EP**2.0*EMCDIF» AL, (Y,l  I  )i 
2LS  2LS  4U 

AL*LP**2.0*H1*X*F^CE  XP ( -AL  AM- ANU )-5.0E-i*AL*EP**2.0*X**2.0<‘FMCEXP( 
-ANU )*FHCDIFIAK,(T,1 ) It 
2LS  3LS  31) 

-5.0t-l*AK*CP**2.0*FMC0lF  I AL, ( Y , 1 1 ) *F MCS I N( Y ) ** ( -1 . 0 ) *FHCCOS I Y 1 t 
2LS  4LS  IU 

-•>.OC-l*EP**2.O*Hl*FMCEXPt-ALAHi*FHC0IFI  AL,(  Y,  1 1  It 
2LS  4LS  21) 

-5.UL-1*AL*EP»*2.0*FHC01F<  AK, ( r,l  I  It 
2LS  4LS  4U 

5.0t-l*cP**2.0*H0*FMCDlFI AL, IY,l) )i 
3LS  3LS  1U 

AK4AL*EP**2.0*X*FMCS IN (YI **2 . 0*FMCEXP I -AL AH  I - AL*EP**2 .0*H2*X* 
FHCSINIYI**2.0*FMCEXPI-ALAM)*5.0E-1*AL*EP4*2.0*X**2.0*FMCSINI Y)** 
2.04FMCEXP(-ALAM»*FHCDIF« AK,IX,1I I-X*FMCSIN< Y I **2 ,0*FMCEXP < -AL AN  1 1 
3LS  3LS  2U 

5.0c-l*AK*EP**2.0*FMCSIN(Y)**2«0*FHCDIF(AL,(Y,l))-FHCSIN(Y)*FMCCOS 

1  Y|t 

3LS  3LS  4U 

AL*tP**2.0*Hl*X*FMCS  INI  Y  J  **2 . 0*FMCEXP  l  -  AL  AH-  ANU )  -5.0E-1  *AL  *EP*«>2 . 0 
*X**2.0*FHCSIN(Y)**2.0*FHCEXP(-ANUI*FMCDIF(AK, <T,1) It 
3LS  4LS  3U 

-5.CE-l*AL*LP**i..0*FMC0IF(  AK,  IT,  II  It 
4LS  4LS  XU 

5.0E-1*AL4EP*42.G*HO*FHCEXPI-ALAM+ANU)*FMCDIF(ANU,IX,1)  )♦ 
5.0t-l*AL*EP**2.0*H2*FHCEXP(-ALAH*ANU)*FHCDIF( ANU,IXf II I-AL*EP** 

2  »0*F  MCt XP I -AL AH  I AFHCDIF (  HI  *  I T  *  1 1 1 *5  »0E-l*AL*EP**2«0*FMCt  XPI-ALAH* 
ANU)*FHC01F(H0t ( X,l I  I »5 .OE- 1  *FHCEXP I  - ALAM+ANU I *F HCDIF ( ANU, I X , III t 

4LS  4LS  2U 

5.0E-l*EP**2  .  G*HO*X**(  -2. 0  I  *FHCEXI’ ( ANU)*FHCD  IF  I  AL,  IV,  X)  It 
4LS  4LS  4U 

-5, uE-i*AL*EP**2.0*Hl *FHCEXP( -ALAH I *FMCI)I F ( ANU, I X, 1 1 1 ♦5.0E-14 AL*EP 
**2.0*FHCDIF(H0,IT,1 lit 
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THE  RICCI  TENSOR,  MIXED  COMPONENTS 
IU  II 

-AL*EP**2.0*H2*X** ( -l .0 ) *FMCEXP( -ALAM ) *FMCDI  F(ALAM,(X,1))- 
2 . 5E-I *AL *EP**2 .0  *H2  *FMCEXP( -AL AM  I *FMCD I F I AL AM, I  X, 1 ) I *F MCD I F ( ANU, ( 

X , 1 ) ) ♦2.5E-1*AL*EP**2.0*H2*FMCEXP(-ALAM)^FMCD IF ( ANU, ( X, 1 ) ) **2 ,0* 

5 ,OE-l*AL*EP**2.0*H2*FMCEXP( - AL AM  I *FMCD I F ( ANU, I X,2> I -2. 0* AL*EP** 
2.0*X**I-I.01*FMCEXP(-ALAMI*FNCOIFIAK,  IX,m»AL*EP**2.0*X**l-i.0l* 
FMCEXP(-ALAM) ♦FMCDIF (H2, ( X , 1  I I»5.0E-l*AL*EP**2#0*FMC£XP(-ALAM)* 

F MCD IF ( AK, (X, 1 ) ) ♦FMCD I F ( AL AM, (X, 1 ) |-AL*EP**2 • 0*F MCEXP ( - AL AM ) ♦ 

FMCDIF ( AK , ( X , 2 ) )-2.5E-l*AL*EP**2.0*FMCEXP(-ALAMI*FMCDIF(ALAM, C  X ,  1  ) 

) ♦FMCD I F ( HO , C  X , L ) I*5.0E-I*AL*EP**2.0*FMCEXPI -ALAM ) ♦FMCD IF ( ANU, I  X , l 
) >*FMCDIF(H0, <X,i ) U2.5E-L*Al*EP**2.0*FMCEXPl-ALAM»*FMCDIF(ANU,(X, 

1 )  l*FMCDIF(H2  » ( X, 1 1 ) ♦5.0E-1*AL*EP**2 .0*FMCEXP(-ALAM) ♦FMCD IF  CHO, I  X, 

2 )  )*5.0E“1*AL*EP**2.0*FMCEXP(-ALAM-ANU)*FMCD!F(ALAM, (  X,  1 ) )  *FMCDI  F  ( 
Hi ,  (T,  I )  )-AL*EP**2.0*FMCEXP( -ALAM- ANU  I  ♦FMCDIF  (HI,  (T,l),(X,l))» 
5.0E-I*AL*EP**2.0*FMCEXPf -ANU )*FMCDIF I H2, ( T, 2  ) I-5.0E- 1*EP*»2 .0*H2* 
X**(-2.0l*FMC$IN(Y)**(-l.0)*FMCC0S(Y)*FMCDlF(AL»(Y»l) 1-5. 0E-1*EP** 
2.0*H2*X**(-2.0)*FMC0IF(AL, ( Y , 2  I ) -X^ ( - 1 . 0 1 ♦ F MCE X P ( -ALAM ) ♦FMCD I F ( 
ALAM, ( X • 1 ) ) -2. 5E~l*FMCEXP( -ALAM) ♦FMCDIF (ALAM, (X,l > ) ♦F MCDI F ( ANU, C  X , 

1  I ) *2 .5E-1*FMCEXP( -ALAM ) ♦FMCD IF ( ANU, ( X, I ) l♦♦2.0♦5.0E- i*FMCEXP( - 
ALAM) ♦FMCDIF! ANU, (X,2) IS 

2U  IL 

-5.0E-l*EP**2»0*H0*X*+(-3.0)*FMCDIF(AL,(Y,I) ) *2 . 5E-1^EP*^2 . O^HO*X 
♦♦( -2.0) ♦FMCDIF ( AL,(Y,1 ) I ♦FMCD IF ( ANU, (X,1))*5.0E-1*EP**2. 0*H2*X** I 
-3.0)*FMCDlF(AL»(Y,l))*2. 5E-1 *EP**2 .0*H2*X** (-2.0) ♦FMCDIF ( AL • ( Y, 1 ) 

) ♦FMCDIF! ANU, (X,l ) ) -5.0E- I *E P**2 . 0*X *♦ ( -2 .0 ) *FMCE XP ( - ANU ) ♦FMCDI F ( 
AL, ( Y,l) ) ♦FMCDIF (HI, (T, 1 ) )-5.0E-I*EP**2.0*X**(-2.0)*FMCDIF(AX,(X,I 
) ) *FMC0 I F ( AL, ( Y, l ) ) *5 . OE- l ♦EP**2 . 0*X** ( -2 .0 ) ♦FMCD I F ( AL , ( Y , L ) I ♦ 
FMCDIF (HO, (X, 1 ) IS 
*U  11 

AL*EP**2.0*H1*X^(-1  ,OI*FMCEXP(-ALAM-ANU)*FMCDIF(  ALAM,  (  X,  1 )  l*AL*EP 
♦♦2.0*Hl*X**(-l.0)*FMCEXP( -AL AM-ANUI ♦FMCD IF ( ANU, (X,l) )»AL*EP**2.0* 
X**(-1.0I*FMCEXP( -ANU ) ♦FMCDIF ( AK , ( T , 1 ) )  - AL*E  P**2 • 0*X**( -1 . 0) ♦ 
FMCEXP(-ANU)*FMC0IF(H2,(T,1) )-5.0E-l*AL ♦EP^^2 .0*FMCEXP( -ANU) ♦ 

F MCDIF(AK,(T,1))*FMCDIF( ANU, ( X, 1 1 ) ♦AL*EP**2. 0*FMCEXP( -ANU) ♦FMCDIF! 
AKf(T,l),(X,lII*5.0E-l*EP**2.0*Hl*X^*l-2.0I^FMCSIN(YI**l-l.0l* 

F  MCCOS ( Y I *FMCEXP ( -ANU) *FMCOI F ( AL , ( Y , 1  I I^.OE-^EP^.O^HIO**!^  .0 
) ♦FMCEXPI -ANU ) *FMCD IF(AL,(Y,2))S 
1U  2L 

-5. OE-1 ♦EP**2. O^HO*X ♦* !-i.O)*FMCEXP( -ALAM) ♦FMCDI F(AL, (Y.l )  )♦ 

2 . 5E-1 ♦EP**2. 0*H0*FMCEXP( -ALAM )*FMCDIF(AL,(Y,1) ) ♦FMCDIF ( ANU, ( X, 1 1 t 
♦  5 .0E-1*EP**2. 0*H2*X** ( -1 .0 1 *FMCEXP( -AL AM ) ♦FMCDI F ( AL , I V , 1  I )♦ 
2.5E-1*EP**2 •©♦H2*FMCEXP ( -ALAM ) ♦FMCDI F ( AL  » ( Y , 1 1 )* FMCD IF (ANU,(X,1) I 
-5 .0E- l*EP**2  «0*F MCEXP ( -ALAM )*FMCDIF(AK,(X,1 ) l*FMCOIF (  AL « ( Y, 1 )  )  + 

5 .OE-1 *EP**2 •0*FHCEXPI-ALAMl*FMCOIF( AL , ( Y , 1 ) I ♦FMCDI F ( HO. ( X , 1 1 )- 
5.0E-1*EP**2.0*FHCEXP(-ALAM-ANU)*FMCDIF( AL, ( Y,1 ) ) ♦FMCDI F ( HI . ( T , 1 ) I 
S 
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THE  RICCI  TENSOR  (continued) 


2U  2L 

-AKfAL*EP**2.0*X**(-2.0)-5.0E-i*AK*EP**2.C,>X**(-2.0)*FMCSlN( YJ**( 
-1.0) *FNCCOS ( Y ) *FHCDIF ( AL, ( Y, 1 )  l-5.0E-i*AK* CP**2. 0*X**( -2 .0)* 
FMCDIFtAL, (Y,2) ) ♦ AL*EP**2 . 0*H2*X** < -2.0 ) *FMCEXP < -ALAM )-5. OE- l *AL* 
EP**2.0*H2*X**<-1.0) ♦FMCEXP(-ALAH)*FMCD1F(ALAM, (X,l) )  «-5  .OE-1  *AL*EP 
**2.0*H2*X** (-1.0)*FMCEXP(-ALAM>*FMCDIF< ANU, ( X.l ) )-2. 0*AL*EP**2. 0* 
X**(-1.0)*FMCEXP(-ALAM)*FMC0IF( AK, (X. 1) 1 ♦5.0E-1*AL*EP**2.0*X** f 
-1.0J*FMCEXP(-ALAM1*FMCDIF(H0, IX, 1) ) »5. OE-1* AL*EP**2. 0*X** ( - 1 . 0) * 
FMCEXP(-ALAM)*FMCDIF(H2, (X,l ) )-AL*EP**2.0*X**(-1.0)*FMCEXP(-ALAM- 
ANU)*FMCDIF(Hl,(T,l> )  +  2 . 5E-l*Al*EP**2. 0*FMCEXP ( -ALAM) *FMCDI F  (  AK,  (X 
,  1 ) I*FMCDIF(ALAM»(X, 1 ) )-2 . 5E- l*AL*EP**2. 0*FMCEXP ( -ALAMI *FMCDI F < AK , 
( X  « 1 ) I *FMCDI F I ANU, ( X • 1 ) I -5 .0E-1*AL*EP**2 .0*F MCEXP ( -AL AM ) *FNCDI F ( AK 
,  <X,2))*5.0E-l*AL*EP**2.O*FMCEXPC-ANU)*FNC0IF{AK,{T,2) I^.OE-iaEP 
**2. 0*H0*X** ( -2.0 ) *FMCDI F I AL  * ( Y « 2 ) )-5.0E-l*EP**2.0*H2*X** ( -2 . 0 »♦ 
FMCDIF ( AL> (Y,2))*X**(-2.0)*FMCEXP(-ALAM)-X**(-2.0I-5.0E-1*X**<-1.0 
)*FMCEXP<-ALAM)*FMCD1F(ALAM, IX,in*5.0E-l*X**(-1.0)*FMCEXP(-ALAHI* 
FMCDIF  (ANU, (X. ID* 
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4U  21 

2.5r-l*fcP**2.0*Hl*FMCEXP(-ALAM-ANi  >  *FMCI)I  F  (  AL  , ( Y, 1 ) ) *FMCD I T ( ALAM ,  ( 
X  » 1 )  )  -2. 5E-1 *EP**2.0*H1*FMCEXP( -AL AM-ANU)*FMCDI  F  (  AL  * ( Y* 1 )  )*FMCDIF ( 
ANU, I  X , 1 ) l-5.0E-l*EP**2.0*FMCEXP(-ALAM-ANU)*FMCDIF( AL,(Vf 1 )  )* 
FMCOIF  (HI , (X,l ) ) *5.0E-1*EP**2.0*FMCEXP(-ANU)*FMCDIF( AK,  (T,  l )  I* 
FMCDIF (AL, (Y, 1 ) ) *5 .0 E-1*EP**2 .0*FMCEXP ( -ANU) *FMC0  IF ( AL , (Y,  l)  I* 
FMCDIF1H2, (T,l))t 
3U  3L 

-AK*AL*EP**2.C*X**(-2.0)-5.0E-l*AK*EP**2.0*X**(-2.0)*FMCSIN(Y)**( 
-1.0I*FMCCOS  <Y)*FMCOIF(AL, ( Y , l >  I -5 . OE-1 *AK*E P**2 . 0*X**( -2 .0 ) * 
FMCDIF (AL, ( Y , 2 ) I ♦ AL*EP**2 . 0*H2*X** ( -2 .0 ) *FMCEXP  (  - AL AM )-5 . 0E-1*AL* 
EP**2.0*H2*X**(-i.0> *FMCEXP(-AIAM)*FMCDIF(ALAM, ( X, 1 1 ) *5 . OE-1 *AL*E P 
**2.0*H2*X**(-l.OI*FMCEXP(-ALAM)*FMCOIF(ANU,  (X,l I )-2.0*AL*EP**2.0* 
X**(-1.0)*FMCEXP(-ALAM)*FMCDIF(AK, ( X,X ) )*5.0E-1*AL*EP**2.0*X**( 
-1.0)*FMCEXP(-ALAM)*FMCDIF(H0, (X, 1  I  ) *5 .OE-1* AL*EP**2 . 0*X** ( -1 . 0 ) * 
FMCEXP  (-ALAM I *FMCD IF ( H2, (X,l ) )-AL*EP**2.0*X**(-1.0)*FMCEXP(-ALAM- 
ANU)*FMCDIF( HI , ( T, 1 ) ) *2 . 5E-1 *AL*EP**2 . 0*FMCE XP ( - ALAM I *FMCD IF ( AK, ( X 
, 1 l)*FMCDIF<  ALAM, (X, 1) ) -2 . 5t-l*AL*EP**2 . 0*FMCEXP ( -ALAM) *FMCOI F ( AK, 

( X,l» )*FMCDIF( ANU, (X,l ) )-5.0E-l*AL*EP**2.0*FMCEXP(-ALAM)*FMCDIF( AK 
,  (X,2)  )+5.OE-l*AL*EP**2.0*FMCFXP(-ANli)*rMCOIF(AK, (T,2) >+5.0E-l*EP 
**2.0*H0*X**(-2.0)*FMCSl',(Y)**(-l.O)*FMCCOS(Y)*FMCDlF(AL»  I Y,l  )  )- 
5.0E-1*EP**2.0*H2*X**(-2.0)*FMCSIN«Y)**(-1.0)*FMCCOS(Y)*FMCDIF(AL, 

( Y,1 ) )*X**(-2.0)*FMCEXP(-ALAM)-X**(-2.0)-5.0E-1*X**(-1.0)*FMCEXP(- 
ALAM)*FMCDIF (ALAM, (X,l) )*5. 0E-1*X**(-1. 0)*FMCEXP( -ALAM) *FMCDIF (ANU 
, (X,l)  )t 
11)  4L 

-AL*EP**2.0*X**(-1.0)*FMCEXP(  -  AL AM ) *FMCDI F ( AK ,  (  T,  1 ) )*AL*EP**2.0*X 
**1-1.0)*FMCEXP(-ALAM)*FMCOIF(H2,(T,1)I*5.0E-1*AL*EP**2.0*FMCEXP(- 
ALAM)*FMCDIF ( AK, (T,l  )  ) *FMCDI F I ANU , (X,l ( ) -AL*EP**2.0*FMCEXP(-ALAM) * 
FMCOIF (AK, (T , l ) , (X,l ) )-5.0E-l*EP**2.0*Hl*X**(-2.0)*FMCSIN(Y)**( 
-l.J)*FMCCOS( Y)*FMCEXP(-ALAM)*FMC01F( AL, (Y, 1 ) I -5. 0E-1*E P**2.0*H1 *  X 
**(-2.0)*FMCEXP(-ALAM)*FMC0IF(AL,< Y,2) It 
2U  4  L 

-2.5E-1 *EP**2.0*H1*X**( -2.0 )*FMCEXP( -ALAM )*F MCOI F ( AL, (Y,l ) )*FMCDIF 
(ALAM,  (X, 1 ) ) *2. 5E-1*EP**2. 0*HI*X** ( -2.0 )*FMCEXP(  -  ALAM )*FMCDIF ( AL, ( 
Y,l) )*FMCOIF( ANU, (X, L ) )  +  5. 0E-i*EP**2 . 0*X**( -2 .0  I *FMCE XP ( -ALAM) * 
FMCDIF (AL, (Y , 1 ) )*FMCDIF(H1, (X, 1 ) ) -5. OE-1 *EP**2. 0*X**  ( -2 . 0 ) *FMCDI F  ( 
AK, (T, 1 ) ) *FMCDIF ( AL, (Y,l) ) -5. OE-1 *EP**2 . 0*X** ( -2 . 0 ) *FMCDI F ( AL , ( Y , 1 
l)*FMC0IF(K2,(T,lf)t 
4U  4L 

AL*EP**2.0*H2*X**(-1  . 0 ( *FMCEXP( -ALAM ) *FMCDI F ( ANU, (X, 1 ) ) -2 . 5E-1 *AL* 
E P**2.0*H2*FMCEXP( -ALAM) *FMCOIF( ALAM, (X,l) )*FMCDIF( ANU, (X,l) )♦ 
2.5E-1*AL*EP**2.0*H2*FMCEXP(-ALAM)*FMCDIFIANU,(X,1))**2.0+ 
5.0E-1*AL*EP**2.0*H2*FMCEXP(-ALAM)*FMCDIFI ANU,  (X,2 ) ) +AL*EP**2.0*X 
*♦(-1.0) *FMCEXP(-ALAM)*FMCDIF (HO, ( X , 1  I ) -2.0* AL*EP**2 . 0*X** (- 1 . 0) * 
FMCEXP (-ALAM-ANU)*FMCOIF(  HI, (T, 1 ) ) -5 . OE- 1*AL *EP**2 .0*FMCE XP( - ALAM ) 
*FMCDIF( AK, I  X, 1 ) ) *FMC0IF ( ANU, (X, 1 ) ) -2 . 5E-1*AL*EP**2. 0*FMCt XP (-AL AM 
)*FMCOIF( ALAM, (X, 1) ) *FMCO I F (  HO, ( X , 1 ) )*5.0E-1*AL*EP**2.0*FMCEXP(- 
ALArt)*FMCOIF( ANU, (X, 1) )*FMCDIF(HO, (X, 1) ) *2. 5E-1*AL*EP**2 . 0*FMCEXP( 
-ALAM )*FMCDIF( ANU, (X,l) ) *FMC0IF ( H2 , ( X , 1 ) » *5. OE-1* AL*EP**2 . 0*FMCEXP 
(-ALAM)*FMCOIFIHO, (X,2I I +5 .0E-1*AL*E  P**2 . 0*F  MCEXP  < -AL  AM-AN1I)  * 
FMCDIF (ALAM, (X,l ) )*FMC0IF(H1, ( T,l) ) -AL*EP**2 . 0*FMCEXP ( -AL AM-ANU) * 
FMCDIFIH1, (T ,1), (X,l ))*AL*EP**2.0*FMCEXP(-ANU)*FMCDIF(AK,( T,2) )♦ 
5.0E-1*AL*EP**2 . 0*FMCEXP( - ANU I *FMCDIF ( H2» (T, 2 ) ) ♦ 5 .OE- 1*EP**2 • 0*H0* 
X**(-2.0)*FMCSIN(Y)**(-1.0)*FMCCOS(Y)*FMCDIF(AL, ( Y, 1 » ) *5. 0E-1*EP** 
2.0*HO*X**(-2.0I*FMCDIF(AL,(Y,2))*X**(-1.0)*FMCEXP(-ALAM)*FMCOIF( 
ANU, ( X , 1 ) ) -2 . 5E-1*FMCEXP( -AL AM) *FMCOIF ( ALAM, ( X, 1 ) ) *FMCDI F ( ANU, I X , 1 
) )+2.5E-l*FMCEXP( -AL AM)*FMCDIFI ANU, ( X, 1 ) ) **2 . 0*5 . OE-l*FMCEXPl -ALAM 
) *FMCD IF ( ANU , I  X  ,  2  )  )  t 
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THE  RICCI  SCALAR 

-2 • 0* AK* AL*EP**2 •  0*X **l -2. 0 ) -AK*FP**2 • 0*X** ( -2 .U  ) *FMCS I N <  Y I ** I -l . 0 

) ♦  FMCCOSI Y)^FMCDIF<  AL, ( Y,  1  )  )-AK*EP**2 .0*X**(-2.0)*FMCDIF Ml  t ( Y(2 ) ) 

♦2.L*AL*tP**2.0*H2*X**l-2.0)*FMCEXP(-ALAM|-2.0*AL*EP**2.PlH2*X**( 
-1.0  (♦FMCEXPl -ALAM I ♦FMCDIF (ALAM, IX, 1 ) > ♦ 2 . 0* AL*E P**2 . 0*H2* X** t - 1. 0  I 
♦FMCEXPl -ALAM  |  ♦FMCD I  F  l  ANU,  (  X,  11  ) -5 . 0E-  1*  AL*t: P**2 . 0*H2*FMCE  XPt-AL  AM 
(♦FMCDIF! ALAM, (X, 1 )) ♦FMCDIF (ANU,(X,1 I ) ♦P.0E-1*AL+EP+*2.C*H<^FMCEXP 

( -ALAM )  *FMCD I F ( ANU, ( X, 1) ) **2 . 0+AL*EP* *2 . 0*H2* FMCEXPl -ALAM ) .-FMCD!  F ( 
ANU, IX, 2 ))-6.C* AL*EP**2.3*X**t-l.O) ♦FMCEXPl- ALAM  I  ♦FMCDIF(AK, ( X,l>) 
♦Z.O*AL*EP**2.0*X**I-1.0»*FMCFXPl-ALAM)*FMCOIF(HO, ( X , 1 ) )  ♦  2  . 0*AL*E  P 
**2.0*X**(-l.0)*FMCEXP(-ALAM)*FMCDIF(H2,IX,l > ) -4. 0*AL *E P**u  . 0*X** ( 
- 1 • U ) ♦FMCfcXP ( -ALAM— ANU) ♦FMCDIF ( HI  , (T , 1 ) ) +AL* E P**2 . 0*FMCE X P ( - AL AM ) * 
FMCDIF(AK,(X,l))*FMC0IF(ALAM,(X,in-AL*CP**2. 0*FMCEXP  (  -  AL  AM  * 
FMCDIF (AK, (X, 1  )  )*FMCDIF(  ANU,  IX,  1)  )~2  .  C*  AL*EP  **2 . 0*FMCEX  P  (  -AL  AM  )  * 
FMCD IF (AK,(X, 211-5.0. -1*AL*EP**2.0*FMCFXP(-ALAM  I  ♦FMCD IF ( AL AM, ( X, 1 ) 

) *FMCD I F ( HO , (X.l) )*AL*EP**2.0*FMCEXPI-ALAM)*FMCOI F( ANU, (X, ! ) )  ♦ 

FMCD I F I  HO, I  X , 1 ) ) *5. 0C-1*AL*EP**2.0*FMCE  XP ( -ALAM I ♦FMCDIF ( A Nil, I  X , 1  I ) 
♦FMCDIF ( H2 , ( X , 1  I )+AL*EP**2.0*FMCtXP<-ALAM)*FMCDlF (HO, <X,2 )  ) ♦ AL*EP 
*♦2 . O^FMCEXP ( -ALA  "'-ANU) *FMCDI F( ALAM, (X,1))*FMCDIF(H1,(T,1)  1-2. 0*AL 
♦  EP**2 .0*FMCtXP(-ALAh-ANU  I  ♦FKCDIFl HI  ,(T,11*(X,1)  l+2.0*AL*E',<,*2.0<' 
FMCCXP ( - ANU) ♦FMCD I F ( AK,(T,2))*AL*EP**2.0*FMCfcXPl -ANU) ♦FMCD I F ( H2 , IT 
,  2  )  )  *£P**2.0*HO*X**  (-2.0)*FMCSIN(Y)**(-1 .  Cl)  ♦FMCCOSI  Y )  ♦FMCD  I  F  (  AL ,  (  Y 
,l))+-EP**2.0*H0*X**l-2.O)*FMCDIF(AL»(Y,2) )-E P**2 . 0*H2 *X ♦♦  ( -2 . 0 ) * 
FMCSlN(Y)+*(-l.G)*FMCCOS(Y)*FMCDIF(AL* (Y,l) )-EP**2.0*H2*X  +  M-  2.0)* 
FMCDIF (AL, (Y,2))*2.G*X*M-2.0)*FMCEXP(-ALAM)-2.0*X**(-2.0)-2.0*X** 
( -1.0 (♦FMCEXPl -AL AM) ♦FMCDIF ( ALAM, (X,l) ) *2.0*X** ( -1 .0  I *FMCE XP ( - AL AM 
)  *FMCD  IF  (  AND,  ( X , l ) I -5.0E-1*FKCEXP( -ALAM) ♦FMCDIF l ALAM, ( X , 1 )  I  *FMCD  I  F 
(ANU, (X,l) 1+0 . Or- 1 ♦FMCEXP( -AL AM) ♦FMCDI F ( ANU, (X,l I (♦♦2.0+FMCFXPl- 
ALAM) ♦FMCDIF ( ANU,  I  X , ? I ) t 
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THF  EINSTEIN  TENSOR,  MIXED  COMPONENTS,  0 1 AGONAL  ELEMENTS 
IU  1L 

AK*aL*EP**2.0*X**(-2.0)*5.0E-1*AK*EP<*2.0*X**(-2.0)*FMCSIN(Y)**( 

-1.0) *FflCCOS ( Y ) *FMCD IF ( AL ,(Y,]))+5.0E-l*AK*EP**2.O*X**(-2.Q)* 
FMCDIFIAL, (Y,2) > -AL* EP**2 . 0*H2*X **( -2 .0 ) *FMC E XP ( - AL AM  )  - AL *EP**2 . 0* 
H2*X**(-1.0)*fMCEXP(  -  AL AM ) *T MCOIF I ANU, ( X, 1 ) ) f AL *E P** 2 . 0*X** ( -I . 0 > * 
FMCEXP(-ALAM)*FNCDIF (AK, IX, 1 ) )-AL*EP**2.0*X** l-i.C)*FMCEXP(-ALAMI * 

FMCOIF(FIO,(X,l))  +  2.C*AL*HP**P,'j*X**(-1.0)*FMCEXPI-ALAM- AN  U)*F  MCOIF 
HU  , IT, I ) )*5.0E-1*AL*EP**?.G*FMCEXP(-ALAM)*FMC01F(AK, < X , 1 ) )*FMCDIF 
(ANU, (X,l) )-AL*EP**2.0*FMCEXP(-ANU)*FMCDIFIAK, (T,2) )- 5.0E-1 *EP** 
2.0*HO*X**(-2.0)*FMCSIN( Y)**(-1.0)*FMCC0S(YI *FMCOIF( AL, (Y,l J I- 
S.OE-l*EP**2.O*H0*X**(-2.O)*FMCniF<AL,(Y,2))-X**(-2.O)*FMCEXP(- 
ALAM)+X**(-2.Ci)-X**(-1.0)*FMCEXP(-ALAM)*FMCDIF(ANU,(X,l) l» 


5.0E-l*AL*EP**2.0*H2*X**(-I.0)*FMCEXP(-ALAM|*FMCDIF(ALAM,(X,ni- 
5.0E-l*AL*EP**2.0*H2*X**(-1.0)*FMCEXP(-ALAM) *FMCD I F < ANU , ( X , l  ) )♦ 
2.5E-1*AL*FP**2.0*H2*FMCEXP(-ALAMI*FMCDIF(ALAM, (X,l) I *FMCD I F  (  ANU , ( 
X, 1 ) )-2.5E-l*AL*EP**2.0*H2*FMCEXP(-ALAMI*FMCDIF(ANU,(X,ll)**2.0- 
5 . OE- 1 *AL  *CP**2 . Q*H2*FMCEXP( -ALAM) *FMCDI F ( ANU , ( X , 2 ) ) *AL*EP  **2 . 0*X 
**(-1.0)*FMCEXPt-ALAM)*FMCDIF(AK,|X,lJ)-5.0E-l*AL*EP**2.0*X**(-1.0 
)*FMCEXP(-ALAM)*FMCDIF(HO, ( X , l ) I -5. OE- 1*AL*E P**2 . 0*X** I -1 . 0 1  * 

FMCEXP ( -AL AM ) *FMCDIF |H2, I  X, 1 ) ) +AL*EP**2.0*X** | -1 .0)*FMCEXP(-ALAM- 
ANU)*FMCDIF(H1, IT, l > ) -2 . 5E- 1 *AL*£P**2 .0*FMCE XP ( -ALAM I *FHCD I F ( AK, (X 
, 1) )*FMCDIF( ALAM, IX, l) ) +2 . 5E-1 ♦ AL*EP**2 . 0*FMCEXP ( -AL AM) *FMCO I F ( AK , 
I  X  ,  1 ) ) *FMC DI F I ANU,  I  X,1 ) ) +5 . OE- 1  * AL*E P**2 . 0*F MCEXP I -AL AM  I *FMCl)  I F I AK 
, (X,2I  )+2.5E-l*AL*EP**2.0*FMCEXP|-ALAM)*FMC0lFIALAM, |X,1) )*FMCDIF( 
HO, (X, I) )-S.0E-l*AL*EP**2.0*FMCEXPI-ALAM)*FMCDIFI ANU, (X,l ) l*FMCOIF 
(HO, I  X , 1 ) 1-2.5E-1 *AL*EP**2.0*FMC EXP (-ALAM) *FMCDIF( ANU, ( X, 1 1 )* 
FMCDIF IH2 , (X , 1 ) ) -5.0E-l*AL*EP**2.0*FMCEXP(  -A  L AM ) *FMCD I F (HO,  I  X ,  2 ) » — 
5.0E-l*AL*EP**2.0*FMCEXP( -ALAM- ANU ) *FMCD IF  I  ALAM, (X,l) I*FMCDIF(H1 , I 
T.l) ) ♦AL*EP**2.0*FMCEXP(-ALAM-ANU)*FMCDIF(FU  ,(T,1),(X,1))- 
5.0E-1*AL*EP**2.0*FMCEXP(-ANU)*FMC0IF(AK,(T,2) ) -5 .OE- 1* AL*EP**2.  0* 
FMCE  XP ( -ANU) *FMCD IF  I  H2,  I  T , 2 ) ) -5 .OE- 1*EP**2. 0*H0*X** I -2. 0 1 *FMCS IN ( Y 
) **(-1.0)*FMCC()S(Y)*FMC0IF(AL, I Y, 1 )  )+5.0E-l*EP**2.0*H2*X**(-2.0)* 
FMCSIN(Y)**(-1.0)*FMCCOSCY)*FMCDIFIAL,(Y,I))+5.0E-L*X**(-1.0)* 

FMCEXP(-ALAM)*FMCDIF (ALAM, (X,l ) )-5.0E-I*X**(  -l.O)*FMCEXP(-ALAM)* 
FMCOI F ( ANU, (X, 1 ) I *2 . 5E- L*FMCEXP( - ALAM ) *FMCO I F ( AL AM, (X , I ) )*FMCDIF( 
ANU, ( X, 1 ) l-2.5E-l*FMCEXP( -ALAM) *FMCOIF( ANU, ( X, 1  I ) **2. 0-5.0E-1* 
FMCEXP (-ALAM) *FMCDIF (ANU, (X,2) )t 
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THE  EINSTEIN  TENSOR  (continued) 

3U  3L 

5.0E-1*AL*EP**2.0*H2*X**  (  - 1 .  (J )  *FMC.  E  XP  ( -AL  AM  )  *FMCD I F I  A L A M ,  (  X  ,  1  )  )  - 
5.0E-1*AL*EP**2.0*H2*X**(-1.0)*FM(,£XP(-ALAM)  ♦FMCD I  F  (  ANU  ,  (  X  ,  1  )  )  ♦ 

2  .bE-l*AL*EP**2.0*H2  *FMCEXP( -ALAM t *FMCDI F ( AL AM, t  X ,1 ) ) *FMCD IF { ANU , ( 
X, 1  I )-2.5E-l*AL*EP**2.0*H2*FMCEXP( - ALAM ) *FMC D  I F ( ANU,  IX,  1>)**2.0- 

S. OE-1*AL*EP**2.0*H2*FMCEXP(-ALAM)*FMCOIFIANU, (X,2))*AL*EP**2.0*X 
** (-1.0) *FMCEXP< - ALAM) *FMCOIF< AK, <  X, 1 ) )-5.0E-l*AL*EP**2.?*X**< -l .0 
) *FMCEXP(-ALAM)*FMCDlr {HO, (X,l ) !-5.0E-l*AL*EP**2.0*X**(-L .0)* 
FMCEXP1-ALAM)  ♦FMCDIF  (H2,  (  X,  1  I  )+AL*FP<=*2.0*X**(-i.0)*FMCf.XP(-ALAM- 
ANU)*FMCDIF( HI , ( r ,1 ) l-2.bE-l*AL*tP*<‘2.0<FMCEXP(-ALAM) *f MCOIFI AK, { X 
,  l  n* FMCD IF ( ALAM, (X, l ) ) ♦ 2 . 5E- 1 * AL * EP**2 . 0*FMC EXP  ( -ALAM) *FMCDIF(AK , 
(  X,I)  )  *F  MCO I F (  ANU  i  I  X  ,  I  I  )+b.0E-l*AL*EP**2.O*FMCEXP(-ALAM)*FMCDIF( AK 
.  (X,2  >  )+2.5E-l*AL<=EP**2.0*FMCEXP(-ALAMH  FMCD  IF  (ALAM,  (  X,  l  )  KFMCDI  F  ( 
HO,  (X,  1  )  l-5.0E-l*AL*EP**2.0*FMC£XP( -ALAM) *FMCDIF( ANU,  (X  ,  1  )  ( <‘F  MOD  I  f 
(  HO,  <  X  ,  1  )  )-2 . 5E-1*AL*EP**2.0*FMCF.XP(  -  AL  AM)*FMCDIF(  ANU,  (  X  ,  1  I  )  * 

FMCD IF (H2, (X,l ) ) - 5 . 0 E-l * AL *EP**2 . D*FMCE XP ( -A L AM ) *FMCD I F  (  HO , ( X , 2 )  )- 
5 .0E-1*AI.*EP**2 . 0*FMCEXP(-ALAM-ANU)*FMCDIF  (  ALAM,  (  X  ,  11  )  *  FMC  D I  F  (  HI  ,  ( 

T,  1  )  )  ♦AL*EP**2.0*FMCEXP(-ALAM-ANU)*FMCUIF(H1  ,(T,1),(X,1>I  — 
5.0E-1*AL*EP**2.0*FMCEXP(-ANU)*FMCDIF( AK, IT, 2  I  )- 5. OE - 1  *AL *E P**2 . 0* 
FMCEXP(-ANU)*FMCDIF(H2,  IT,  2)  ) -5  .  GE- 1  <•  EP**2. 0*H0*X**  ( -2.  0  I  *  FMCD  IF  ( 
AL, ( Y,2) >+5.0E-l*EP**2.0*H2*X**!-2.0 ) *FMCDI F ( AL,(Y,2))«-5.0E-1*X**( 
-1 .0  )*FMCEXP(-ALAM)*FMCDIF( ALAM,  (X, 1) )-5.0E-l*X**(-i.OI*FMCEXP(- 
ALAM)*FMCDIF( ANU,  (X,  1)  )  ♦  2 . 5E-  1*FHCF.X  P  I -ALAM )  *FMCD[F  I  ALAM,  (  X,  l  >  )* 
FMCDIFUNU, I  X, 1 ) )-2.5E-l*FMCEXP(-ALAM)*FMCD( F( ANU, (X, 1)  )**2.0- 
5.0E-1*FMCEXP(-ALAM) *FMCDIF( ANU, (X,2 )  It 

4U  4L 

AK*AL*EP**2.0*X<-*(-2.0)+5.0E-l*AK*EP**2.0*X**(-2.0)*FMCSIN(YI**( 
-1.0) *FMCCOS I Y) < FMCD  IF  I AL , I Y, 1 ) ) +5 .OE- 1  * AK *EP**2 . 0*X* *  I  -2 . 0 ) * 
FMC0IF(AL,(Y,2> ) - AL* EP**2 . 0*H2*X** ( -2 . 0 ( *FMC E XP < - AL AM ) ♦ AL *EP**2 . C* 
H2*X**(-1. 0)*FMCEXP(-ALAM)*FMC DIF (ALAM, (X,l )  ) +3. 0*AL *EP**  2 . 0*X** ( 
-1.0)*FMCEXP(-ALAM)*FMCCIF(AK, IX, 1 ) ) - AL *EP** 2 . 0*X** < -  1 . 0  I *FMCEXP ( - 
A LAM)* FMCD IF  I H2 , ( X , 1 ) ) -5 . OE- 1 *AL *E P* *2 . 0*FMC E XP ( -AL AM ) *FMCDI F ( AK , I 
X , 1 ) ) *FMCD I F ( ALAM, ( X ,1 ) ) * AL *EP**2 .0*FMC EXPI - A L AM ) *FMCDI  F l AK, ( X ,2 ) ) 
♦  5.0E-1*EP**2.0*H2*X**(-2.0)*FMCSIN<  Y)*< I -1 . 0  )  *FMCCQS I Y ) *FMCO I F I AL 
» I Y. 1 1 )+5.0E-l*EP**2.0*H2*X**< -2.0)*FMCDIF( AL  ,(Y,2))-X**(-2.0)* 
FMCEXP(-ALAM) *X** I -2 .0  I ♦ X** I -1 .0 ) *FMCEXP ( -AL AM)*FMCDIF( ALAM,  (X, l )  I 
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C .  ^  RIQtANN ,  the  Program  for  Calculating  Rlemann  Tens ora 
1.  Input  far  RUMANH 

RIEMANN  1b  the  program  which  calculates  the  covariant  components, 

Rahcd»  of  the  RieMimi  tensor.  The  Input  for  RTEMANN  i8  Identical  to  that 
for  EINSTEIN;  The  user  must  supply  ATjfalC  end  DEPEND  statements,  an  EPTERM 
declaration,  and  LET  definitions  of  the  metric  components  Gl(l,j)  and  C2(l,J). 
(See  5  III.B.l.)  A  listing  of  RIEMANN  showing  where  to  insert  the  required 
cards  follows. 
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oooooonn  o  00000000  o  00000000  000 


4  lll.c.2 


2.  Complete  Listing  of  MBMMW 


S I BFMC  TWO  NOOECK 

SYMARG 

RIEMANN,  A  PROGRAM  TO  CALCULATE  THE  RIEMANN  TENSCR 

ATOMIC  X,Y,Z,T,£P,ANU, ALAM  \  User:  Replace  these  cards 

DEPEND  I  ANUt  ALAM/X )  J  with  your  own  Input. 

ATOMIC  AND  DEPEND  STATEMENTS  MUST  BE  SUPPLIED  BY  USER 
IMMEDIATELY  PRECEC1NG  THIS  COMMENT. 

THE  ATOMIC  STATEMENT  MUST  LIST  ALL  FUNCTIONS  AND  INDEPENDENT 
VARIABLES  USED  IN  A  PARTICULAR  PROBLEM. 

THE  DEPtNO  STATEMENT  DEFINES  ALL  THE  FUNCIIGNAL  DEPENDENCIES 
FOR  THE  PROBLEM,  l.fc.  WHICH  FUNCTIONS  DEPEND  ON  WHICH  VARIABLES 

REAL  MCS ( 4  *  4, 4 1 
LOGICAL  U 

DIMENSION  CARD! 12),LINE(i2)«CCS(4»4,4),Gl(4,4),G2<4,4) 

INTEGER  EPTERM 

IF  EPTERM  =  0  TERMS  CONTAINING  ALL  POWERS  OF  EXPANSION 
PARAMETER  WILL  BE  RETAINED. 

IF  EPTERM  .NE.  0  TERMS  CONTAINING  THE  EXPANSION  PARAMETER 

TO  A  POWER  GREATER  THAN  2  WILL  BE  DISCARDED. 

A  CARD  CONTAINING  THE  INTEGER  VALUE  UF  EPTERM  MUST  FOLLOW  THIS 
COMMENT. 

.  User:  Replace  this  card 

EPTERM  *0  J  vith  your  own  input. 

00  1  I  =  1,4 
DO  l  J  =  1,4 
LET  G1CI.J)  =  0. 

LET  G2II.J)  =*  0. 

1  CONTINUE 

ALL  ELEMENTS  OF  THE  COVARIANT  METRIC  TENSOR,  G1II.JI,  AND  THE 
CONTRAV AR 1  ANT  METRIC  TENSOR,  G2(I,JI,  HAVE  BEEN  SET  EQUAL  TO 
ZERO.  . „„„ 

USER  MUST  SUPPLY  •  • LET  •  •  STATEMENTS  DEFINING  ALL  NON  ZtRO 
ELEMENTS  OF  THESE  TENSORS  IMMEDIATELY  FOLLOWING  THIS  COMMENT. 


LET  GIll.ll 
LET  G 1 1  2 , 2 ) 
LET  G1I3.3I 
LET  G1 1 4,4 ) 
LET  G2  < 1 , 1 ) 
LET  G2I2.2I 
LET  G2 ( 3 , 3 ) 
LET  02(4,41 


-FMCEXPI  ALAMJ 
-X**2 


-(X*FMCSIN(Y )  1**2 
FMCEXPI ANU 1  ) 

-FMCEXPI-ALAM) 

-L/X**2 

-1 / ( X*FMCS I N I Y  1  1**2 
FMCEXP ( -ANU  1  J 


User:  Replace  these  cards 
with  your  own  input. 
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c 

c  the  covariant  metric  tensor  --  OUTPUT 
c 

501  F0RMAT(651H1C0VAR1ANT  metric  tensor  non  zero  ELEMENTS  I  INPUT)  ) 
DO  2  I  =  1.4 
00  2  J  *  I  » 4 

LET  Q  =  MATCH  ID*  GllltJI.O. 

IFIQ)  GO  TO  2 
MR  I TE ( 6 (800 ) I • J 

800  FORMAT l LH  I5.1HL  I5.1HL  ) 

BEGIN  =  0. 

46  LET  BEGIN  *  BCOCON  G1 (  I  * J ) , L I NE , 12 

WRI TE (6,200 ) (LINEIL) ,L=2, 12) 

IFIBEGIN.NE.O. )G0  TO  46 
2  CONTINUE 
C 

C  THE  CONTRAVARIANT  METRIC  TENSOR  —  OUTPUT 

C 

700  FORMAT^ 54H0C0NTRAVAR I  ANT  METRIC  TENSOR  NON  ZERO  ELEMENTS  (INPUT)) 
DO  47  1=1.4 
00  47  J  *  1.4 

LET  Q  =  MATCH  ID.  G2(I.J)*0. 

I F ( Q ) GO  TO  47 
MR  ITE ( 6. 500 )  I.J 
500  FORMAT ( 1H  I5.1HU  I5.1HU  ) 

BEGIN  =  0. 

48  LET  BEGIN  =  BCDCON  G2 ( I . J • . L INE, 12 
MRITE16, 200)1 LINEIL) ,L=2,12> 

IFIBEGIN.NE.O. )G0  TO  48 

47  CONTINUE 
C 

C  THE  COVARIANT  CHRI STOFFEL  SYMBOLS 

C 

MR  I TE 1 6 . 201 )  c  . 

201  FORMAT I34H1THE  COVARIANT  CHRISTOFFEL  SYMBOLS  ) 

DO  3  I  =  1.4 

DO  3  J  =  1.4 

DO  3  K  =  1.4 

LET  T 1  =  0. 

GO  TO  (31,32.70,71),  I 

31  LET  Tl  =  FMCDIFIG1I J.KI.X.l) 

GO  TO  33 

32  LET  Tl  =  FMCDIFIG1I  J.lO.Y.l) 

GO  TO  33 

70  LET  Tl  =  FMC0IFIG1I J,K),Z,1) 

GO  TO  33 

71  LET  Tl  =  FMCDIFIG1I J,K),T,1) 

33  CONTINUE 
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}  III.C.2 


LET  T2  =  0. 

GO  T0(35,36,72,73)»J 

35  LET  T 2  *  FMCDI F ( Gl ( I , K ) ,X, l ) 

GO  TO  34 

36  LET  T2  *  FMCOI F( Gl ( I ,K ) » Y» 1 ) 

GO  TO  34 

72  LET  T2  =  FMCOI F ( G1 C I , K > , l , 1 ) 

GO  TO  34 

73  LET  T2  «  FMCOI F I Gl (I . K  ) , T , 1 ) 

34  CONTINUE 

LET  T 3  =  0. 

GO  TO  < 38, 39,74.75  )  ,K 
38  LET  T3  =  FMCOI F(  G1 ( I . J  1 1 X , i  ) 

GO  TO  37 

33  LET  T3  *  FMCDI F ( Gl I  I  * J > . V » 1 ) 

GO  TO  37 

74  LET  T3  =  FMCDI F ( Gl ( I , J  ) » Z , 1 > 

GO  TO  37 

75  LET  T3  =  FMCOIFI Gl ( I » J  )  *T » l ) 

37  LET  CCSU.J.K)  =  0.5*IT1+T2-T3> 

LET  CCS(J.I.K)  =  CCSU.J.K) 

LET  Q  =  MATCH  ID»  CCSU.J.K), 0. 

I  Ft  U ) GO  TO  3 

WRITEI6, 203)1  . J.K 

203  FORMAT  1 1H  15.  2HLS  I3.2HLS  I4.1HL) 
LET  ANS  =  EXPAND  CCS(I.J.K) 

LET  ANS  =  ORDER  ANS.  INC.FUL 
BEGIN  =  0. 

5  LET  BEGIN  =  BCDCON  ANS, LINE, 12 
WRITE ( 6,200 M L INE( L) ,L*2» 12) 

200  FORMAT ( IH  5X  12A6) 

IFIBEGIN.NE.O. )G0  TO  5 
ERASE  ANS 
3  CONTINUE 

ERASE  T1.T2.T3 

THE  MIXED  CHR l STOFFEL  SYMBOLS 


205  FORMAT (  30H1 THE  MIXED  CHRI STOFFEL  SYMBOLS  ) 
00  6  I  -  1.4 
DO  6  J  =  1,4 
DO  6  K  =  1.4 
LET  MCS(l.J.K)  =  0. 

DO  7  L  =  1.4 

LET  TT  =  G2(K,L)*CCS( I ,J»L) 

LET  TT  =  EXPAND  TT 
IFIEPTERM1506, 507,506 
507  LET  SUM  *  TT 
ERASE  TT 
GO  TO  514 
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506  CONTINUE 

LET  SI  *  COEFF  TT,EP**0,R 

LET  S2  “  COEFF  TT,EP**1,R 

LET  S3  *  COEFF  TT,EP**2,R 

LET  SUM  «  SI  ♦  S2*EP  ♦  S3*EP*EP 

ERASE  TT 
ERASE  SI t  S2 1 S3 
514  CONTINUE 

LET  MCS ( I  *  J  »K )  =  MCS ( I «  J , K )  ♦  SUM 
ERASE  SUM 

7  CONTINUE 

LET  MCS(I.JiK)  =  EXPAND  MCSII.J.KJ 
LET  MCS I J 1 1 |K )  =  MCS ( I  *  J »  K ) 

LET  0  =  MATCH  IDt  MCS (  I , J . K ) , 0. 

I F I Q ) GO  TO  6 
WR I TE ( 6, 769 )  I,J,K 

769  FORMAT ( 1H  I5.2HLS  I3.2HLS  14.1HU  ) 
LET  ANS  *  ORDER  MCSU , J,K I , INC.FUL 
BEGIN  =  0. 

8  LET  BEGIN  *  BCOCON  ANS, LINE, 12 

WRITE (6,200)  (LINE(L),L*2, 12) 
IFIBEGIN.NE.O. )G0  TO  8 
ERASE  ANS 
6  CONTINUE 


THE  RIEMANN  TENSOR 


WRITE(6,215) 

215  FORMAT (41H1THE  RIEMANN  TENSOR,  COVARIANT  COMPONENTS  ) 

DO  15  1  *  1,4 

ll  =  in 

iF(LL.GT.4)LL=4 
DO  15  J  =LL  ,4 
DO  17  K  =  1,4 
DO  17  L  =  J ,  4 
LET  SUM  =  0. 

DO  18  MU  =  1,4 

LET  TT  =  ( MCSI J ,K ,MU>  *CCS( I ,L« MU )-MC S ( J , L , MU ) *CCS 1 1 ,K,MU) 

1  ) 

LET  TT  =  EXPAND  TT 
IF(EPTERM)520,521,520 
521  LET  SUM  =■  SUM  +  TT 
GO  TO  18 
520  CONTINUE 

LET  SI  *  COEFF  TT,EP**0,R 

LET  S2  -  COEFF  TT,EP**1,R 

LET  S3  =  COEFF  TT,EP«"*2,R 

ERASE  TT 

LET  SUM  =  SUM  ♦  SI  ♦  S2*EP  ♦  S3*EP*EP 
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ERASE  S1,S2,S3 
18  CONTINUE 

LET  Ul  *  0. 


GO 

TO 

(41, 42, 80, 81), K 

41 

LET 

01 

*  FMCOIFICCSI J,L, I  I ,X, 1) 

GO 

TO 

40 

42 

LET 

01 

*  FMCOIF(CCS(J,L,I) ,Y,1> 

GO 

TO 

40 

80 

LET 

01 

=  FMCDIF(CCS(J,L,l),Z,i) 

GO 

TO 

40 

81 

LET 

01 

=  FMCOIF(CCS< J,L,I),T,1) 

40 

LET 

02 

x  0. 

GO 

TO 

(44,45,82,83) ,L 

44 

LET 

02 

*  FMCDIF(CCS(J,K,I) , X , 1 ) 

GO 

TO 

43 

45 

LET 

02 

*  FMCOIF(CCS(J,K,ll,y,l) 

GO 

TO 

43 

82 

LET 

02 

=  FMC0IF(CCS(J,K,1),Z,1) 

GO 

TO 

43 

83 

LET 

02 

=  FMCDIF(CCS( J.K.I ) ,T,1) 

A3  CONTINUE 

LET  DIFF  *  Dl-02 

ERASE  01,02 

LET  RMT  =  DIFF  ♦  SUM 

ERASE  SUM.OIFF 

LET  ANS  *  EXPANO  RMT 

ERASE  RMT 

LET  Q  *  MATCH  10,  AN$,0. 

IFIQIGO  TO  17 
WRITEI6, 206)1, J,K,L 

206  FORMAT ( 1H  5X  I5.1HL  I5.1HL  I5.1HL  15,IHL  ) 
LET  ANS  =  ORDER  ANS.INC.FUL 
BEGIN  »  0. 

25  LET  BEGIN  *  BCDCON  ANS, LINE, 12 
WRI TEI6.200) (LINE! IL) , IL>2, 12) 
IFIBEGIN.NE.O. )G0  TO  25 
ERASE  ANS 
17  CONTINUE 
15  CONTINUE 
STOP 
ENO 
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3.  Output  far  RIBTAHN 


The  printed  output  for  RBMANN  includes:  the  metric  components,  end 
g1"1;  the  Chrlstoffel  symbols  (eqs.  (2)  end  (3))}  end  the  covmrlent  components 
of  the  Rleomnn  curvature  tensor 


Rabed  “  rbc  rsdl  "  Fbd  r»d  +  rbda,c  “  rbca,d  *  '8' 

The  format  Is  similar  to  that  for  the  program  EINSTEIN,  and  error  messages  are 
Identical  to  those  for  EINSTEIN.  The  computer  prints  only  the  nonzero  ccopo- 
nents  of  the  Rlemann  tensor;  and  of  those  components  which  can  he  obtained  from 
each  other  by  symmetry  operations,  It  prints  only  one.  For  example,  R^^  la 

printed  if  it  Is  nonzero;  but  Rgljl*'  R12h3»  R2l43»  R34l2»  R4312»  R3421»  *nd 
1^312  are  never  printed  because 

R1234  *  "  R2134  *  "  R1243  *  +  H2l43  ”  +  B34l2  “  '  *  “  R3421  “  +  B4312* 


As  an  example  of  output  from  a  computation  by  RIBUNN,  we  give  below  the 
results  for  the  static,  spherically-synaetric  metric  of  equation  (l).  The 
Input  vhlch  generated  these  results  was  given  in  5  IU.B.l  and  also  In  the 
listing  of  RUMANN  In  }  III.C.2. 

METRIC  TENSOR  COMPONENTS  AND  CHRISTOFFEL  SYMBOIfl 
-  output  Identical  to  that  from  EINSTEIN  as  given  on  page  23. 


(HE  MEMANN  TENSOR.  CCVARI  ANT  COMPONENTS 
1L  2L  U  20 

-5.0E-1*X*FMCC1F(ALAH, IX • 1 ) )< 

1L  30  U  3L 

-5.0E-1*X*FMCSIN(Y  » **2 .C»  FFCC l F ( ALAR, <X,l » )  f 
1L  AO  10  AO 

2 . 5E- 1*FFCEXP( ANO I * FFCC  IF ( ALAN, I X. 1 ) l*FHCDIFI ANlt. (X.1II-2.5E-1* 
FMCEXP(ANU*FMCCIF(ANO,  CX,inee2.O-5.0E-ieFMCEXP(AN0»AFMCCIFtANO.  ( 
X.2IIS 

2L  3  L  20  30 

X**2.C*FHCSIN(Y)**2.0*FFCEXP(-A0AM)-X**2.Q*FPCSlN(Y)ee2.0l 
20  AO  20  A0 

-5.0E-l*xeFMCEXP<-AOAM  +  ANUaFMCClFUNU,«X,H)$ 

30  AO  30  A0 

-5.0E-l*X*FFCSIN(Y)A*2.G*FPCEXP(-AOAMANU)eFPCCIF(ANO»JX.l)l$ 
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Bennie  Prhbla—  for  RIEMAHH 

Input  and  output  for  on*  saaple  problem,  the  static  spherically  eymnetric 
■strlc  of  equation  (l),  vers  given  In  ths  last  tvo  sections.  As  a  second 
enable  consider  the  metric  (7)  for  a  star  in  even-parity,  nonradial  pulsation. 
The  input  to  RXSMAHK  far  this  metric  is  identical  to  the  input  to  EIKSTSIN 
(see  page  25  ).  The  output  frcei  RIB4A1W  is  partially  identical  to  that  tram 
khowhth  (pages  26  -  28  are  identical))  hut  instead  of  outputting  Ricci  and 
Einstein  tensors  aft'ir  the  Christoff  el  systole,  RXEMAKK  outputs  Rlenann  tensor 
components,  of  which  we  present  only  five: 


THE  R1EMANN  TENSOR,  COVARIANT  COMPONENTS 

5.0E-i*AK*AL*EP**2.0*XeFMCDlFtALAM,IX,l»l-AL*EP**2.0*XeFMC0IFIAK,I 
X,in+5.0E-l*AL*fcpee2.0*X*FMCDIFIH2,U,l>  )*2.5E-1*ALSEP**2.0*X** 
2.0*FMCUIF(AK,(X,1)I*FMC0IFIALAM,(X,1I )-5.0E-l*AL*EP**2.0exe*2.0* 

FMCDIF(AK,(X,2) I -5.0E- l*EPe*2.0*H2*FMCEXP (ALAM)*FMC0IF(AL,  CY,2)  )~ 
5.0E-1*X*FMCDIFIALAM,(X,1) It 

1L  2L  1L  4L 

-5.0E-ieEP**2.0*Hl*X**(-1.0I*FMCOIFI AL, IY, 1 1 )-2.5E-l*EP**2»0*Hl* 
FMCD1FI AL.IY.1I )*FMCOIF< ALAN, C X, 1 1 1 *2. 5E-1*EP**2 . 0*H1*FMCD I F I AL  ,  (Y 
, HjeFMCDIFIANU, C X, 1 1  I -5.0E-l*EP*e2.0*FMCEXP I AL AM  I *FMCDIF l  AL.CY,  1) 

1*FMC0IF(H2. <T,l))*5.0E-l*EP**2.0*FMCDlFI  AL.IY.ll I*FMCDIFIH1,IX,1) 

It 


5.0E-l*AL*EP**2.0*Hl*X*FMCEXPI-ALAM)eFMCDIFI ANU, (X,l) •♦5.0E-1*AL* 

EP**2.0*X*FMCDIF(AK» ( T ,1 ) l-5.0E-l*AL*EP**2.0*X*FMCCIF(H2, 1  T,1 1 I- 
2. 5E-1*AL*EP**2.0*X**2.0*FMCDIFI AK,  I  T, II )*FMC0IFC  ANU, ( X, 1 1  !♦ 
5.0E-l*AL*EP**2.0*Xe*2 .0*FMC0IFC  AK,(T,1),(X,1) I ♦5.0E-l*EP**2 .0*H1* 
FMC0IF(AL,(Y,2llt 

-5.OE-l*AL*EP*»2!o*Hl*X*FMCEXPI-ALAM)*FMCOIFUNU,(X,lI)-5.0E-l*AL* 

EP**2 . 0*X*FMCD IF  I AK , (T , U ) *5. OE-l*AL*E P**2. 0*X*FMC0 1 F ( H2»  I  T,  1 1 )♦ 
2.5E-l*AL*EP**2.0*xe*2.0*FMCDlF(AK»IT,l)l*l  MCO IF  C  ANU,  I  X,  1 1  1- 
5.0E-1*AL*EP**2.0*X**2 .0*FMCDl F ( AK,  I T, 1 ) , (X , 1 1 )-5.0E-l*EP**2.0*Hl* 
FMCDIF ( AL , ( Y,2 ) )  t 

5.0E-l*AK*AL*EP**2.0*X*FMCSINm**2.0*FMCDIFIALAM,  <X,l)  »-AL*EP** 
2.0*X*FMCSIN(YI**2.0eFMCDIF(AK»lX,l) I *5.0E- 1 *AL *EP**2. 0*X*FMCS IN l Y 
)**2.0*FMCDIF(H2,IX,nK2.5E-l*AL*EP**2.0*xe*2.0*FMCSIN(YIPe2.0* 

FMCDIF < AK, <X,1) »*FMCDIFI ALAM, I X, II !-5.0E-l*ALeEP*e2.0*X**2 . 0* 
FMCSINI Y I **2.0*FMC01  F I AK, I X, 2  I  1-5.0E-1*EP**2 . 0*H2*FMCS INI  Y ) *FMCC0S 

I YI*FMCEXPIALAMI*FMCDIF(AL,IYtll l-5.0E-l*X*FMCSINCYI**2.0*FMCDIF( 

ALAM, (X, lilt 


hi 


§  m.D.i 


D. 


MOTION,  the  Program  for  Calculating  the  Divergence 


of  the  Stress-] 


1»  Input  for  MOTION 

MOTION  Is  the  program  which  calculates  the  divergence  of  the  stress-energy 
tensor,  W#b  The  input  far  MOTION  is  identical  to  that  for  EINSTEIN 
(page  9)  except  for  this;  One  must  also  supply 

6.  "LET"  statements  defining  all  nonzero  mixed  components  of  the 
stress-energy  tensor,  W^,  in  terms  of  atomic  variables. 
is  denoted  in  MOTION  by  V(l,j)  -  i.e.,  the  first  index  is  co¬ 
variant;  the  second  is  contra  variant . 

A  listing  of  MOTION  showing  where  to  insert  the  required  cards  follows: 


k2 


oooooooo  o  oonoooon  r»  onooooooo  r>  o  o  o  o  o 


i  m.D.2 


2,  Complete  Listing  of  MOTIOH 
IBFMC  EOH  NODECK 

MOTION*  A  PROGRAM  TO  CALCULATE  THE  DIVERGENCE  OF  THE 
STRESS-ENERGY  TENSOR 


SVMARG 

ATOMIC  X,Y,Z,T,EP»ANU,ALAM,P,RHO  \  User:  Replace  these  cards 
OEPENO  (ALAM,ANU,P,RHO/XI  f  with  your  own  Input. 


ATOMIC  AND  OEPENO  STATEMENTS  MUST  BE  SUPPLIEO  BY  USER 
IMMEDIATELY  PRECEDING  THIS  COMMENT. 

THE  ATUMIC  STATEMENT  MUST  LIST  ALL  FUNCTIONS  AND  INDEPENDENT 
VARIABLES  USED  IN  A  PARTICULAR  PROBLEM. 

THE  DEPEND  STATEMENT  DEFINES  ALL  THE  FUNCTIONAL  OEPENOENCIES 
FOR  THE  PROBLEM,  I.E.  WHICH  FUNCTIONS  OEPENO  ON  WHICH  VARIABLES. 

REAL  MCS I  4,4,41 
LOGICAL  Q 

D I  MENS  ION  CARD!  12  I , L I NE ( 1 2 ) ,CCS I  4 , A ,4 I ,G1 1  4 ,  A 1  ,G2 14 ,4 » 

DIMENSION  W(4,4I 

INTEGER  EPTERM 

IF  EPTERM  =  0  TERMS  CONTAINING  ALL  POWERS  OF  EXPANSION 
PARAMETER  WILL  BE  RETAINED. 

IF  EPTERM  .NE.  0  TERMS  CONTAINING  THE  EXPANSION  PARAMETER 

TO  A  POWER  GREATER  THAN  2  WILL  BE  DISCARDED. 

A  CARD  CONTAINING  THE  INTEGER  VALUE  OF  EPTERM  MUST  FOLLCW 
THIS  COMMENT. 


EPTERM  =  0 

DO  1  I  =  1,4 
DO  1  J  *  1,4 
LET  GKI.JI  =  0. 
LET  G2 (  I ,  J  )  =  0. 

1  CONTINUE 


}User:  Replace  tbla  card 
vlth  your  own  input. 


ALL  ELEMENTS  OF  THE  COVARIANT  METRIC  TENSOR,  Gl(I.J),  AND  THE 
CONTRAVAR I  ANT  METRIC  TENSOR,  G2(I,J),  HAVE  BEEN  SET  EQUAL 
TO  ZERO. 

USER  MUST  SUPPLY  "LET"  STATEMENTS  DEFINING  ALL  NON  ZERO 
ELEMENTS  OF  THESE  TENSORS  IMMEDIATELY  FOLLOWING  THIS  COMMENT. 


LET  Glll.ll 
LET  01(2,21 
LET  01(3,11 


-FMCEXPI ALAMI 
-X**2 

-(X*FMCSIN(Y) 1**2 


LET  01(4,41 
LET  G2I1.1I 
LET  02(2,2) 
LET  G2  (  3,  3 1 
LET  G2 1  4, 4 ) 


FMCEXP ( ANU I 
-FMCEXPI -ALAMI 
- 1 /X**2 

-1/(X*FMCSIN(Y» )**2 
FMCEXP ( -ANU I 


l  User:  Replace  these  cards 
[  vith  your  own  Input. 


J 


1*3 


o  o  o  o  o  o  o 


}  III.D.2 


WR ITE ( 6, 50 1 ) 

501  FORMAT  1  5 1H  1C0V  AR  I  AN  T  METRIC  TENSOR  NCN  ZERU  ELEMENTS  I  INPUT  )  I 
DO  2  I  =  1,4 
DO  ?  J  *  1,4 

LET  0  =  MATCH  ID,  Gl(I,J),0. 

IF(Q)  GO  TO  2 
WRITE(6, 800)1, J 

800  FORMAT (  1H  [5,1HL  I5,1HL  ) 

BEGIN  =  0. 

40  LET  BEGIN  =  8CDCUN  G 1 ( I , J ) , LINE, 12 
WRITE ( 6, 2001 1L INEIL ) ,L=2,12) 

IFIBEGIN.NE.O. )G0  TO  40 
2  CONTINUE 

WK I TE I  6 , 700 ) 

700  FORMAT (54H0C0NTRAVAR I  ANT  METRIC  TENSUR  NON  ZERO  ELEMENTS  (INPUT)) 
DO  42  I  =  1,4 
DO  42  J  =  1 ,4 


LET  Q  =  MATCH  ID, 
IFIQIGO  TO  42 

G2< l,J),0. 

WR  I  TE 1 6 , 500 )  I , J 

500 

FORMAT ( 1H  I  5 , 1HU  1 
BEGIN  =  0. 

5,  1HU  1 

41 

LET  BEGIN  =  BCDCON 

G2  (  I , J ) , LINE, 12 

WK I T  E ( 6 , 200 1 ( L I NE 

(L)  , L  =  2 , 12) 

IFIBEGIN.NE.O.  IGO 

TO  41 

42 

CONTINUt 

DO  20  I  =  1,4 

DO  20  J  =  1,4 

LET  W( I, J)  =  0. 

20 

CONTINUE 

208 


ALL  ELEMENTS  OF  THE  STRESS-ENERGY  TENSUR,  MIXED  COMPONENTS 

W(l,J),  ARE  SET  EQUAL  TO  ZERO.  NOTE - THE  FIRST  INDEX  OF  W 

IS  DOWN  (COVARIANT)  AND  THE  SECOND  INOEX  IS  UP  (CONTRAVARIANT ) . 
USER  MUST  SUPPLY  "LET"  STATEMENTS  DEFINING  ALL  NON  ZERO 


ELEMENTS 

I 

LET 

W( 1,1  ) 

- 

LET 

W  (  2 , 2  ) 

X 

LET 

W  (  3 , 3  ) 

= 

LET 

W ( 4 ,4 ) 

= 

WR  I  TE ( 6 , 208 ) 

-P 

-P 

-P 

RHO 


User:  Replace  these  cards 
with  your  own  input. 


FORMAT ( 6 5H IS TRESS- ENERGY  TENSOR  NON  ZERO 
IS  (INPUT)  ) 


ELEMENTS,  MIXED  COMPONENT 


i  III. 0.2 


00  16  1*1.4 

00  16  J  *  1,4 

LET  Q  =  MATCH  ID,W( 

IFIQIGO  TO  16 
WR 1 TE  ( 6 ,505  )  I  •  J 

505  FORMAT! 1H  15,  1HL  14,  IHU  ) 

BEGIN  =  0. 

LET  ANS*W ( 1  ,  J  I 

96  LET  BEGIN  =  BCOCON  ANS  , LINE, 12 

WRITE (6,200) (L INE  I L  )  ,L*2, 12) 

IF ( BEG  IN. NE . 0. I  GO  TO  96 
ERASE  ANS 
16  CONTINUE 

WR 1 TE ( 6,87fl  ) 

878  FORMAT!  35H 1  THE  COVARIANT  CHRISTOFFEL  SYMBOLS  > 
00  3  I  =  1,4 

00  3  J  =  1,4 

DO  3  K  *  1,4 

LET  T 1  =  0. 

GO  TO  (31,32,70,71),! 

31  LET  T 1  *  FMCOI FIGIIJ.KI.X, 1) 

GO  TO  33 

32  LET  T 1  *  FMCDIFIGIIJ.KI ,Y, II 

GO  TO  33 

70  LET  T 1  =  FMCDIF(Gl(J,K),2,l) 

GO  TO  33 

71  LET  T 1  =  FMCOI F(G1(J,K),T,1) 

33  CONTINUE 

LET  T 2  =  0. 

GO  TOI35, 36, 72, 731, J 

35  LET  T 2  =  FMCOI F ( G 1  <  I  ,K ) , X , 1 ) 

GO  TO  34 

36  LET  T 2  =  FMCOI F C G 1 1  I , K I , Y , 1 1 

GO  TO  34 

72  LET  T 2  =  FMCOIF (G1(I,K),Z,1) 

GO  TO  34 

73  LET  T2  =  FMCOIF I G1  (I  ,K  I , T , 1 1 

34  CONTINUE 

LET  T3  =  0. 

GO  TO  (38, 39, 74, 75), K 

38  LET  T3  =  FMCOIF ( G1 II , J ), X , 1 ) 

GO  TO  37 

39  LET  T3  *  FMCOI F ( G  1C  I . J ) , Y , 1 ) 

GO  TO  37 

74  LET  T3  =  FMCOI F ( G1 ( I , J ) , l , 1 ) 

GO  TO  37 

75  LET  T3  *  FMCOIF IGI C I , J) ,T , l I 

37  LET  CCS ( I , J , K )  =  0.5*1 T 1*T 2-T3 ) 

LET  CCSU,  I,K)  =  CCS  ( I ,  J ,  K  1 
LET  0  =  MATCH  10,  CC S ( ! , J »K I , 0 . 

I F ( 0 ) GO  TO  3 
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WRITEI6.203) I ,J,K 

203  FORMAT ( 1H  I5.2HLS  I3.2HLS  I4.1HL) 

LET  ANS  =  EXPAND  CCS(I,J,K) 

LET  ANS  -  OROER  ANS,INC,FUL 
BEGIN  =  0. 

5  LET  BEGIN  *  6CDC0N  ANS, LINE, 12 
WRITE(6,200ICLINEILI,L=2,12) 

200  FORMAT! 1H  5X  12A6) 

IFIBEGIN.NE.O. )G0  TO  5 
ERASE  ANS 
3  CONTINUE 

ERASE  T1.T2.T3 
HR  I TE ( 6 »  209 ) 

209  FORMAT! 1H0) 

HR  I  TE ( 6 , 210 ) 

210  FORMAT ! 28H0ALL  OTHER  COMPONENTS  VANISH  I 

HR  I TE ( 6, 205  I 

205  FORMAT!  30H1 THE  MIXEO  CHRI STCFFEL  SYMBOLS  ) 
OQ  6  I  -  1,4 
DO  6  J  =  1,4 
00  6  K  *  1,4 
LET  MCS ( I , J , K )  *  0. 

00  7  L  *  1,4 

LET  TT  =  G2(K,L)*CCS(I,J,L) 

LET  TT  =  EXPAND  TT 
LET  SUM  =  TT 

IFIEPTERM.EO.OIGO  TO  502 
LET  SI  *  COEFF  TT,EP**0,R 
LET  S2  *  COEFF  TT,EP**1,R 
LET  S3  *  COEFF  TT,EP**2,R 
ERASE  TT 

LET  SUM  »  SI  ♦  S2*EP  ♦  S3*EP*EP 
ERASE  S1.S2.S3 

502  LET  MCSU.J.tO  =  MCSIl.J.K)  ♦  SUM 
ERASE  SUM 

7  CONTINUE 

LET  MCSIl.J.K)  *  EXPAND  MCSIl.J.K) 

LET  MCS(J.I.K)  *  MCS ( I , J.K ) 

LET  0  *  MATCH  10,  MCSI 1 , J , K > ,0. 

1  F  1 0 ) GO  TO  6 
WRITE(6,204) I , J.K 

204  FORMAT ! 1H  15.2HLS  I3.2HLS  I4.1HU  ) 

LET  ANS  *  OROER  MCS (  I , J.K ) , I NC, FUL 
BEGIN  -  0. 

8  LET  BEGIN  *  BCOCON  ANS, LINE, 12 

WRITE  I  6,200) (L INEI L  )  ,L«2, 12) 

IFIBEGIN.NE.O. )G0  TO  8 
ERASE  ANS 
6  CONTINUE 

00  350  I  *  1,4 
00  350  J  ■  1,4 
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i 


r>  r>  o 


{  iii.d.2 


00  350  K  =  1,4 
ERASE  CCS(  I  • J,K ) 
350  CONTINUE 

WR I TE ( 6, 209  I 
MRITE(6,210) 


THE  EQUATIONS  OF  NOTION  ENIMUI 
UR  I  TE 16, 207  ) 

207  F0RMAT(47HITHE  EQUATIONS  OF  NOTION,  COVARIANT  CCNPCNENTS  ) 
DO  21  NU  =  1,4 
LET  SI  *  0. 

DO  22  NU  «  1,4 
LET  T 1  *  0. 

GO  TO  (23,24,25,261  ,NU 

23  LET  T 1  *  FMCDIF  <  W(  NU,  NU I , X , 1 1 

GO  TO  27 

24  LET  T 1  *  FMC01F(W(NU,NUI,Y,1) 

GO  TO  27 

25  LET  T 1  *  FMCDIF(H(NU,NUI,2,1) 

GO  TO  27 

26  LET  T 1  =  FHCOIF(H(NU,NU),T,lI 

27  CONTINUt 

LET  SI  *  SI  ♦  T1 

LET  SI  *  EXPAND  SI 

22  CONTINUE 

LET  S2  «  0. 

LET  S3  *  0. 

00  28  I  -  1,4 
DO  28  NU  «  1.4 

LET  T 1  *  MCS( NU, I , NU I *U  I  NU, I  I 
LET  T 1  «  EXPAND  T1 
IF!EPTERM)508,509,508 

509  LET  S2  *  S2  ♦  T1 

ERASE  T1 
GO  TO  512 
508  CONTINUE 

LET  SO  *  COEFF  Ti,EP**0,R 
LET  S4  «  COEFF  T1,EP**1,R 
LET  S5  *  COEFF  T1,EP**2,R 
ERASE  T 1 

LET  S2  «  S2  ♦  (SO  ♦  S4*EP  ♦  S5*EP*EP) 

ERASE  SO , S4 , S5 
512  CONTINUE 

LET  S2  «  EXPAND  S2 

LET  T2  -  HCS(NU,NU,II4W(  I,NU) 

LET  T2  *  EXPANO  T2 
IFIEPTERM»511,510,511 

510  LET  S3  «  S3  ♦  T2 

ERASE  T2 

LET  S3  ■  EXPANO  S3 
GO  TO  28 

511  CONTINUE 
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LET  SO  =  COEFF  T2,EP**0,R 
LET  S4  =  COEFF  T2,EP**1,R 
LET  S5  =  COEFF  T2,LP**2,R 
ERASE  T 2 

LET  S3  =  S3  ♦  (SO  ♦  S4*EP  ♦  S5*EP*EP)  s 

ERASE  S0.S4.S5 
LET  S3  =  EXPAND  S3 
2b  CONTINUE 

LET  EM  =  S1+S2-S3 

ERASE  S1.S2.S3 
HR  I TE ( 6 , 206 )MU 
206  FORMAT ( 1  HO  15,  1HL  I 
563  FORMAT (151 

LET  ANS  =  EXPAND  EM 
ERASE  EM 

LET  ANS  =  ORDER  ANS.1NC.FUL 
BEGIN  =  0. 

29  LET  BEGIN  =  BCDCON  ANS, LINE, 12 
WRITE (6,200) (L  INE( I L  » ,  IL*2.  12  I 
562  FORMAT! 12A6) 

IFIBEGIN.NE.O.  ) GO  TO  29 
ERASE  ANS 
21  CONTINUE 
2500  CONTINUE 
STOP 
END 
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§  III.D.3 


3.  Output  for  MOTION 


The  printed  output  for  MOTION  includes:  the  metric  components,  and 
,  end  the  stress-energy  tensor,  —  all  of  which  are  input  by  the  user  — ; 
the  Christoff  el  symbols  (eqs.  (2)  and  (3));  and  the  divergence  of  the 
stress-energy  tensor  ("equations  of  motion") 


(ai)a 


r  w  ^ 

Ji 


0) 


The  format  is  similar  to  that  of  EINSTEIN;  and  error  messages  are  identical 
to  those  for  EINSTEIN. 

As  an  example  of  output  from  a  computation  by  MOTION,  we  give  below  the 
results  for  the  interior  of  a  static,  relativistic  stellar  model.  The  metric 
is  the  static,  spherically  symmetric  metric  of  equation  (l),  and  the  stress- 
energy  tensor  is 


Here  p  (pressure)  and  p  (density  of  mass-energy)  are  functions  of  r  ■  X 
only.  This  metric  and  stress-energy  tensor  are  precisely  the  ones  included 
in  the  listing  of  MOTION  in  the  last  section.  The  output  follows: 
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METRIC  TENSOR  COMPONENTS  (INPUT) 

—  Identical  to  those  given  for  this  metric  on  page  23. 


iTR  ESS-ENERGY  TENSOR  NON  ZERO  ELEMENTS,  MIXED  COMPONENTS  (INPUT) 
-Pi 

21  2U 

-Pi 

3L  3U 

-Pi 

41  All 
RHCi 


CHRISTOFFEL  STMBOIB  (OUTPUT) 

-  identical  to  those  given  for  this  metric  on  page  22. 


THE  tUUATIONS  OE  MOTION,  COVARIANT  COMPONENTS 
11 

-S.OE-i+PeFMCCIHANO,  (X.I  )  )-5.0E-I*RHO*FMCOIF 

I  x ,  1 ) )  i 

a 

O.Oi 

31 

O.Oi 

a 

O.Oi 


ANU, (X, 1  )  l-FMCC I F I P , 
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4.  Sample  Problems  for  MOTION 

Input  and  output  for  one  sample  problem,  the  Interior  of  a  static  relativ¬ 
istic  stellar  model,  were  given  In  the  last  two  sections.  As  a  second  example, 
consider  a  star  in  even-parity,  nonradlal  pulsation.  The  metric  for  such  a 
star  was  given  in  equation  (7),  and  the  stress-energy  tensor  is  as  follows 
(cf.  Thorne  and  Campolattaro  1967): 

W11  -  W22  -  W33  -  -  p  +  r"2  e"^2  W  P((dp/dr)  -  rp(Wn), 

\k  «=  p  -  r"2  e'^2  W  ?t(bp/br)  +  (p  +  p)(&i/n) 

-  (p  +  p)  e"v  Pf  ^  -  r-2  e^2  (dw/St ) J  , 

-  (p+p)  r‘2  e'V2  Pf(3w/dt), 


m  (p+p)  e'v  (dPj/doHdv/dt), 
\2  -  -  (P  +  P)  r'2  (dPf/de)(dv/dt). 


(10) 


Here  p  (density),  p  (pressure),  and  7  (adiabatic  index)  are  functions  of  r 
only;  An/n  is  defined  by 

Wn)  ■  £-  r"2  e'^2(dw/dr)  -  l(/+l)  r-2  V  +  K  +  H^J  P^;  (ll) 

l  la  the  integer  order  of  P^(cos8);  V  and  V  are  functions  of  r  and  t  associated 
with  the  displacement  of  the  fluid  from  equilibrium;  and  the  remaining  quanti¬ 
ties  were  discussed  on  pages  2b  and  26. 
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To  prepare  input  for  computing  the  divergence  of  the  8 tree 8 -energy 
tensor,  we  sake  the  following  changes  to  F0HMAC  notation  In  addition  to 
those  discussed  on  page  25  : 

i  *  AJ,  p  ♦  P 

V  ♦  V,  p  ♦  RH0 

V  -*•  WF,  7  *  GAM 

Also,  because  we  need  answers  only  to  first  order  in  the  perturbation  from 
equilibrium,  we  attach  the  expansion  parameter  EP**2  to  the  perturbation 
quantities  AK,  BO,  HI,  H2,  V,  and  WF,  wherever  they  appear.  The  resultant 
input  and  output  are  as  follows: 
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nmm 


ATOMIC  X,Y,Z,T,EP,AL,AJ,ANU,ALAM,AK,HO,Hi,H2,V,WF,GAM,P,RHO 
DEPEND  (ANU, ALAM,GAM,P,RHQ/XI » I AL/Y ) 

DEPEND  (HO.Hl,H2.V,WF,AK/X,T) 


EPTERM  «  3 


LET  GUI. II 
LET  G1U. 41 
LET  G1 1 4.  1 ) 
LET  GII2.2) 
LET  GU3.3I 
LET  GH4.4I 
LET  G2( 1. 1 1 
LET  G2(l,4) 
LET  G2I4.1) 
LET  G2I2.2I 
LET  G2I3.3) 
LET  G2 ( 4.4 ) 


-<-AL*EP**2*H2U.)*FMCEXP< ALAMI 
AL*EP**2*Hi 
GUI. 4) 

-( 1 .-AK*AL*EP**2  )*X**2 
- ( 1 .-AK*AL*EP**2  )*X**2*FMCSIN(YI**2 
( AL*EP**2*H0*1. )*FMCEXP(ANU) 

-I AL*EP**2*H2*i. )*FHCEXP(-ALAM) 
AL*EP**2*H1*FMCEXPI-ALAM-ANU) 

G2( 1*4) 

-(AK*AL*EP**2U. I*X** 1-2 1 
-(AK*AL*EP**2+l. )*IX*FMCSIN( Y)  I** 1-2  I 
I-AL*EP**2*HQU. )*FMCEXP(-ANU) 


1 

2 

3 

1 

2 


1 


l 

1 

1 

2 
3 


LET  Wll.ll 

LET  Mil, 4) 

LET  W  ( 2 . 2  I 
LET  W ( 2  *4  I 

LET  W  (  3 , 3  I 
LET  W(  4 » 1  I 

LET  W 1 4 . 2  I 

LET  W ( 4  *4  I 


-l-AJ*lAJU.I*V*X**('-2)*AK*H2/2.-X**(-2)* 
FMCEXPI-ALAM/2.  )*FMCDIF<WF,  (X,  I)  )  )*AL*EP**2'i'GAM* 
P+AL*EP**2  *WF*X**(-2 I *FMCEXP (-AL AM/2.  I* 
FMCDIFIP,  IX.lll-P 

AL*EP**2*Hl*(P*RHOI*FMCEXPI -ANU )-AL*EP**2* 
IP+RHO)  *X**(-2)*FMCEXPI ALAM/2.-ANLI* 

FMCOIFI WF. ( T. 1 1  I 
Mil. II 

EP**2*<P*RHOI*FMCEXPC-ANU)*FMCDIFIAL,l  Y,  1) I* 
FMCDIFIY.IT.il) 

MU. II 

AL*EP**2*(P+RHQ)*X**I-Z)*FMCEXPI-ALAM/2. I* 
FMCDIFIWF,  (T.l) I 

-EP**2*( P+RHO) *X**l-2 I *FHCD I F I AL . ( Y.  1)1* 
FMCDIFIV.IT, 111 

(-AJ*IAJ+l.)*V*X**(-2)*AK+H2/2.-X**l-2l* 
FMCEXPI-ALAM/2. )*FMCDIF(WF, (X.l) I )*AL*EP**2* 
IP+RHO)-AL*EP**2*WF*X**(-2)*FMCEXP(-ALAM/2. I* 
FMCOIFIRHO, (X, l)  U-RHO 
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COVARIANT  METRIC  TENSOR  NON  ZERO  ELEMENTS  (INPUT) 

1L  1L 

-(-AL*£P**2,0*H2*l«C)*FMCEXP(ALAM)l 
1L  AL 

AL  *EP**2 . 0*H1 i 
2L  2L 

-(-AK*AL*EP**2,0+1,CI*X**2.0S 
3L  3L 

-(-AK*AL*EP**2.0+l.C)*X**2.Q*FMCSIN(Y)**2.0S 
At  AL 

(  AL*EP**2.0*H0+1.0)*FMCEXP( ANU) S 

CON  TRA VAR  I  ANT  METRIC  TENSOR  NON  ZERO  ELEMENTS  (INPUT) 
1U  1U 

-( AL*EP**2.0*H2*I.0)*FMCEXP<-ALAMH 
1U  AU 

AL*EP**2. 0*H1*FMCEXP ( - ALAM- ANU I  t 
2  U  2U 

-(  AK*AL*EP*'*2.0+1.0)*X**<-2.0)  S 

3U  3U 

-(  AK*AL*EP**2.0+1.0)*X**(-2.0I*FMCSIN(Y)**(-2.0)S 
AU  AU 

(  -AL*EP**2.O»H0+l.O)*FMCEXP(-ANU) S 


STRESS-ENERGY  TENSOR  NON  ZERO  ELEMENTS,  MIXED  COMPONENTS  (INPUT) 

1L  1U 

-(-AJ*(AJ+1.0)*V*X**(-2.0)*AK«-H2*5.0E-1-X**(-2.0)*FMCEXP(ALAM*( 
-5.0E-1I )*FMCDIF(WF,(X,1)  )  )  *AL*EP**2.0*GAM*P*AL*EP**2.0*HF*X**( 
-2.0)*FMCEXP(ALAM*(-5.0E-l))*FMCCIFIP, !X,1))-PS 
1L  AU 

AL*EP**2.0*Hl*(P+RHC)*FMCEXP(-ANU)-AL*EP**2»0*( P+RHO) AX**  (  -  2,0 ) * 
FMCEXPI ALAM*5.0E-i-ANU) *FMCDIF(WF,  (T,l  ) >S 
2L  2U 

- ( -A J  * ( AJ*1.0)*V*X**l-2.0)+AK+H2*5.0E-l-X**(-2.0)*FMCEXP(ALAH*( 
-5.0E-1) )*FMCDIF(WF ,(X,1) ) ) *AL*EP**2. 0*GAM*F*AL*EP**2.0*MF»X**( 
-2.0)*FMCEXP( ALAM*(-5.0E-1) >*FMCDIF(P,  (X,l) )-PS  1 

2L  AU 

EP**2.0*(P  +  RH0)*FMCEXP(-ANU)*FMCDIF(AL,(Y,1 ) ) *FMCC IF ( V,  ( T ,  1 ) ) S 
3L  3U 

-  ( -AJ*( AJ<-1.0)*V*X**(-2.0)+AK  +  H2*5.0E-1-X**(-2.0)*FMCEXP(  ALAM*( 
-5.0E-1) l*FMCDIF(WF,(X,ll  I  ) *AL*EP**2.  OeGAM*PAALAEP**2.0*WF*X** ( 
-2.0)*FMCEXP(ALAM*(-5.0E-1 ) )*FMCDIF  (  P,  < X,  1 1  )-P* 

AL  1U 

AL*EP**2.0*(P+RHO)*X**(-2.O)*FMCEXP(ALAM*(-5«OE-l))*FMCCIF(WF»(T,l 
)  IS 

AL  2U  ! 

-£P**2.0* ( P+RHO )*X**(-2.0) *FMCDI F( AL*  (Y,l) )*FMCDIF(V, (T,l)  )S 
AL  AU 

(  -  A J  * ( A J* 1 • 0 )  *V*X** ( -2 • 0 ) ♦ AKaH2*5,0E-1— X** ( -2,0 ) *FMCEXP ( ALAMA( 
-5.0E-I) )*FMCDIF(WF,(X,1) ) ) *AL*EP**2. 0* I P  +  RHO I -AL*EP**2.0*WF*X**(  j 

-2.0)*FMCEXP(ALAM*(-5.0E-1))*FMC0IF(RH0,(X,1 ) l+RHCS  j 

i 

CHRISTOFFEL  SYMBOLS  f 

—  Identical  to  those  on  pages  27-29,  1 

54  I 


§  in.D.4 


THE  EQUATIONS  OF  MOTION,  COVARIANT  COMPONENTS 

1L 

-2.0*AJ*At*EP**2.0*GAM*P*V*X**(-3.0)+5.0E-l*AJ*AL*EP**2.0*GAM*P*V* 
X**(-2.0I *FMCOIF( ANU, ( X,1 ) ) *A J* AL*EP**2 . 0*GAM*P*X** ( -2 . 0  I *FMCD I F  ( V 
, (X, L ) 1*AJ*AL*EP**2.3*GAM*V*X**(-2.0)*FMC01F(P» (X, 1 1  I +5.0E- l*AJ*AL 
*EP**2.0*P*V*X**(-2.0I *FMCD IF ( ANU, (X, 1) 1 ♦ AJ*AL *EP* *2 . 0*P* V*X**  I 
-2.0)*FMCDIFIGAM, (X, 1) I +5 .  Ot- 1* A J* AL*EP**2 . 0*RHO*V*X** I -2 . 0 1 * 
FMCDIF I ANU, I  X, 1) ) -2 . 0* A J**2 .0* AL *EP**2 . 0*GAM*P* V*X** ( -3 . 0 ) ♦ 
5.0E-1*AJ**2.0*AL*EP**2.0*GAM*P* V*X**{ -2.01 *FMCDIF( ANU, (X, l ) )*AJ** 
2.0*AL*EP**2.0*GAM*P*X**(-2.0I*FMCDIFI V, I  X , 1 > 1 + A J**2 . 0* Al*E P**2 .0* 
GAM*V*X**(-2.0)*FMCCIFIP,(X,11 1 *5 . OE-1 *A J**2 . 0* AL*EP**2 .0*P*V*X**  ( 
-2.0)*FMCDIF(ANU,(X,1) 1 +A J**2 • 0* AL*E P**2. 0*P*V*X**  I -2.0I*FMCDIF( 
GAM, (X, 1) )+5.0E-l*AJ**2.0*AL*EP**2.0*RHO*V*X**(-2.0)*FMCCIF(ANU,IX 
, 11 1-5. OE-1 *AK*AL*EP**2.0*GAM*P* FMCDIF (ANU, (X, 1 1 1 - AK*AL *EP**2 . 0* 
GAM*FMCOIF( P, ( X, I  1 1-5. OE-1* AK*AL*EP**2.0*P*FMCDIFIANU, I  X, 1 1 1 -AK*AL 
*EP**2.0*P*FMC0IF(GAM, (X, 1) 1 -  5 . OE- 1*AK*AL*EP**2 . 0*RHO*FMCD I F ( ANU, ( 
X, 1)1-2. 5E-1*AL*EP**2.0*G AM *H2*P*FMCDIFIANU,(X,1 1 1-5 .OE- 1*AL*EP»* 
2.0*GAM*H2*FMC0IF(P, IX, II )-2.0*AL*EP**2.0*GAM*P*X**(-3.0)*FMCEXP(- 
5.0E-1*ALAM)*FMC01F(WF,(X, 1 ) ) -5 . OE- 1*AL*EP**2 . 0*GAM*P*X** l -2 .0 1  * 
FMCEXP(-5.0E-1*ALAM)*FMCDIFIALAM, IX, l) 1 *FMCG I F  1  WF , I  X , 1 1 )♦ 
5.0E-1*AL*EP**2.0*GAM*P*X** 1-2.0 )*FMCEXP( -5.0E- 1* AL AM  1 *FMCD I F I ANU, 
( X.ll )*FMC0IFIWF,(X,1))+AL*EP**2.0*GAM*P*X**I-2.0)*FMCEXP(- 
5.0E-1*ALAMI*FMCDIF(WF,(X,2 )>-AL*EP**2.0*GAM*P*FMCClF(AK,(X,l) 1- 
5.0E-1*AL*EP**2.0*GAM*P*FMCD1F(H2, IX, 1 1  1 ♦ AL*EP**2. 0*GAM*X** ( -2.0  I ♦ 
FMCEXP«-5.0E-l*ALAM)*FMC0IFIP,(X,l)  I *FMCD I F I WF , I  X , 1 1 1 -2 . 5E- 1* Al*EP 
**2.0*H2*P*FMCDIF( ANU, IX, 1 1 1-5 . 0E-1*AL*EP**2. 0*H2*P*FMCC I F  I  GAM, (X, 
1 1 1-2. 5E-1* AL*EP**2 .0*H2*RH0*FMCDIF I ANU, I  X, 1 1 1  *5.  OE- 1  * AL*EP**2. 0*P 
*X**(-2.0)*FMCEXPI-5.0E-1*ALAM)*FMCDIFI ANU,(X,1) )*FMCDIF(WF, IX, 1) I 
♦AL*EP**2.0*P*X**(-2.0I*FMCEXP( -5. OE- 1*AL AM  1 *FMCC IF ( GAM , I  X , 1)1* 
FMCDIF I WF , ( X , 1 1 |-AL*EP**2.0*P*X**(-2.0)*FMCEXP(5.0E-1*ALAM-ANU>* 
FMCDIFIWF, IT ,2) )+AL*EP**2.0*P*FMCEXP(-ANUl*FMCDIFIhl, I T , 1 1 1 - 
5.0E-1*AL*EP**2.0*P*FMCDI  F I  HO, I  X , 1 1  1  +  5 .0E-1*AL *EP**2 . 0*RH0*X** 1 
-2.01*FMCEXP(-5.0t-l*ALAM|*FMCOIFI ANU, (X, II )*FMCDIF! WF, IX, 1 1 l-AL* 
EP**2.0*RHQ*X**(-2.0)*FMCEXP(5. OE-l *AL AM-ANU )*FMCDIFIWF,(T,21 l  +  Al* 
EP**2.0*RH0*FMCEXPI-ANU)*FMCDIFI Hi, IT, 1 1  )-5.0E-l*AL*EP**2.0*RH0* 
FMCDIF (HO, (X, 11 >-2.0*AL*EP**2.0*WF*X**(-3.0)*FMCEXP(-5.0E-l*ALAM)* 
FMCDIF (P, ( X, 1 1 )-5.0e-l*AL*EP**2.0*WF*X**(-2.0)*FMCEXP(-5.0E-l*ALAM 
)*FMC0IF(ALAM,(X,1)  1  *FMCD I F I P, ( X , 1 1  1  +5 . 0E-1*AL*EP**2 . 0*WF*X** I -2. 0 
)*FMCEXP(-5.0E-1*ALAM)*FMCDIF(ANU, (X, 1 1 )*FMC0I FI P « I  X , 1 1 1  + 
5.0E-1*AL*EP**2.0*WF*X**(-2.0)*FMCEXP( -5 . OE-1* AL AM  1 *FMCC IF ( ANU, (X, 
1 1 )*FMCDIFIRHO, I  X, 1 1 1 *AL*E P**2 . 0*WF*X** I -2. 0 1 *FMCEXP ( -5.0E- 1* AL AM  1 
♦  FMCDIF IP, IX, 2  I |f-AL*EP**2.0*X**(-2.0)*FMCEXP(-5.0E-l*ALAM)*FMCDIF( 
P,IX,1) )*FMCDIF(WF, (X,1)I-5.0E-1*P*FMCDIF(ANU, I  X, 1 1 )-5.0E-l*RH0* 
FMCDIF (ANU, <X,1) 1-FMCDIFIP, (X,l)  )» 
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THE  EQUATIONS  OF  MOTION  (continued) 

2L 

A J*EP**2.0*GAM*P*V*X**<-2 .0  l*FMCUIFt  AL , ( Y, 1 ) ) +AJ**2.0*EP**2 . 0<-GAM* 
P*V*X**<-2.0)*FMCUIF(AL,<Y,1) ) - AK*E P**2 . 0*GAM*P *FMCD [ F  I  AL,  ( Y, 1 ) )- 
5.0E-1*EP**2.0*GAM*H2*P*FMCDIF( AL, ( Y, 1 ) ) +EP**2 • 0*GAM*P*X**  < -2.01* 
FMCfcXP(-5.0E-l*ALAM|*FMCDIF( AL, t  Y  .  I ) ) *FMC DI F ( WF , I  X, 1  I ) - 5. OE- 1 *E P** 
2 . 0+HO*P*FMCG I  FI AL , ( Y» 1 ) I -5 .0E-i*fcP**2 . 0*H0*RHO*FPC0IF ( AL  ,  ( Y ,  1 )  )♦ 
EP**2.0*P*FMCEXP(-AMI)*FMCCIFCAL,  < Y , 111 *FMCD l F ( V , ( T,  2 » ) +EP**2 . 0* 
RHO*FMCEXP(-ANU)*FMCDIF(ALi (Y,1)|*FMCDIF(V,(T,2)  ) +EP**2 . 0*WF*X** ( 
-2.0l*FMCEXP(-5.0fc-l*ALAM)*FMCOtF(  AL  ,  < Y , 1 ) ) *FMCD I F  (P,  (  X , 1  ) )» 

3L 

0  •  0 1 

AL 

-AJ*AL*SP  +  *2.0*P*X**  <-2.0)*FMCDIF(  V, IT, 1  ) )-AJ*AL*EP**2.0*RH0*X**| 
-2.C)*FMCDIFIV»ir,l) I-AJ**2.0*AL*EP**2.0*P*X**(-2.0)*FPCCIF(V»( T, l 
) )-AJ**2.0*AL*l-P**2.O*RHO*XP*(-2.O»*FMCOIFI V, ( T,  1  l)+5.0E-l*AL*EP** 
2.0*P*X**(-2.O)*FHCeXPI-5.Ot-l*ALAM)*FMCDIFC ANU. ( X, 1 1 ) ♦FMCD I F ( WF , ( 
T  » 1  I  1+5.0 c-l+AL*EP**2.Q*RHl)*X**(-2.0)*FMC£XP(-5.0E-l*ALAP)*FMCUlFI 
ANU.  (  X,  1  (  )  ♦FPCD I  F  ( WF ,  <r,in+AL*EP**2.0*X**(-2.0)*FPCEXP(- 
5.0E-1*ALAP)+FPC0IF(P,(X,1)  ) *FMCO  I  F { WF , (T,l))-EP**2.0*P»X**(-2.0)* 
FMCS!M( Y)**(-l .U)*FNCCOb< Y I *FMCL I F  (  AL  ,  I Y , 1 ) I *FMCD  l  F  <  V » ( T,1 ) )-EP** 
2.0*P*X**t-2.0)+FPCLlF(AL ,  I Y , 2 > ) ♦FMCU I F I  V , I T , 1 ) 1 -EP**2 . 0*RHQ*X** ( 
-2.0)*FMCSIN(Y I  *  * (  —  L . 0 ) ♦FMCCOS ( Y ) *FMCO I  FI AL , ( Y ,  1  )  )*FMCDIF<  V. ( T,  1 )  ) 
-EP**2.0+KHU*X*+(-2.0 ) *FMC(ilF(ALi(Y,2)|*FMCCIF(Vt(Til))t 
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E.  SUBSTITlfrE,  tne  Prograa  for  Maying  Changes  of  Variables 

1.  Input  for  SUBSTITUTE 

SUBSTITUTE  is  the  program  which  takes  long,  complicated  expressions  and 
makes  changes  of  variables  In  them.  The  Input  for  SUBSTITUTE  Is  handled  in 
a  slightly  different  manner  than  that  for  EINSTEIN,  RUMANN,  and  MOTION: 

Four  statements  are  supplied  by  the  user  In  the  body  of  the  program,  and  the 
remaining  input  information  is  supplied  on  "data  cards"  at  the  end  of  the 
program* 

The  statements  supplied  in  the  body  of  the  program  are  these: 

1.  ATOMIC  and  DEPEND  statements.  In  these  statements  must  be  listed 
all  atomic  variables  for  the  original,  unsubstltuted  expressions 
and  for  the  substitutions  which  are  to  be  made.  For  example,  if 
A(x)  in  the  expression 

is  to  be  replaced  by 

A  =  -  In  (1-2VX),  (13) 

then  A,  X,  m  DU6t  all  be  declared  in  the  ATOMIC  statement;  and 
A  and  m  must  both  be  declared  in  DEPEND. 

2.  ERTERM  statement.  EETERM  mist  be  defined,  as  in  the  program 

EINSTEIN,  to  tell  the  computer  whether  or  not  to  truncate  the  ; 

answers  at  second  order  in  EP.  (cf.  page  9.) 

3.  A  "LBL  PARAM"  statement.  This  statement  contains  a  list  of  the 
names  of  all  atomic  variables  for  which  substitutions  are  to  be 
made.  For  example,  in  the  situation  described  above,  A  appears 

I 
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In  the  LBL  PARAM  statement.  The  entiles  In  the  LBL  PARAM  list  must 
appear  In  precisely  the  same  order  as  Is  used  for  the  input  data 
(cf.4b  below).  The  form  of  the  LBL  PARAM  statement  Is  as  follows: 

col. 3  col. 7 

+  + 

LBL  PARAM  (NAME1,  AA(l)),  (NAME2,  AA(2)),  - 

(NAMEN,  aa(n)). 

Here,  MAME1,  — ,  NAMEN  are  the  names  of  the  atomic  variables  for 
which  substitutions  are  being  made.  No  more  than  50  changes  of 
variables  can  be  made  unless  the  user  redimensions  AA. 

4.  The  data  supplied  by  the  user  at  the  end  of  SUBSTITUTE  is  punched 
on  cards  with  the  first  entry  of  each  card  in  column  1,  and  the  last 
entry  before  or  in  column  72.  The  data  which  must  be  supplied  are 
as  follows: 

a)  Two  integers,  on  a  single  card,  giving  the  number  of 
substitutions,  "NSUB",  to  be  made,  and  the  number  of 
complicated  expressions,  "NEXP",  to  be  expanded  with 
the  substitutions.  The  first  of  these  integers,  "NSUB", 
must  end  in  column  5  and  the  second,  "NEXP",  must  end  in 
column  10.  In  this  card,  and  this  card  only,  the  data 
need  not  start  in  column  1. 

b)  Next  come  the  NSUB  substitutions  —  i.e.,  expressions  for 
the  old  atomic  variables  in  terms  of  the  new  ones.  These 
are  read  in  one  at  a  time,  each  preceded  by  a  card  listing 
a  label  or  title  of  the  quantity  being  read  in,  (these 
labels  can  be  72  characters  in  length),  and  each  followed  by  a 
blank  card.  The  substitution  expressions  must  be  read  in 
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with  the  sane  order  as  Is  used  to  list  them  In  the 
LBL  PARAM  statement  (see  above).  Each  substitution 
expression  can  be  many  cards  long  -  with  each  card's 
entries  beginning  in  column  1  and  ending  before  or  in 
column  72.  The  computer  identifies  the  end  of  a  sub¬ 
stitution  expression  by  a  dollar  sign,  $,  with  which 
the  user  must  terminate  each  substitution, 
c)  The  last  quantities  to  be  read  in  are  the  NEXP  complicated 
expressions  which  are  to  be  expanded  in  terms  of  the  new 
atomic  variables .  These  are  read  in  one  at  a  time,  each 
preceded  by  a  label  card,  each  terminating  in  a  dollar 
sign,  and  each  followed  by  a  blank  card. 

A  listing  of  SUBSTITUTE,  including  sanple  data  which  illustrate  the 


above  instructions,  follows. 


n  n  n  n  nnnnnnnn  n  n  n  n  n  n 


{  m.E.2 


2.  A  Complete  Listing  of  SUMITIWR 


SIBFMC  ONE  NODECK 
SYMARG 

PLACE  ATOMIC  AND  DEPEND  STATEMENTS  FOLLOWING  THIS  COMMENT 

atomic  X  *  Y  *  Z  *  T  *alam,and*FM*PHI  \  Users  Replsce  these  cards 
DEPEND  (ANUtALAM«PHI  .EM/X)  J  with  your  own  Input. 


DIMENSION  CARD! 600 ) ,AA( 50) 
DATA  BLANK  /6H  / 


INTEGER  EPTERM 


IF  EPTFRM  =  0  TERMS  CONTAINING  ALL  POWERS  OF  EXPANSION 
PARAMETER  WILL  BE  RETAINED. 

IF  EPTFRM  .NF.  0  TERMS  CONTAINING  THE  EXPANSION  PARAMETER 

TO  A  POWER  GRFATER  THAN  2  WILL  BE  DISCARDED. 

A  CARD  CONTAINING  THE  INTEGFR  VALUE  OF  EPTERM  MUST  FOLLOW  THIS 
COMMFNT. 


EPTFRM  =  0 


) 


User:  Replace  this  card 
with  your  own  input. 


WR  I  TE ( 6 .202 1 

202  FORMAT (57H1SUBST I  TUT  IONS  TO  RE  MADE  IN  CURVATURE  TENSOR  EXPRES5I0 
INS  ) 

WRITEI6.2O0) 

200  FORMAT (1H0) 

RE AD ( 5 , 1 01 1NSURS.NEXP 

101  F  ORMAT (215) 

DO  10  II  =  1.NSIJPS 

READ!  5  *102  MCARDI  J)  »J=1 .1?) 

102  F  ORMAT ( 1 2  A  6 ) 

WRITE (6*105) (CARD! J) *J=1 ,12) 

105  F0RMATI1H  5X  1 2 A 6 ) 

I  =  0 

1  CONTINUE 
JJ  =  1+1 
JJJ  =  1+12 

READ! 5,102  M  CARD(J ) , J  =  JJ.  JJJ) 

IF  (CARD!  JJ)  .EQ.RLANOGO  TO  B 
WRITE (6*106) (CARD( J) ,J=JJ,JJJ) 

106  FORMAT110X.12A6) 

1=1+1? 

GO  TO  1 
3  K  =  0 

LET  AA(II)  =  ALGCON  CARD,K 
10  CONTINUE 


PLACE  THE  LBL  PARAM  STATEMENT  FOLLOWING  THIS  COMMENT 


C 

LBL  PARAM  (  ALAM,AA< 1  )  )  *  ( ANU , A A ( 2 ) ) 
C 


,  User:  Replace  this  card 
'  with  ycur  own  input. 
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DO  20  II  =  ltNEXP 
WR I TE ( 6  » 1 04 ) 

104  FORMAT  < 1 H 1  ) 

WRITF<6,200) 

RFAD<  5.222 ) ( CARD ( J ) ,J=1 .1?) 

222  FORMAT ( 1 2 A  6 ) 

WRITF (6*223) (CAPD(J) « J= 1  *  1 2 ) 

221  FORMAT ( 5X  *  1 2 A6 ) 

I  =  0 

1?  CONTINUE 
JJ  =  I  +  1 
JJJ  =1+12 

RF AD ( 5  *  I  02 ) (CARD  (  J ) »J=JJ*JJJ) 
IF(CARO(JJ).EQ.BLANK)GO  TO  2 
WR I T  E ( 6  « 1 05 ) (CARD!  J)  *J  =  JJ*JJJ> 

I  =  I  +  12 
GO  TO  12 
2  K  =  0 

LET  EXP  =  ALuCON  CARD « K 
LET  ANS  =  SUBST  EXP  »LRL 
IF(EPTFRM.EQ.O)GO  TO  32 
LET  ANS  =  EXPAND  ANS 
ERASE  EXP 

LET  ZO  =  COEFF  ANS,EP**O.PP 
LET  Z1  =  COEFF  ANS,FP*»1*PP 
LET  Z2  =  COEFF  ANS,EP**2.PP 
ERASE  AMS 

LFT  EXP  =  FP**0*Z0  +  EP**1*Z1  +  FP**2*72 
GO  TO  33 

32  CONTINUE 

LET  EXP  =  ANS 

33  CONTINUE 

LET  EXP  =  EXPAND  FXP 
WR  I  TE ( 6  » 1 07 ) 

107  FORMAT (34H0THE  EXPRESSION  AFTER  SUBSTITUTION  ) 

WR I TE (6*200) 

S  =  0. 

11  LET  S  =  BCDCON  EXP  ,  CARD  ♦  12 
WRITE (6. 105) <CARD(  J) *J  =  2»12  > 

IF (S.NF.O. )G0  TO  11 
20  CONTINUE 
STOP 
END 
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Drt»  far  3UB6TITVIE 


(User:  Replace  these  data  cards  by  your  own  data  cards.) 


2  4 

ALAM 

—  FMCLOG ( 1 -2#£M/X )  $ 

ANU 

2*PH  I  * 

1U  1L  EINSTEIN  TENSOR  COMPONENT 

-X**(-2)»FMCEXP(-ALAM)  +  X**<-2)  -  X** ( -1 ) *FMCEXP I - ALAM ) *FMCD I F ( ANU.X ♦ 1 > 
* 


2U  2L  EINSTEIN  TENSOR  COMPONENT 

1,/<X«2. )*FMCFXP(-ALAM)  *  FMCHIF ( ALAM.X .1 )  -  l./(X*2.) 
*FMCEXPI-ALAM)*FMCD1F< ANU.X.l )  +  FMCFXP I-ALAM) *FMCDI  F ( ALAM ,X . 1 >  /4. 

•FMCDIFIANU.X.l >  -  FMCEXP ( -AL AM  )  *FMCD I F ( ANU . X . 1 ) **2 / 4.  -  FMCEXP ( -ALAM i 
*FMCD I F ( ANU »X  *21/2.  $ 

31J  3L  EINSTEIN  TENSOR  COMPONENT 

l./I X«2. > *FMCEXP I-ALAM) *FMCDIF ( ALAM.X. 1 )  -  1. / (X*2. ) *FMCEXP I -AL AM ) 
•FMCDIF I ANU.X.l )  +  FMCEXP I -AL AM > «FMCD I F I ALAM , X , 1 ) *FMCD I F I ANU . X . 1 > /4 . 
-FMCEXP (-ALAM)*FMCD IF ( ANU.X.l ) **2/4.  -  FMCEXP ( -ALAM ) *FMCDI F ( ANU .X ,2 ) 

/?»  t 

4U  4L  EINSTEIN  TENSOR  COMPONENT 

-1«/(X*X)*F MCE XP(— ALAM)  +  l./(X*X!  +  1 . /X*FMCEXP I -AL AM ) *FMCD I F I  AL AM ,X . 1  ) 
S 
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3.  Output  of  aUBBTITOTE 

The  output  for  SUBSTTTifTE  begins  with  a  reproduction  of  all  the  substitu¬ 
tions  (changes  of  variable)  to  be  made.  Then  follow,  in  the  order  in  which 
they  were  input,  the  various  coaiplicated  expressions  to  be  expanded.  These 
are  printed,  first  In  their  original  form,  and  lmedlately  thereafter  expanded 
in  terns  of  the  new  atonic  variables. 

As  an  example,  here  is  the  output  generated  by  the  data  which  were  given 
in  the  last  section  (N8UB  «  2,  NEXT  -  4)l 

SUBSTITUTIONS  TO  BE  MADE  IN  CURVATURE  TENSOR  EXPRESSIONS 
ALAM 

-FMCLOGI 1-2*EM/X)  * 

ANU 

2*PHI  1 


1U  1L  blNSTblN  TENSOR  COMPONENT 

-X**I-2)*FMCEXPI-ALAM)  ♦  X**(-?)  -  X** I - 1 ) *FMCEXP I -AL AM )*FMCDIFI ANU» X, 1 ) 
t 

THE  EXPRtSSION  AFTER  SUBSTITUTION 


(-3.0) *2.0*EM*X**(-2.0)*FMCDIF( PHI , IX, 1 ) )*4.0*X**(-l .0)* 
FMCOI F (PHI,(X,1) I  *  I -2 .0 ) t 
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2U  21  EINSTEIN  TENSOR  COMPONENT 
l./IX*2.)*FMC£XPI-ALAM)  *  FMCD I F I AL AM , X , 1 )  -  l./(X*2.) 

♦  FMCEXP ( -ALAM )*FMCDI FI ANU,  X,l)  ♦  FMCEXP I -A L AM ) *FMCDI F I AL AM, X , I >/4. 

+  FMCD I F ( AND,  X , 1 )  -  FMCEX P ( -AL AM ) *FMCDl F I ANU, X , 1)  **2/4 .  -  FMCEXPI -ALAM) 
*FMCDlF(ANU,X,2)/2.  * 

THE  EXPRESSION  AFTER  SUBSTITUTION 


-EM*X**(-3.0)+EM*X**I-2.0)*FMC0IFIPHl,(X,l))*EM*X**<-1.0)*FMCDIF( 
PHl,IX,l))**2.0*2.0+EM*X**l-1.0>*FMCDIF(PHI, CX,2!)*2.lHX**(-2.0>* 
FMCDIF(EM,IX,m*X**l-l.O)*FMCDIF(EMt  IX,1)  >*FMCDI F I  PH  I ,  I  X ,  1 >  > -X** I 
-l.O)*FMCDIF(PHI,IX, IH-FMC0IFIPHI,(X,1I)**2.0-FMCDIFIPHI,  tX,2l) * 


3U  3L  EINSTEIN  TENSOR  COMPONENT 
l./(X*2. )*FMCEXPI-ALAMl*FMCOIF(ALAM,X,l)  -  1 ./(X*2. ) *FMCE XP I -ALAM ) 
♦FMCDIFIANU.X.l)  ♦  FMCEXP ( -ALAM > *FMCD I F ( AL AM , X , 1) *FMCDI F  (  ANU, X , 1 ) /4. 
-FMCEXP(-ALAM1*FMCDI FI ANU, X, 1 )**2/4.  -  FMCEXP(-ALAM)*FMCD!FIANU,X,2> 

/2.  * 

THE  EXPRESSION  AFTER  SUBSTITUTION 


-EM*X**I-3.0)*EM*X**C-2.0)*FMC0IFIPHl, I x, 1))+EM*X**(-1.0)*FMCDIF I 
PHI , <X,1)I**2.0*2.0»EM*X**(-1.0)*FMCDIF(PHI, I  X , 2  > ) *  2 . 0*X* *  I -2 . 0 > * 
FMCOIFIEM, (X,ll >*X**I-1.0)*FMCDIFIEM, IX, l 1 >*FMCOIF(PHI, (X, 1 ) ) -X** I 
- 1 . 0 ) ♦ FMCD I F I  PHI  ,  I  X, 11 l-FMCOIFI PHI , I  X , 1  I ) **2 . 0-FMC DI F 1  PHI , ( X, 2 ) ) ( 


4U  4L  FINSTEIN  TENSOR  COMPONENT 

-l./(X*XJ*FMCEXP!-ALAM»  ♦  l./IX*X>  +  1 . / X*FMCEXP( - AL AM) *F MCU I T I AL AM, X , 11 

* 

THE  EXPRt-SSION  AFTER  SUBSTITUTION 


X**(-2.0)*FMCI)IFIrM, IX, l >1*2.0* 
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4.  Sample  Problems  for  SUBSTITUTE 

Input  and  output  for  one  sample  problem  were  given  in  the  last  two  sec¬ 
tions.  As  a  second  example  we  give  input  and  output  for  a  problem  in  which 
limited  core  storage  actually  forced  us  to  use  SUBSTITUTE: 

Chitre,  Hartle,  and  Thorne  (1967)  have  been  studying  recently  the  funda¬ 
mental  modes  of  pulsation  far  slowly  rotating,  relativistic  stellar  models. 

p 

To  second  order  (e  )  in  the  rotation  and  first  order  (5 )  in  the  pulsation, 
the  metric  for  such  a  star  can  be  put  into  the  form 

gii  ■  e*  t1  +  2  *2  r_1  *^*0  +  ¥a*  +  2  5r"1  *\  ■  +  ¥2)} 

gi4  “  g4i  “  * €  ^  ¥2 

822  -  r2  {1  +  2  €2  (v2  -  h2)P2  -  €2t  KgPg} 

833  “  1-2  Bin2°  l1  +  2  «2  (v2  -  b2)P2  -  «2£  K2P2) 
g34  "  g43  *  p2  sln2®  -  °)  -  eC.lx) 

«44  “  ‘  *V  (I  +  2  2  (hQ  +  hgPg)  -  «2  e‘vr2  ain28  (3  -  n)2 

+  (»0  +  ¥2)}  (14) 

1  2  3  U 

where  x  *  r,  x  «  Q,  x  ■  <p,  x  -  t,  and 

Pg  -  i  (3  cos2a  -  1).  (15) 

In  this  metric  X,  v,  w,  hp,  hg,  nig,  nig,  and  Vg  are  functions  of  r  only; 

Jx,  Hq,  Hg,  Nq,  Ng,  and  Qg  are  functions  of  r  and  t;  and  e,  5,  0  are  constants. 

This  metric  is  so  complicated  that  EINSTEIN  was  unable  to  calculate  the 
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corresponding  alxed  Christoff  el  symbols  or  Ricci  tensor  or  Einstein  tensor. 
An  error  number  7337  (co^mter's  memory  exhausted)  vas  obtained  when  the 
computation  was  attempted. 

To  get  around  this  difficulty,  a  simplified  metric,  equivalent  to 
expression  (it)  hut  containing  fewer  terms,  was  put  Into  EINSTEIN: 


gll 

- 

^  {1  +  «2  B]  , 

glb 

- 

«bl  *  -  e2  5q  » 

®22 

- 

r2  [1  +  e2  C}  , 

*33 

- 

r2  sln20  {1  +  c2  C]  , 

Hk 

a 

843  -  -  «(r2  sln2e)  V, 

- 

-  ew  (1  +  e2  <33  +  e2  (r2 

u 

B 

m 

(1  -  c2  flj  , 

It 

8 

- 

g*1  -  .  c2  e”v  Pe 

22 

8 

- 

r-2  (1  -  €2  a)  , 

33 

8 

- 

(r  BlnQ)"2  {1  -  e2  CJ  -  e: 

3b 

8 

= 

g^3  »  -  c  e’v  ?lf. 

bb 

8 

c 

-  e*v  (l  -  e*aj , 

(16a) 


(16b) 


Here  the  new  variables  9^,  9e,  and  tf  are  functions  of  r  and  t;  while  a,  B, 
C,  and  5^  are  functions  of  r,  0,  and  t.  EINSTEIN  was  able  to  compute  the 
Bled  and  Einstein  tensors  for  this  s  1^)11  fled  metric  to  order  €2  vithout 
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| 

i 

exhausting  the  computer's  memory.  The  resultant  Ricci  and  Einstein  tensors 
vere  then  put  through  SUBSTITUTE,  which  made  the  following  changes  of  variable 
-  accurate  to  first  order  in  {  -  to  return  them  to  the  original  notations 

\  -  e*  (1  +  2  Cr*1  e*  aQ)  , 

■  e"X  I*  “  2  tr'1  e*  mq)  , 

sjf  .  0  -  <3  +  , 

a  -  2  ^“^(<5  -  u)r2  ein28  +  2  hQ  +  £NQ 

(IV) 

+  (2  ^  +  5N2)  }  (3  cos2®  -  1)  , 

B  -  2  e*  {n^j/r  +  (fflg/r)  $  (3  cos28  -  l))  {1-2  Ce^Vg/r] 

-  5  {Hq  +  h2  i  (3  cos2e  -  1)) 

a  -  {2  (v2  -  hg)  -  SK2)  i  (3  cos2e  -  1)  , 

yq  -  t  «2  i(3  cos2e  -  i). 

In  making  these  changes  of  variable  by  means  of  SUBSTITUTE,  ve  exhausted  the 
computer's  memory  on  a  few  of  the  components  of  the  Einstein  tensor.  When 
that  happened,  we  broke  the  difficult  components  into  pieces  and  performed 
the  substitutions  piece  by  piece. 

As  a  guide  to  the  user  of  SUBSTITUTE  who  has  a  long  problem  to  do,  ve 
reproduce  here  the  input  to  and  the  output  from  SUBSTITUTE  for  part  of  the 
above  problem.  Mote  that  SUBSTITUTE  does  not  truncate  in  powers  of  { 

(denoted  EZ  below),  ao  that  truncation  must  be  done  by  hand. 
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not 


ATOMIC  X  .  Y  *Z  »  T  .EP.EZ 

ATOMIC  AL AM* ANU«PI X .PIE.A.B.C.FQ.WR 

ATOMIC  AMUO.AJ1.HO.H2.CHO.CH7.FNO.FN2.EMO,EM2.AK2.V2.Q2.0M,OMR 
DEPEND  (ALAM.ANU,OMn.HO*H2.FMO*EM2*V2/X) 

DEPEND  IP1K*PIE*WR  * AMUO  * A J1 »CHO  *  CH2  *EN0*EN2*AK.2»Q2/X»T ) 

DEPEND  (A.B.C.FO/X.Y.T) 

(5  HI  .2.1.1) 

EPTERM  -  0  (j  ni. 2.1.2) 


LBL  PARAM  IPIK.AAIl)).  IP  IE . AA (2 II »  I WR  .  AA ( 3 ) ) * 
1  (B.AAI5I).  ( C  *  AA  (6 )  )  •  (  FO*  AA  (71) 


! A  * AA ( 4 ! ) * 

(ft  m.E.i.s) 


PIK  l  2  (ft  m.2.i.4  («)) 

FMCEXPI ALAM )  *  (1.  +  2.  *  FZ  *FMCEXP  I  ALAM )  *  AMUO/XIS 

PIE  ■ 

FMCEXPI -ALAM)  *  t 1 .-2.*EZ*FMCEXP I ALAM ) *AMUO/X )  $ 

WR  a 

OM  -  0*B  +  EZ  *  AJ1  *  (5  UZ .2.1.4  (b)) 

A  « 

?.  *  EZ  *  FMCFXPI-ANU)  *  AJ1  *  (OM  -  OMR)  *  X  *  X  *  FMCSIN(Y)  **2+2. 

*  HO  ♦  EZ  *  ENO  +  (2.  *  H2  +  EZ  *  EN2I  *  0.5  *  (3.  «  FMCCOS(Y)  ** 

2  -  1.1  5 

R  = 

2.  *  FMCEXPI ALAM )  *  I  EMO/X  ♦  EM2/X  *  0.5  *  (3.  *  FMCCOS(Y)  **2  -  1.)) 

*  (1.  -  2.  *  EZ  *  FMCEXP(ALAM)  *  AMUO/X)  -  EZ*(CH0+CH2*  0.5* 
I3.*FMCC0SIY)**7  -  1.))  s 


c  = 

(7.  *  I V7  -  H2>  -  FZ  *  AK2 )  *  0.5  *  (3.  *FMCCOMY>  **  2  -  1.)  S> 


FO  = 

EZ  *  02  *  0.5  *  (3.  *  FMCCOSIY)  **?-!.)$ 


1  S  RICCI  TENSOR  COMPONFNT 

5.0E-I*EP*PIE*X**2.0*FMCSIN(Y) **2.0*FMCEXP  I -  ANU) *FMCDI K ( WR . ( T  *  1 ) . ( 
X.l)  )-2.5E-l*EP*PIF**2.0*X**2.0*FMCS!M  Y )  **2. 0*FKCEXP ( -ANU ) *FMiCDl  F 
IPIK.IT.l) ) *FMCD I  FI WR .  I  X.l ) ) S 


(ft  m.2.1.4  (c)) 


4  3  RICCI  TENSOR  COMPONENT 

-2.0*FP*PIE*X*FmCF.NIY)**2.0*FmCFXD(-ANU)*FMCDIF('a'R,  (X.l  )  )  + 
7.5F-1*EP«PIF*X**?.0*FmcsIN(Y)**?.0*FMCFXP(-ANU)  *F,'COIF(  aa-j.IX.  1  )  ) 
*FVC0IF(VR,(X,1 ) )-5.0F-l*FP«PIF*X**2.0*FMCSIN( Y)**7.0#FMCFXPI-AMU) 
*FvCOIc(WR» (X .? ) )-7.5F-l*FP*PIF**2.n*X**2.0*F"C5iM(Y)**7.0*FHCEXP( 
-A»CUI*FMCDIF  (PI<-  .(X.l)  )*F.'«COir  CWR.IX.I  )  )-5.nE-l  *EP*X**?.0J  F.'iCSI  NI  Y 
)**?.0*FvcrxP|-ANUI*FMCI)ir(PIF  ,  (X.l  I  MFVCOIFC.-'R,  (  X  ,  1  )  )  % 
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4  3  RICCI  TENSOR  COMPONENT 

-2.0*EP*PIE*X*FMCSIN(Y)**2.0*FMCEXP(-ANU)*FMCCIF(WR,  (X, 1))* 
2.5E-l*eP*PlE*X**2.0*FMCSIN<Y)**2.0*FMCEXP(-ANU)*PMC0IFI ANU,  (X,  III 
♦FMCD I F  I  WR  , I  X,  1  I l-5.0E-L*EP*PIE*X**2.0*FMCSIN(Y)**2.0*FMCEXPI-ANU) 
♦FMCDIFIWR, ( X  *  2  I ) -2. 5E-1 *fc P*P I E **2 . 0*X **2 . 0*FMC S IN  (  Y I **2.0*FMCEXPt 
-ANUI*FMCDIF(PIK, <  X,  1  I )*FMCDIF I WR, I  X , 1  I  ) -5. CE- 1*E P*X**2 . 0*FMC$ IN ( Y 
)**2.0*FMCEXPI-ANU)*FMC01F(PIE, <X,1) )*FMCDIFIWR,  ( X , 1 1 ) t 

THE  EXPRESSION  AFTER  SUBSTITUTION 


AMUO*EP*EZ*X*FMCSIN< Y I **2 . 0*FMCE XP I  - ANU  I  *FMC0 I F I ANU, (X, 1 ) >*FMCDIFI 
OMB, I  X, I ) )*5.0E-l-AMU0*£P*EZ*X*FMCSINt Y  )  * *2 . 0* FMCE XP I - ANU I *FMCD I F I 
OMB, IX, 2) )+AMU0*EP*EZ*FMCSIN(Y)**2.O*FMCEXP(-ANU)*FMCDlFIOMB,IX»ll 
)*I-3.49999998)fAMU0*EP*£Z**2.0*X*FMCSIN(Y)**2.0*FMCEXP(-ANUI* 
FMCDIFI AJ1 , ( X,l» ) *FMCO I F I  ANU ,  I  X , 1 ) )*I-5.0E-1 I  <-AMUO*EP*EZ**2.0*X* 
FMCSINI Y)**2.0*FMCEXP( -ANU )*FMCD IF  I AJl ,  (  X , 2  I  I + AMU0*EP*EZ**2. 0* 
FMCSINI Y)**2.0*FMCEXP(ALAM-ANU>*FHCDIF< AMUO,  I  X  ,  l  I  I  *F MCD IF  I OMu, | X , 1 
)  )*(-2.0M-AMU0*tP*EZ**2.0*FMCSINI Y I  * *2 . 0*FMCEXP I -ANU ) *FMCD  IF  I AJl ,  I 
X,l) )*3.4999999b*AMU0*fcP*EZ**T.0*FMCSIN(Y)**2.0*FMCEXPI AL AM-ANU ) * 
FMCDIFI AJl, ( X, 1 ) ) *FMCDIF< AMUO,  I  X , 1) I *2 . 0+ AMU0**2. 0*EP*EZ**2. 0*X**I 
-1.0)*FMCSIN(Y)**2.0*FMCEXP(ALAM-ANU)*FMCDIF(OMB,(X,  1 ) )*2.0*AMU0** 
2.0*EP*EZ**2.0*FMCS1N(Y)**2.0*FMCEXPI A L  AM-ANU )*FMCCIFIALAM, IX, II )* 
FMCUIFIOMB, I  X, I  I  I *1-3.0) *AMU0**2 .0*EP*EZ **3 . 0*X** I -1 .0)*FMCSIN( Y ) 
♦*2.0*FMCEXP( ALAM-ANU )*FMC0I FI AJl, (X, 1  )  I  ♦ I -2 .0  I +AMU0**2 .0*EP*EZ** 
3.0*X**(-1.0I*FMCSIN1Y)**2.0*FMCEXP{ AL AM*2 . O-ANU I *FMCU I F I AMUO, I  X, l 
I  I  *FMCO I F I OMB , (X, 1  I  I *2 . 0»AMU0**2 .0*EP*EZ**3.0*FMCSIN(Y)**2.0* 
FMCEXPI ALAM-ANU) ♦FMCDIFI AJl, IX,  1  I  I *FMC D  I F I AL AM , ( X ,  1 1  I *3 .0* AMUO** 
2.0*EP*EZ**4.0*X**(-1.0)*FMCSINIY)**2.0*FMCEXP(ALAM*2.0-ANU)* 
FMCDIFI AJl ,  ( X, 1  I )*FMCQ|F( AMUO,  I  X, 1 1 ) *  I -2.01 ♦ AMUO** 3. 0*EP*E Z**3.0*X 
♦*I-2.0)*FMCSIN( Y)**2.0*FMCEXP(  ALAM*2 . O-ANU ) *FMCD I F  I  UMB, IX, 1  I  1*1 
-2.0)fAMUU**3.0*EP*EZ**3.0*X**(-1.0)*FMCSINIY)**2.0*FMCEXP(ALAM* 
2.0-ANUI*FMC0IF( ALAM, (X,l I  I *FMCO I F I OMb ,  I  X , 1 1) *4.0  + AMU0**3. 0*EP*fcZ 
**4.0*X**(-2.0)*FMCSIN(Y)**2.0*FMCEXP(ALAM*2.0-ANU)*FMCDIFIAJ1,1X» 

1)  )*2.0*AMUD+*3.0*tP*EZ**4.0*X**l-1.0) *FMCSINI Y)**2.0*FMCEXP(ALAM* 
2.0-ANUI*FMC0IFI AJl, I  X , 1  I ) *FMCU I F I AL AM ,  I  X , 1  I  I  *  I -4 .0 ) +EP*EZ *X* 
FMCSINI Y)**2.0*FMCEXP< -ALAM-ANU >*FMCDl FI  AJl, (X,l I  I ♦  I -2 .0 ) ♦ EP*EZ*X* 
FMCSINI YI**2.0*FMCEXP(-ANU)*FMCDIF(AMU0,IX,1) )*FMCCIFIOMB, (X,l) >*< 
-5.0E-1)+EP*EZ*X**2.0*FMCSINIY)**2.0*FMCEXP(-ALAM-ANU)*FMCDIFIAJ1, 
IX, 1 1 )*FMCDl FI ALAM, I  X, 111  *2 . 5E- 1  *EP*fc  Z  *X**2 .0*FMCS  IN  I  Y  I  **2 . 0* 
FMCEXPI-ALAM-ANUI *FMCOIF(AJi,  (X, 1) )*FMCCIF(ANU, IX,  1  I  I *2. 5E-1 *EP*EZ 
*X**2.0*FMCSIN(Y)**2.0*F  MC EXP  I- ALAM- ANU  )  *FMCDIF(AJ1,  (  X  ,  2  I  )  *  I 
-5.0E-l)+EP*EZ**2.0*X*FMCSINI Y  I  **2 .0*FMCEXP I -ANU ) *FMCDI F I AJl, (X, 1) 
)*FMCDIF( AMUO, IX, 1  I  I *5 . OE- 1 *EP*X*FMCS I N I Y I **2. 0*FMCEXP ( -AL AM-ANU I  * 
FMCDIFI OMB, | X, 1 1 ) *2.0f£P*X**2.0*FMCSINI Y  )  **2 . 0*FMCEXP I -ALAM-ANU ) * 
FMCDIFI ALAM, (X,l))*FMCDIF(OMB,(X,l))*(-2.5E-l)+EP*X**2.0*FMCSIN(Y) 
**2.Q*FMCEXPI-AIAM-ANU)*FMCDIF< ANU, I  X, 1  ) ) *FMCD I F I OMB , (X,l I  )*( 
-2.5E-1 I ♦  EP*X**2.0*FMCSIN<  Y 1**2 .0*FMCEXP( -ALAM-ANU )*FMCCIF 10MB, I  X, 

2) 1*5. 0E-1J 


G8 


ft  HI. £.4 


1  3  RICCI  TENSOR  COMPONENT 

5.0E-1*EP*PIE*X**2.0*FMCSIN( Y ) ♦*2.0*FMCEXP ( -ANU ) *FMCDIF ( WR, I T ,  1 ).  ( 
X,l) )-2.5E-l*EP*PIE**2.0*X**2.0*FMCSINI Y ) **2.0*FMCEXP I -ANU > *FMCD IF 
IP  IK, ( T,ll I *FMCDI F ( WR  ,  (X,  1 ) ) $ 

the  expression  after  substitution 


-AMU0*EP*E2**2.0*X*FMCSIN( Y) **2 . 0*FMCEXP( -ANU ) ♦FMCDIFI AJl, ( T,1 ) , (X 
,11  >+AMuO*E?*EZ**2.0*FMCSlN(Y » **2. 0*FMCEXP( ALAM-ANU > *FhOOI FI AMUO, I 
T , 1 1 J*FMCOIF(OMB,IX,l> ) *  ( -2 . 0 ) ♦AMU0*EP*EZ**3.0*FMCS I N I V )4*2. 0* 
FMCEXPI ALAM-ANU)*FMCDIF< AJl, (X, 1 ) ) *FMCDIF ( AMUO, ( T , 1 ) ) *2.0 ►AMUO** 
2.0*EP*EZ**3.0*X**(-1.0)*FMCSIN(YI**2,0*FMCEXP(ALAM*2.0-ANU>* 
FMCOIFI AMUO, I T , 1 ) I ♦FMCOI F I  OMB, I  X , 1 1) *2. 0*AMU0**2. 0*EP*EZ**4. 0*X**( 
-1.0)*FMCSIN( Y»**2.0*FMCEXP( ALAM*2 . O-ANU) *FMCD I F I AJl, !X,1> )*FMCDIF 
(AMUO, <  T,1 ) > *1-2.0  H-EP*EZ*X*FMCS1N< Y I **2. 0*FMCEXP I -ANU > *FMCD I F I 
AMUO,  ( T,l i )*FMCDIF(UMB,IX,1I ) *5 . 0E-1*EP*EZ*X**2. 0*FMCS INI Y ) **2 . 0* 
FMCEXP(-ALAM-ANU)*FMCDIF(AJ1,(T»1J,(X,1))*5.0E-1*EP*EZ**2.0*X* 
FMCSINt Y)**2.0*FMCEXP(-ANUI*FMCDIFIAJ1, (X,l > ) *FMCC IF  I AMUO, (T,l) )♦( 
-5.0E-1 I » 
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IV.  A  BRIEF  GUIDE  TO  FORMAC  PROGRAMMING* 

A.  The  FORMAC  System 

In  writing  a  program  using  FORMAC  it  is  helpful  to  consider  the  FORMAC 
compiler  as  an  extension  of  the  facilities  of  FORTRAN  IV  which  allows  for  the 
manipulation  of  algebraic  symbols.  The  FORTRAN  IV  features,  such  as  integer 
and  floating  point  arithmetic,  type  statements,  data  statements,  and  logical 
operations,  all  exist  and  can  be  used  in  FORMAC  programs. 

When  a  FORMAC  program  is  submitted  for  compilation  the  machine  Imnediately 
produces  a  purely  FORTRAN  program,  and  the  FORMAC  commands  are  rendered  as 
comnent  statements. 


For  further  details  see  the  FORMAC  mnual  (IB4  1962). 
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B.  FORMAC  Executable  Statements 


5  rv. 


Because  algebraic  quantities  rather  than  numbers  are  usually  being 
manipulated  by  a  FORMAC  program,  input  data  are  considered  as  alphanumeric 
Information  and  are  read  into  the  machine  by  the  FORTRAN  IV  "A"  format.  A 
thirty-six  bit  computer  word  in  core  storage  can  contain  6  alphanumeric 
characters.  Thus  72  columns  on  an  I  EM  data  card  represent  12  computer  words. 
The  program  SUBSTITUTE  illustrates  how  an  algebraic  expression  up  to  50  cards 
in  length  can  be  read  into  memory. 

Once  an  expression,  say  ALGEXP,  has  been  entered,  it  must  be  converted 
to  a  form  that  can  be  manipulated  conveniently  by  the  compiler.  This  is  done 
by  the  FORMAC  executable  statement  ALGC0N 
K  =  0 

LETT  A  =  ALGC0N  ALGEXP,  K 

It  must  be  emphasized  that  all  variables  which  appear  in  expressions  that  are 
read  in  as  data  must  be  declared  in  the  ATOMIC  and  DEPEND  statements. 

When  a  quantity,  such  as  A  above,  has  been  put  into  proper  form  by  means 
of  an  ALGC^N  statement,  it  can  then  be  operated  on  by  other  FORMAC  executable 
statements.  The  following  executable  statements  are  available  in  FORMAC. 
(Note  that  some  of  these  statements  have  been  described  in  earlier  sections.) 

AUTS1M  Used  to  evaluate  elementary  functions  that  have 

constants  or  FORTRAN-type  expressions  as  arguments. 

BCDCjiN  To  be  described  below. 

CENSUS  Used  to  count  the  number  of  computer  words,  terms, 

or  factors  in  a  FORMAC  expression. 
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C0EFF  Used  to  obtain  the  coefficient,  in  an  algebraic 

expression,  of  an  atomic  variable  or  an  atomic  variable 
raised  to  a  specified  power. 

ERASE  Used  to  erase  a  FORMAC  quantity  in  core  storage. 

EVAL  Used  to  evaluate  a  FORMAC  expression  using  supplied 

values  of  the  variables.  This  results  in  a  numeric 
quantity. 

EXPAND  Used  to  simplify,  collect  terms,  and  remove  parentheses 

in  algebraic  expressions. 

FIND  Used  to  establish  a  function  relation  between  a  FORMAC 

expression  and  prescribed  atomic  variables. 

LET  Used  to  define  algebraic  expressions. 

MATCH  See  below. 

PART  See  FORMAC  manual. 

SUBST  See  below. 

jtoBER  Used  to  specify  the  sequence  in  which  factors  in  products, 

and  terms  In  sums  are  to  appear  In  expressions . 

These  cotoman>  .  are  described  in  more  detail  in  the  FORMAC  annual,  and 
most  of  them  were  used  in  the  ALBERT  programs. 

It  is  useful  to  note  that  most  operations  performed  on  numeric  quantities 
have  their  FORMAC  analogs.  For  example,  logical  operations  can  be  performed 
on  FORMAC  variables  Just  as  they  can  be  on  FORTRAN  variables:  One  algebraic 
expression  can  be  coopered  to  another  by  the  FORMAC  function  MATCH.  The 
sequence 

ATOMIC  X,  Y 
LOGICAL  Q 
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LETT  C  =  X**2  +  Y**2 
LET  Q  =  MATCH  ID,  B,  C 
IF(Q)  GO  TO  2 
CONTINUE 

would  compare  B  to  C;  and  if  they  were  equivalent,  it  would  eet  the  value  of 
q  to  the  FORTRAN  logical  TRUE  and  then  transfer  to  statement  number  2  for  its 
next  instruction. 

Once  a  desired  algebraic  quantity  has  been  constructed,  it  can  be  evalu¬ 
ated  numerically,  or  substitutions  can  be  made  that  will  convert  it  to  yet 
another  symbolic  quantity.  The  FOFMAC  operation 
LET  Z  a  SUBST  B,  LIST 

will  result  in  a  new  algebraic  quantity,  Z.  LIST  is  a  listing  of  the  symbols 
in  B  and  the  names  of  new  expressions  that  are  to  take  their  places.  The 
actual  formtion  of  LIST  is  done  by  the  FOHMAC  function,  PARAM.  For  example, 
consider  the  following: 

ATOMIC  X,  Y  ALPHA,  BETA 

LET  C  =  X**2  +  Y**2 
1ST  A  -  ALPHA**2 
1ST  B  «  BETA**2 
LIST  PARAM  (X,A),  (Y,B) 

LET  Z  =  SUBST  C,  LIST 

Z  now  has  the  value  ALPHA**4  +  BETA**4 .  Notice  that  the  statements  on  cards 
begin  in  column  7  but  the  word,  LIST,  is  placed  in  columns  1  ♦  4. 

Assuming  a  computation  has  been  completed,  the  information,  other  than 
numeric,  must  be  put  into  a  form  suitable  for  printing.  This  is  accomplished 
by  the  FORMAC  command  BCDC|6n. 


73 


§  IV. B 


DIMENSION  CARD  (l 2) 

ASS  m  0 

LET  AJB  =  BCDC0N  Z,  CARD, 12 

will  cause  alphanumeric  Information  to  be  stored  In  the  12  computer  words  of 
the  array  CARD.  This  can  In  turn  be  written  out  In  FORTRAN'S  "A"  fornmt. 
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What  has  teen  attempted  with  the  ALBERT  programs  Is  to  make  them  as 
simple  and  easy  to  use  as  possible.  As  mentioned  earlier,  FORMAC  on  the 
IBM  7094  suffers  from  a  severe  lack  of  core  storage.  The  program  EINSTEIN 
plus  system  overhead  took  In  the  neighborhood  of  5/7  of  the  32K  memory 
capacity  of  the  machine.  Reading  the  input  metrics  In  on  data  cards  would 
conserve  some  core  space;  but  for  clarity  we  chose  to  define  the  metric  by 
means  of  LET  statements  In  the  body  of  the  program.  We  have  circumvented 
the  care  problem  to  some  extent  by  simplification  of  the  metrics  themselves, 
by  punching  the  output  from  EINSTEIN  on  cards,  and  by  expanding  that  output 
in  terms  of  new  variables  using  the  SUBSTITUTION  program  (cf.  §  IH.E.4). 

FORMAC  provides  far  the  use  of  an  overlay  feature,  which  might  be 
another  alternative  for  solving  core  problems.  Yet  another,  is  to  purge 
from  the  system  portions  of  the  user's  program  after  they  have  been  executed 
(see  Clemens  and  Matzner  1967).  However,  the  only  truly  satisfactory  solu¬ 
tion  to  the  core  problem  will  be  to  go  to  a  computer  with  a  larger  memory  — 
e.g.  the  IBM  360  Berles. 
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